Georgios has slipped the 2 into the wrong place; you want

polarplot([2*sin(theta),theta,theta=0..2*Pi]);

&* can be used as a non-commutative operator. You can have this and similar operators and assign them various properties. See ?neutral and ?define

you have z31x*31 where you mean z31*x31 and z32x*21 where you mean z32*x21.

When I correct those, I get 346 solutions.

Cheers,

David.

I'm not a mathtype user, but the web site suggests you can input in latex. Since Maple can generate latex with the latex command, in principle this should work. Perhaps someone else has tried this?

I agree that it is not modifiable; I do not know how to get that.

Frequently at conferences you have to submit your powerpoint presentation to the conference organizers for use on their computers. The I avoid Mathtype since it is not loaded on most conference machines, and then your output is "font nightmare". I have even been to conferences where even the simple Word equation editor was not installed. So I avoid modifiable in favour of the always displayable. This is especially the case since one can use slide years later, forgetting that they were created with some special software in mind.

If you just want a nice equation in powerpoint, as opposed to getting the formula into mathtype, then you can export the maple worksheet as *.rtf, open it in Word, and then copy the equation from Word to Powerpoint

(at least in Windows with Maple 10, this works). Exporting as html and then copying the gif files give much poorer output quality.

Cheers,

David.

See the resultant command. According to the documentation, this uses the subresultant algorithm for low degree polynomials. So you could steal the code in `resultant/subresultant/univar`. I'm not familiar with this algorithm, but perhaps printing out v gives you the sequence you want?

subres:=proc(p,q,x) option`Copyright (c) 1992 by the University of Waterloo. All rights reserved.`;
 local c,d,u,v,r,g,h,du,dv;
 userinfo(2,resultant,`using the subresultant algorithm`);
 u:=p;
 du:=degree(u,x);
 v:=q;
 dv:=degree(v,x);
 c:=1;
 g:=1;
 h:=1;
 userinfo(3,resultant,`degree of inputs is`,du,dv);
 while 0<dv do
   d:=du - dv;
   c:=c*(-1)^(du*dv);
   r:=prem(u,v,x);
   u:=v;
   v:=r;
   print(v);
   du:=dv;
   dv:=degree(r,x);
   userinfo(3,resultant,`degree of pseudo remainder = `||dv);
   divide(v,g*h^d,'v');
   g:=coeff(u,x,du);
   if d=1 then
     h:=g
   else
     h:=iquo(g^d,h^(d - 1))
   end if
 end do;
 c*iquo(v^du,h^(du - 1))
end proc:

and then

subres(3*x^4+x^2+5,x^3+3*x-7,x);

gives the output:

5-8*x^2+21*x

-343+673*x

6171072

96423

If you don't care whether or not your matrix is symmetric or Hermitian, then you can use any nonsingular matrix, R to transform your diagonal matrix of eigenvalues.

R.D.R^(-1);

So you can generate R using RandomMatrix, but you should check that it is nonsingular. This is still not the most general way, but is useful for most purposes (the matrices generated are diagonalizable, with diagonal Jordon forms).

I don't have v11, so I can't see everything on your worksheet, but the following idea should work. If your global "constant" is r and the functions defining your parametric plot are something like

alpha:=x->r*sin(x);
beta:=x->r*cos(x);

then you can plot curves with several r values by

plots[display](seq(plot([alpha,beta,0..2*Pi]),r=[1,3,5]));

If you want more individual control over the plots, e.g. different colors, then you can store the plots in variables and then use display to put them together .

David.

Eval(y,x=0);

gives

Eval(y,x=0);

(though the Maple tag seems not to display it well in the post). This is the mathematical equivalent of your statement in words.

interface(displayprecision=3);

should do the trick.

Try phaseportrait or DEplot in the DEtools package.

There is a missing line in your output - a list of integers which can be converted to known mersenne primes. My output has the following line after the first "end if":

w := [2, 3, 5, 7, 13, 17, 19, 31, 61, 89, 107, 127, 521, 607, 1279, 2203, 2281, 3217, 4253, 4423, 9689, 9941, 11213, 19937, 21701, 23209, 44497, 86243, 110503, 132049, 216091, 756839, 859433, 1257787, 1398269, 2976221, 3021377, 6972593, 13466917, 20996011, 24036583, 25964951, NULL];

So the routine has no smarts... Try showstat if you can't see it on another try.

If you think about eval([x,y],{x=0,y=-1}) giving [0,-1], then you can work this up with map2 to give what you want:

q:={x = 0, y = -1}, {x = 3/4, y = -7/16}:

map2(eval,[x,y],[q])[];

gives

[0,1],[3/4,-7/16]

as required. This works for any number of solutions returned by solve.