i got some trouble when i tried to build large matrix. in my case, notification error out of bound appear when looping stop at 9 from 24 repeatation. 

and this is my looping command:

the result of the script was:

now i feel so desperate so finish my final project because the error, please help me

Hello,
I try to visualize this formula in a 3D graphic.
However, I get again and again this error message displayed and unfortunately knows no more solution.
I am a total beginner with Maple, I hope here can help me.

plot3d(100*(0,15+0,035*x)+100*(0,15+0,035*9), 100*(35*9))+0,05*(100-(100*(0,15+0,035*x)))*(0,15+9*y), x = 2 .. 9, y = 0, 1 .. 0, 45, axes = boxed);

Error, (in plot3d) unexpected options: [5*(100+(0, -1500, -3500*x))*(0, 15+9*y), x = 2 .. 9, y = 0, 1 .. 0, 45]

Hi everyone

I am stuck in my code. I have a multivariate polynomial and I am trying to define an order on the variables such that I can write my polynomials in a desired normal form. I tried with the Groebner package but that's quite different from what I want.

I want to define the order on the variables as u[i]=v[i] and u[i]<u[i+1].

Lets say my polynomial is f then,

input f= u[1]^2+u[2]+v[3]

output f = v[3]+u[2]+u[1]^2.

If there is a clash between u and v then, either of them can come first, for example,

input f= u[1]^2+u[2]+v[3]+u[3]^2

 output f=v[3]+u[3]^2+u[2]+u[1]^2 or f=u[3]^2+v[3]+u[2]+u[1]^2

I hope my question is clear. Thank you for helping me out. Your time is much appreciated :)
EDIT 1: Comment towards the end of the answer is a solution for this question.

Dear all,

I have the following problem: Maple does not simplify the denominator in the following example:

which gives

However, the result should be B. If only the denomiator is expanded it works: 

gives

which equals the nominator except for the B...

How can I use simplify in order to yield the desired result? 

Thanks a lot!

Hello maple users. 

I am learning to use this incredible software, but I am not getting to obtain all binary 4x4 matrices of rank 3 and 4?

Anyone knows how to do that ?

Thanks
My best regards!

Gustavo

Hello

Dear I want to solve the system of ODEs in attached file for different values of m1 i.e., 3,4,5 how I will use the value of m1 in dsolve. I am waiting your positive response thanks in advance.

Help.mw

Hi, My Error is this:

DEtools; DEplot(MODEL, VARS, DOMAIN, RANGE, [IC1, IC2], stepsize = .1, arrows = THIN, linecolor = COLORS);
%;
Error, (in DEtools/DEplot) invalid range for independent variablehow can i remove error?

Hi,

Why pdsolve is not correctly solving the following

where I get    

but it is correctly solving the next one

where I get ???

I see that it is due to the fact that the variable x is multiplied by a constant, but why is not maple able to manage that?

Thanks for your help,

Javier

Hi Mapleprimes,

I wrote some quick maple code -

for a from 1 to 10 do
if isprime(a)=false then print(ifactor(a)) else print('1',a) end do;
end do

I want to use the factors command to reproduce oeis.org/A027750/

Can anyone help me with this?

Regards,

Matt

Good day!

As part of an exercise I've calculated the length of a hypotrochoid numerically. To check my result I repeated the calculation in Maple, but received a different result. When double checking using WolframAlpha I got the same result as with my numerics. Maybe someone of you can tell me where I made a mistake.

Thanks in advance.
Sören

Link to WolframAlpha calculation: http://www.wolframalpha.com/input/?i=x%28t%29+%3D+%281-0.6%29+cos%28t%29+%2B+0.8+cos+%28+%281-0.6%29%2F0.6+*+t%29,+y%28t%29+%3D+%281-0.6%29+sin%28t%29+-+0.8+sin+%28+%281-0.6%29%2F0.6+*+t%29+from+t%3D0+to+6*pi

Download hypotrochoid.mw

I have a complex function for both the electric and magnetic fields for 2 laser pluses colliding they are only 2 dimensions z,t;

Where a,,  are real constants and Z is the initial offset for the pluses,

I would like to plot them together with there real and imaginary parts on two axes and then extended along the z direction, if possible I would like to animate them. To hopefully get a moving plot like this except it's a laser pluse not a continuous wave,

I want write function get value a1, a2, ..., an to b1, b2, ..., bn.

Example get 1, 3, -9, 5 to a, b, c, d. This is my code

