Maple Questions and Posts

These are Posts and Questions associated with the product, Maple


This may sound like a dumb question, but i'm seeking a procedure to do it better.



with(combinat, setpartition) :
P := [$2..5] :

Tours := setpartition(P);M:=nops(Tours)

[[[5], [2, 3, 4]], [[2], [5], [3, 4]], [[3], [5], [2, 4]], [[4], [5], [2, 3]], [[2], [3], [4], [5]], [[2, 3, 4, 5]], [[2, 5], [3, 4]], [[2], [3, 4, 5]], [[2, 4], [3, 5]], [[3], [2, 4, 5]], [[2, 3], [4, 5]], [[4], [2, 3, 5]], [[3], [4], [2, 5]], [[2], [4], [3, 5]], [[2], [3], [4, 5]]]





number of elements in each 'group'


2, 3, 3, 3, 4, 1, 2, 2, 2, 2, 2, 2, 3, 3, 3



i need to add 1 to each 'subgroup' : These are the first two:


[[[1, 5], [1, 2, 3, 4]], [[1, 2], [1, 5], [1, 3, 4]]]



I need to add 1 to each 'subgroup' in a more automatic way.





Please sir,

I am trying to solve this problem and is giving me the error code "Error, (in thaw) numeric exception: division by zero
". Anyone with useful informations please.

Thanking you in anticipations for your response. 


h := proc (x) options operator, arrow; 1+x*tan(theta)/delta+phi*sin(2*Pi*x) end proc:

A := M*(1+K):

theta[1] := sqrt((-B+sqrt(-4*A*C+B^2))/(2*A)):

xi := H[a]^2+(1+K)/K[p]:

L := xi[1]*cosh(theta[1]*h(x))+beta*xi[1]*theta[1]*sinh(theta[1]*h(x))*sinh(theta[2]*h(x))-xi[2]*cosh(theta[2]*h(x))-beta*xi[2]*theta[2]*sinh(theta[2]*h(x))*sinh(theta[1]*h(x)):

F := proc (x) options operator, arrow; (2*xi[1]/theta[1]-2*xi[2]/theta[2])*sinh(theta[1]*h(x))*sinh(theta[2]*h(x))-2*L*h(x) end proc:

w := S*(sinh(theta[1]*y)*sinh(theta[2]*h(x))-sinh(theta[2]*y)*sinh(theta[1]*h(x)))/(L*xi):

u := S*(xi[1]*cosh(theta[1]*y)*sinh(theta[2]*h(x))-xi[2]*cosh(theta[2]*y)*sinh(theta[1]*h(x))-L)/(L*xi)-1:

psi := S*(xi[1]*sinh(theta[1]*y)*sinh(theta[2]*h(x))/theta[1]-xi[2]*sinh(theta[2]*y)*sinh(theta[1]*h(x))/theta[2]-L*y)/(L*xi)-y:

q := int(u, y = -h(x) .. h(x)):


Error, recursive assignment


"Deltap:=(∫)[0]^(1)(∂P)/(∂x) ⅆx"



Dear all;

I have a function  C(S,t)  and I would like to show that this function is the solution of Partial differential equation ( denotes by PDES). I introduce all my function but when I try to simplify my PDES, I can not get zero as result. 

Many thanks for your help to simplify the PDES to get PDES=0. So, that C(S,t) is the solution of my PDEs.


I am trying to construct a random set of angles measured in degrees with pairs that are the equivalent in radians.  for example [0,0],[90,Pi/2],[-30,-Pi/6],etc...  but completely random.

so far I have the following:

a := rand(-90 .. 90);
seq([a(i), convert(a(i)*degrees, radians)], i = 1 .. 4);

the trouble is it appars it is randomizing a both times it is called in the sequence.  I would like to have it randomize for each iteration in the sequence.

thanks in advance.

I want to solve a pde equation:

equa1 := diff(u(x,y), x, x)-y(1+x) = 0;

# with codition:

con:=u(0,y) = 0, (D(u[x]))(0,y) = 0;

the anwer must be :    u(x,y)= y(x2/2  + x3/6)
How can i solve that with maple?

Please excuse my bad English

I am plotting a pair of lists of points using

