Maple 2015 Questions and Posts

These are Posts and Questions associated with the product, Maple 2015

I'm doing an interactive component and i need to plot a function and some points (which can change).
My question is if is it possible display the graphic with a value of automatic Plot Ranges.
For example the graphic would started with the first point and finisched with the last, like in this image (but i did it manually).

I want a sort of autozoom.

I am looking for persons with experience in setting up maple for the blind. The setting is a 21-year old blind
 student who wants to do math. He uses Mac OS/X and the built-in VoiceOver screen reader, plus refreshable Braille display (one-line), a Braille printer, and Nemeth, the Braille extension for math symbols.

Initially, we are working from no experience whatsoever, trying ideas, mostly failure. Getting some advice from persons who are blind and using maple would be wonderful. If anyone knows of an instance, a person who is blind and regularly uses maple, then we would both would be delighted to receive advice on what to install, how to set it up, and how to use maple on a daily basis.

One trouble for the VoiceOver screen reader is verbosity of maple output. Any methods to reduce the amount of output would be welcome. To uderstand the issue, look at the output of these two commands when interface(prettyprint=0) has been set to help out the VoiceOver screen reader:

matrix([[1,2],6,7]]); Matrix([[1,2],[6,7]]);

The first outputs the same text, which VoiceOver will read accurately. The second dumps out two lines of extraneous and confusing information, only 20% of which is useful. The verbosity is not restricted to matrix packages linalg or LinearAlgebra. It is a problem wtih nearly all packages.

-Grant Gustafson, Salt Lake City, Univ of Utah Math Dept



I'm doing this interactive component and for coordinates of points I use one TextArea for each coordinate.

a:= Do(%CB);

x1:=Do (%x1);
y1:=Do (%y1);
x2:=Do (%x2);
y2:=Do (%y2);
x3:=Do (%x3);
y3:=Do (%y3);
x4:=Do (%x4);
y4:=Do (%y4);
x5:=Do (%x5);
y5:=Do (%y5);

#x6:=Do (%x6);
#y6:=Do (%y6);
#x7:=Do (%x7);
#y7:=Do (%y7);
#x8:=Do (%x8);
#y8:=Do (%y8);
#x9:=Do (%x9);
#y9:=Do (%y9);
#x10:=Do (%x10);
#y10:=Do (%y10);





if a=2 then F2:=y=LinearFit([1, x], vx2, vy2, x);
end if:

if a=3 then F3:=evalf[2](LinearFit([1, x], vx3, vy3, x));
end if:

if a=4 then F4:=LinearFit([1, x], vx4, vy4, x);
end if:

if a=5 then F5:=LinearFit([1, x], vx5, vy5, x);
end if:

I was browsing on maple and i sow this:

I would use a table to insert coordinates of point like this image, but i don't know how can I do.

Can you help me plese?

Thanks in advance

Dear friends:

I am facing two problems

1. one is to get solution of the below system of ODE for L=100 (highlited as red) and

2. the other is I want the graph in the form of solid line not poit, asterisk etc.


restart; epsilon := .1; Pr := 1; beta := .1; Sc := 1; S := 1; L := 20;
for i from -L while i <= L do;
a[i] := 1.0*i/L;
end do;
for i2 from -L while i2 <= L do;

fw := a[i2]; 

Eq1[i2] := eval(diff(F(eta), eta, eta, eta)+F(eta)*(diff(F(eta), eta, eta))-(diff(F(eta), eta))^2+S*(epsilon-(diff(F(eta), eta)))+epsilon^2);
Eq2[i2] := eval((diff(G(eta), eta, eta))/Pr-G(eta)*(diff(F(eta), eta))+F(eta)*(diff(G(eta), eta))); 
Eq3[i2] := eval(diff(H(eta), eta, eta)+Sc*(F(eta)*(diff(H(eta), eta))-beta*H(eta)));
IC[i2] := F(0) = a[i2], (D(F))(0) = 1, (D(F))(L) = epsilon, G(0) = 1, G(L) = 0, H(0) = 1, H(L) = 0;
dsys1[i2] := {Eq1[i2], Eq2[i2], Eq3[i2], IC[i2]};
dsol1[i2] := dsolve(dsys1[i2], numeric, output = listprocedure, range = 0 .. L);
dsol1x[i2] := subs(dsol1[i2], diff(F(eta), eta, eta));
dsol1y[i2] := subs(dsol1[i2], G(eta));
dsol1z[i2] := subs(dsol1[i2], H(eta)) end do;

for j from -L while j <= L do; 
g[j] := eval(-dsol1x[j](0)) end do;

g6 := pointplot({seq([n/L, g[n]], n = -L .. L)}, symbol = asterisk, symbolsize = 15, color = red);

Please see the problem and correct as soon as possible. I am waiting your positive respone.

Muhammad Usman

School of Mathematical Sciences 
Peking University, Beijing, China


How to pdeplot Curl(A(x,y,z))=0 in maple

it has A1(x,y,z) , A2(x,y,z) and A3(x,y,z)


How can I produce the same output with something like alphadeg := alpharad * 180/Pi?

T := proc (p, q, R) alpharad = evalf[5](arccos((1/2)*(p^2+R^2-q^2)/(p*R))), alphadeg = 180*evalf[5](arccos((1/2)*(p^2+R^2-q^2)/(p*R)))/Pi end proc;
T(3, 4, 3.75);
           alpharad = 1.2515, alphadeg = 71.70566805

When I try alphadeg = alpharad * 180/Pi I get
Error, illegal use of an object as a name

