Maple Questions and Posts

These are Posts and Questions associated with the product, Maple


restart; _local(gamma); _local(I); _local(pi)



Warning, The imaginary unit, I, has been renamed _I


# Set up numerical values for all problem parameters
  params:=[       gamma=0.142,        tau=0.112,      mu[1]=0.4e-3,
                  beta[1]=0.081,  b[h]=10,  psi=0.011, phi=0.05
            epsilon=0.2e-2,     rho=0.5e-1, beta[2]=0.092, beta[o]=0.034
                    q=.5,      eta=0.213,   M[h]=100
              delta=0.021,       alpha=0.57e-1,   p=.2,   beta[k]=0.025
               omega=0.056,      mu[c]=0.0019,  (mu)^(2)[1]=0.091

Error, missing operator or `;`


# Define main function
# R := sqrt((omega+mu[1]+eta)*mu[1]*mu[c]*psi*beta[o]*M[h]^2*(delta+mu[1]+eta+phi[1])*(`ε`+mu[1])*omega*beta[k]*   (1/((omega+mu[1]+eta)*(`ε`+mu[1])*mu^2)))

# Compute "all" derivatives and evaluate numerically.








  Y'a t-il des suggestion de résolution de ce qui suit? merci.

a) On pose F(-1,0) et F`(1,0) les foyers de l'ellipse E, on veut déterminer son équation implicite f(x,y)=0, 
    des points A(x,y) sur l’ellipse en utilisant la description bifocale et sans utiliser Maple.
  b) Utilisez la commande implicitplot de Maple et le résultat de a) pour tracer plusieurs ellipses sur
    le même graphe en variant la valeur de la constante C,(la somme des distances aux foyers F, F’ est
une constante qui est C).


  a) We put F (-1,0) and F` (1,0) the foci of the ellipse E, we want to determine its implicit equation f (x, y) = 0,
     A (x, y) points on the ellipse using the bifocal description and without using Maple.
  b) Use the command implicitplot of Maple and the result of a) to draw several ellipses on
     the same graph by varying the value of the constant C, ( FA +F'A = C ). 

     Give at the end a Maple procedure that generates this ellipse.

Hello there,

I'm trying to write something like: dsolve=(myOde,numeric,method=rkf45,abserr=1e-7,relerr=1e( -(7+x) ) ),

where x is a variable that updates every loop of my for loop (say for example x = 2). I'm basically trying to improve my accuracy each iteration of the loop.

Anyways, Maple doesn't like how there is a variable in the expression for relerr - it returns the command just as I typed it, not, for example, as

dsolve=(myOde,numeric,method=rkf45,abserr=1e-7,relerr=1e( -(9) ) ).

How can I achieve my goal?


with(plottools):F := proc (N) local a, b, L; L := NULL; L := sort([op({seq(seq(a/b, a = 0 .. b), b = 1 .. N)})]); return L end proc; F(6);
            [   1  1  1  1  2  1  3  2  3  4  5   ]
            [0, -, -, -, -, -, -, -, -, -, -, -, 1]
            [   6  5  4  3  5  2  5  3  4  5  6   ]
Ford6 := proc (i) local d, k, n, r; k := i; n := numer(F(6)[k]); d := denom(F(6)[k]); r := (1/2)/d^2; return [n/d, r], r end proc; nops(F(6));
for i to 13 do C || i := Ford6(i) end do;

display(circle(C1), circle(C2), circle(C3), circle(C4), circle(C5), circle(C6), circle(C7), circle(C8), circle(C9), circle(C10), circle(C11), circle(C12), circle(C13), axes = normal, scaling = constrained, color = blue, size = [800, 800]);

I need an example programme please help me i m beginner


I've used maple to solve a bunch of differential equations symbolically and the equations are pretty massive and so I don't particularly want to type by hand again on MS Word so if someone can tell me how to copy to word so it appears in an editable equation format I'll be very grateful.


Kind Regards



I am new to Maple.

As a student, how do I install Maple on my PC? 

I need it to solve PDEs and ODEs.

Many thanks for your kind reply.


Hello there,

I have created a MAPLE document in a slideshow format.

Is there a way to print a copy of my presentation with header and footer included and in a format that fits on letter format paper or pdf.

Thank you for your help.


I have 2 questions here

1st I want to extend this set Ug3 := {1, 19, 23, 29} by repeatively subtracting 6 fron its elements. first pass should give

Ug3 := (Ug3-~ 6) union Ug3  Ug3 := {-5, 1, 13, 17, 19, 23, 29}

Now I dont want to add any elenments less than 6 so -5 need to go

Second pass produces

Ug3 := (Ug3-~ 6) union Ug3                Ug3 := {-11, -5, 1, 7, 11, 13, 17, 19, 23, 29}

Obiviously no elements less than 6 to be added on second pass either.

At this point there is no further elements that could be added to the set.

I am looking for an efficient method. I can use For loops and if then checks. The sets grow in length quickly. Repetition lengths are 6, 30 210, 2310,... I use sets because the numbers get ordered automaticially. Are sets efficient or should I use some form of table?

I understand that Maple 2018 is now able to solve 3 independent variable PDE & BC problems in bounded domains through separation of variables by product and eigenfunction expansion.

My solution domain is (x,y,t), (i.e three independent variables) but I would like to need to use numerical integration.  Are there plans to make numerical integration for PDEs with three independent availables, and if so when is that facility likely to be available?




restart: with(plots):
bisect:=proc(f, intbeg::numeric, intend::numeric, ss::float, del::float)
   local Num_points,limit1,approx,i,x1,x2,a,b,flag1,pet,A,B,NewApprox,CheckIt,newDelta,
         interval_begin,interval_end,step_size,delta, y_approx, y_A, y_New;


for i from interval_begin by step_size to interval_end do
   if x1=0 then x1:=-10^(-5)fi:
   if x2=0 then x2:=-10^(-5)fi: 
   a := evalf( f(x1) ):
   b := evalf( f(x2) ):
    if a=0 then a:=-10^(-5)fi:
    if b=0 then b:=-10^(-5)fi:

   if evalf(a*b) <= 0 then 
     approx := (x1+x2)/2:
     y_approx := evalf( f(approx) ):
     if abs(y_approx) < delta then 
        print([approx, y_approx]):
        while pet=0 do
           y_A := evalf( f(A) );
           y_New := evalf( f(NewApprox) ); 
           CheckIt := y_A * y_New:
           newDelta:=abs( y_A - y_New ):

           if newDelta <= delta then
              print([NewApprox, y_New]):
           elif CheckIt > 0 then

           if flag1 > limit1 then
              print(`*** Exceeds limits****`);    


print(`Number of Roots Found = `,Num_points):
plot(f(x), x=interval_begin..interval_end, color=red, thickness=1);


I'm working on several worksheets where I solve numerical DE's. I would like to compare the solutions by plotting in one graph So I though it would be nice to have a seperate worksheet for this.

My approach:

Statements from my original sheets:
> solution:=dsolve({ODES,ICS},maxfun=0,numeric,range=0..10,output=listprocedure):
> X__br:=eval(xbr(t),solution):

I can plot these variables in my original sheet with plot(X__br(t),t=0..10).
However, when I transfer it to a different sheet by:

> save X__br, "filename.mpl"

And then load it into a new sheet by
> read("D:/.../filename.mpl");

There seems to be no possible way to plot the X__br(t) anymore. How can I solve this?
My error when I try plot(X__br(t),t=0..10) is:

Warning, unable to evaluate the function to numeric values in the region; see the plotting command's help page to ensure the calling sequence is correct


Just a detail : How to get the power write in title ( without displaying the symbol ^ ) ?



How do I change the color of the title of my graph? Thanks


Could you explain me where the third output comes from?

print~([a, b]);
   [ ]

Thanks in advance

5 6 7 8 9 10 11 Last Page 7 of 1519