myfunc := proc(oldVars,newVars)

    local i;

    if nops(oldVars) = nops(newVars) then

        for i from 1 to nops(newVars) do

            newVars[i] := oldVars[i];

        end;

    fi;

end proc:

But not work! Can someone help me. Thank you very much

Hello,

In the maple.ini, i defined some folders where i have put my packages.
It seems that, when I made some altercations in the maple.ini, they are not taken into account immediately.
Do you have ideas about when the update of the maple.ini is taken into account ?
And how can I force the update ? Should I restart my computer so that the update of the file is taken into account ?
Thanks a lot for your feedback.

Hey guys,

I have the following occurence:

ii_inf:=x^(2-s)*(x^(-s)*GAMMA(3-s)*GAMMA(2-2*s)/(GAMMA(2-s)*GAMMA(3-2*s))+x^(s-2)*GAMMA(3-s)*GAMMA(2*s-2)/GAMMA(s))/(2-s)+(1/2)*(2*s*x-x+1)*(x+1)/(((x+1)^s)^2*(2*s^2-3*s+1))+x^(1-s)*(x^(-s)*GAMMA(2-s)*GAMMA(-2*s+1)/(GAMMA(1-s)*GAMMA(2-2*s))+x^(-1+s)*GAMMA(2-s)*GAMMA(2*s-1)/GAMMA(s))/(1-s)+(x+1)/((2*s-1)*((x+1)^s)^2)+x/((-1+s)*x^s)-(x+1)/((-1+s)*(x+1)^s);

ii_inf=simplify(ii_inf);

asympt(ii_inf,x,3);

Multiseries:-asympt(ii_inf,x,1);
gives different results...the last one however seems to be the correct one...

What is happening here?

I want to solve the problem described below. I tried using two methods as shown below, each method has been runing for days without solving it. I will really appreciate your help.

Det1 := (1/256)*(Aiso*(c+t)^2*(a^2+b^2)*(mu-1)*Pi^2-4*a^2*b^2*c*Gc)*(16*Aiso^2*Do^2*(c+t)^4*(a^2+b^2)^6*Pi^12+10*(a^2+b^2)^5*((c+t)^2*Aiso+4*Do)*Gc*(c+t)^2*c*a^2*Do*b^2*Aiso*Pi^10+(a^2+b^2)^4*((c+t)^2*Aiso+4*Do)^2*Gc^2*c^2*a^4*b^4*Pi^8-(1024/81)*a^6*Aiso^2*b^6*Tcr^2*(c+t)^4*(a^2+b^2)^2*Pi^4-(2560/81)*a^8*Aiso*b^8*c*Gc*Tcr^2*(c+t)^2*(a^2+b^2)*Pi^2-(1024/81)*a^10*b^10*c^2*Gc^2*Tcr^2)*(Aiso*(c+t)^2*(a^2+4*b^2)*(mu-1)*Pi^2-4*a^2*b^2*c*Gc)*(Aiso*(c+t)^2*(a^2+b^2)*(mu-1)*Pi^2-a^2*b^2*c*Gc)*(16*(c+t)^4*(a^2+(1/4)*b^2)^3*Do^2*(a^2+4*b^2)^3*Aiso^2*Pi^12+(10*(a^2+b^2))*((c+t)^2*Aiso+4*Do)*Gc*(c+t)^2*(a^2+(1/4)*b^2)^2*c*a^2*Do*(a^2+4*b^2)^2*b^2*Aiso*Pi^10+((c+t)^2*Aiso+4*Do)^2*Gc^2*(a^2+(1/4)*b^2)^2*c^2*a^4*(a^2+4*b^2)^2*b^4*Pi^8-(1024/81)*(c+t)^4*(a^2+(1/4)*b^2)*Tcr^2*a^6*(a^2+4*b^2)*b^6*Aiso^2*Pi^4-(2560/81)*a^8*Aiso*b^8*c*Gc*Tcr^2*(c+t)^2*(a^2+b^2)*Pi^2-(1024/81)*a^10*b^10*c^2*Gc^2*Tcr^2)*((mu-1)*(c+t)^2*(a^2+(1/4)*b^2)*Aiso*Pi^2-a^2*b^2*c*Gc)/(b^20*a^20*(c+t)^16) = 0;

# method 1;
EQN := RootOf(Det1, Tcr);

EQN_2 := allvalues(EQN);

# method 2;

EQN := solve(Det1, Tcr);