pointplot(Listap, Listbp, symbol = point, symbolsize = 1, size = [1200, 1200])


How could is do this with plot so I can add colours? Along the lines Listap(i)^2+Listbp(i)^2 =R, R is in the range 0..1,then colour =R*256 or any other imaginative way of adding colour.








z := (m+I*n)/(p+I*q)



g := proc (z) options operator, arrow; (z-I)/(z+I) end proc;

proc (z) options operator, arrow; (z-I)/(z+I) end proc


bz := simplify(evalc(Im(z)));



a := simplify(evalc(Re(g(z))));



b := simplify(evalc(Im(g(z))));



"  r:=15;   Lista:=Vector();  Listb:=Vector();  j:=1;  for m from -r to r do   for n from -r to r do   for p  from -r to r do   for q from -r to r do  if p<>0 and q<>0 and m^2-2 m q+n^2+2 n p+p^2+q^2<>0 and bz>=0 then  Lista(j):=a; Listb(j):=b;  j:=j+1;  end if;  end do:  end do;  end do;  end do:  j; :"








pointplot(Lista, Listb, symbol = plottools:-point, symbolsize = 1, size = [1200, 1200])








On a plot, how do I change the size of the probe info text box ? I cannot read the very small text displayed when the probe info is active. My laptop screen size is only 13 inches.

Thank you

So I have the following:

And I want to create an animation similar to something like:

Not necesarily 3D. Just the interative slider part. 

I want to create a slider that adjusts the value of the exponent k1 (from the red curve) from 1 to say, 1000. 

I also want it to show the corresponding value of  G[1] on the application.

I'm sure its pretty simple, it's just that I have no idea where to start.


Hello I have the following small piece of code.


which exports sin(theta)=.25, how do I force the 0 to display.



why does the following fail:

sol:=solve([cot(x)=2,x>=0,x<2*Pi],x,AllSolutions, Explicit);

How do I get a system solve to dispaly an ordered triple and not {x=ans, y=ans, z=ans}

for instance

sys := { z = 4, x+y= 10, x-y = 5 };
L:=solve(sys, {x,y,z} );


{x=7.5, y=2.5, z=4}

I want 


I have a piecewise that exists on a large intervals, say [-5..15], however, I know need the function to represent [1..9].

There are two approaches I know of: using convert(f,pwlist) or using op(n,f), see attachments.    and  PiecewiseTruncate.pdf

However, these two methods are not elegant enough, I have to reassemble the pieces in the new function.

How to make the solution neat?


I really want to use if condition inside of eval and sum, example as below:

sum(eval(y=x^(j),(if j=1 then x=2 else x=3 fi)),j=1..2)

The reason is: the value of x to be evaluated depends on the value of j which differs inside of the sum (in the context of B-spline functions).

Any lights? Thanks,

i tried to solve a nonlinear ode with numerical method but maple can't solve it and this error occur:

Error, (in dsolve/numeric/bvp) initial Newton iteration is not converging

my maple codes are attached below:

can any help me?

Hi everybody,

I face the following problem :
I develop a Maple simulation code with its own maplets-based interface, and I would like to draw up an "execution report" that presents many informations in a smart and synthetic way :

  1. a summary of the data the user provided through the interface (numerical values, files, ...)
  2. the results of some preliminary operations (ie. the construction of the equations to solve)
  3. some plots of the numerical solution (typically of an ode system)

This execution report should mix texts, plots, charts, ...
At least, I would like the report to be written in LaTeX.

I use to use the software R which provides some facilities to build reports of this kind  (packages knitr or Sweave for those who would know) , and I am wondering if there exists something comparable in Maple ?
Could the DocumentTools package (I know by name only) help me ?
More generally, if you were faced with such a problem, how would you deal with it?

As always any suggestion or advice will be greatly appreciated

Thanks in advance

I need to formulate the follow relationship in proper math symbols:

Differential in A (last days value minus today's) has a tendency to reach the Differential in B (today's)

I though this could be expressed with

AΔ -> BΔ

BUt I guess there are more elegant and mathematically correct ways to do this in Maple?

thank you!

1 2 3 4 5 6 7 Last Page 2 of 1253