Thanks! Les

Is there a method to relate groebner bases with monomials ideals

How to plot magnetic field of maxwell equations like vector field

I am debugging a Maple program (and I cannot contact its author). Allegedly, the program ran under some old version of Maple, but it does not run under my Maple 2015 or Maple 9. I am trying to make it run under Maple 2015. I found a place where the program does not do what it is expected to do, but don't know how to make it do what I want:-)

I inserted a DEBUG statement at some place in the program, and the debugger window opens when the program comes to that statement. I execute  the following statement in the debugger:  traperror(Groebner:-NormalForm(g,AA,matermorder)); if I get some benign result, such as 0, I push "continue" in the debugger, and the program runs until it gets to the DEBUG statement again. I repeat these two operations several times until the traperror statement gives the following: polynomials must be members of the algebra" (my comment: only one quotation mark).

So I would like to amend the program to process the exception differently from how it is done now, but I don't know how to get access to the result of the traperror statement. For example, when I am trying to execute the following operators in the debugger: adebug:= traperror(Groebner:-NormalForm(g,AA,matermorder));adebug; , I only get the following diagnostics: Warning, extra characters at end of parsed string  So I don't understand what adebug contains and how I can introduce some logic based on the contents of adebug to handle the exception differently.

I suspect I am missing something simple due to my lack of experience with Maple. I tried to put the result of traperror into a parse statement, but did not get any satisfactory results.


Thank you 

I have to calculate a sort of trend line, but I don't know how it works with maple.
I tried with Excel.

Do you think the method is right? How can I replicate it on maple?

Thanks in advance

Dear friends! I am facing problem to solve the below system of ODEs numerically please find the mistake and correct it.

alpha := -1; R := 2; m := 2; Pr := 7; Le := 1.25; Nt := .2; Nb := .2; g := .5; K1 := .1; Q := .5

Eq1 := eta^3*(diff(F(eta), eta, eta, eta, eta))+alpha*(eta^4*(diff(F(eta), eta, eta, eta))+eta^3*(diff(F(eta), eta, eta))-eta^2*F(eta))-2*eta^2*(diff(F(eta), eta, eta, eta))+3*eta*(diff(F(eta), eta, eta))-3*(diff(F(eta), eta))+eta*R*(diff(F(eta), eta))^2-3*eta*R*F(eta)*(diff(F(eta), eta, eta))+3*R*F(eta)*(diff(F(eta), eta))+3*eta^2*R*F(eta)*(diff(F(eta), eta, eta, eta))-eta^2*(diff(F(eta), eta))*(diff(F(eta), eta, eta))-M^2*(eta^3*(diff(F(eta), eta, eta))-eta^2*(diff(F(eta), eta))); Eq2 := eta*(diff(G(eta), eta, eta))+alpha*Pr*eta^2*(diff(G(eta), eta))+R*Pr*F(eta)*(diff(G(eta), eta))+Nb*eta*(diff(G(eta), eta))*(diff(H(eta), eta))+Nt*eta*(diff(G(eta), eta))^2+diff(G(eta), eta)+Q*Pr*eta*G(eta) = 0; Eq3 := eta*(diff(H(eta), eta, eta))+alpha*Le*Pr*eta^2*(diff(H(eta), eta))+R*Le*Pr*F(eta)*(diff(H(eta), eta))+Nt*eta*(diff(G(eta), eta, eta))/Nb+Nt*(diff(G(eta), eta))/Nb+diff(H(eta), eta)-g*Le*Pr*eta*H(eta)-Le*Pr*K1*eta = 0;

IC1 := F(0) = 0, F(1) = 1, (D(F))(0) = 0, (D(F))(1) = 0, (D(G))(0) = 0, G(1) = 1, (D(H))(0) = 0, H(1) = lambda; dsys1 := {Eq1, Eq2, Eq3, IC1}; dsol1 := dsolve(dsys1, numeric, continuation = lambda, range = 0 .. 1);

dsol1x := subs(dsol1, F(eta));

dsol1y := subs(dsol1, G(eta)); dsol1z := subs(dsol1, H(eta));

With my best regards and sincerely.

Muhammad Usman

School of Mathematical Sciences 
Peking University, Beijing, China

Hi everybody, 

The sequence 

keeps printing the result on the screen

How is it possible to redirected the output of showstat(Myproc) into a file ?

Thanks in advance

how to convert a nested for loop to iterative version with stack

#my testing for wildcard to one
#after testing, it loop a very long time and not stop
ppp := [[0,0,0,x],[0,0,1,0],[0,1,0,0],[1,0,0,0]]:
check1 := [seq(0,ii=1..nops(ppp))];
ttt1 := [[0,0,0,1],[0,0,1,0],[0,1,0,0],[1,0,0,0]]:
mmmeaght1 := [seq(0,ii=1..nops(ppp[1]))]:
bbb1 := [seq(0,ii=1..nops(ppp[1]))]:
emap := [(xx) -> if [xx < 0 assuming x > 0] then 0 else 1 end if, (xx) -> evalf(1/(1+exp(xx)))]:
MM(ppp, ttt1, mmmeaght1, bbb1, check1, emap);

when test wildcard variable for input, would like to assume x > 0 then

i try assuming x > 0 , got error


I want to calculate the intersection between three circles.
I know that in this case i can calculate intersection of only the first and second equation, but I need this for a interactive component.

The command "intersection"[GEOMETRY] work only with 2 circles.

I did this but it doesn't work.


3 4 5 6 7 8 9 Last Page 5 of 44