Maple 2016 Questions and Posts

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

This worksheet is a modification to Kitonum's excellent

The mod adds the ability to display labelled contours for expressions in x and y defined parametrically.

Your comments are welcome.

How to improve menu speed in Maple 2016?

I am using Maple 2016, when I right click the equation, it really takes a long time to show the menu. It seems that the CPU usage and Memory usage is low, it should be fast as it is only a showing menu action, what's wrong with that?

The worksheet below rolls an ellipse along the y axis with constant energy.

How can the physics be enhanced to roll the ellipse along a non-linear curve (e.g. a sine curve) with constant energy?

Help create file Excel in ExcelTools, but error row 564?

Please help me? 


I am looking for the table containing the value m,vv and mm. and these values I can easily copy and paste in excel sheet. 

Plz suggest commands 



The conditional trigger in the code below does not cause a halt although odeplot shows x dipping to a value below 2.

What is the correct coding of this event to cause a halt when x becomes less than 2?

Does any Maple documentation show examples of working discrete events with conditional triggers?


DE1 := diff(x(t), t, t) = -x(t)-(diff(x(t), t))+y(t);

DE2 := diff(y(t), t, t) = y(t)-(diff(y(t), t))+x(t);

ICs := x(0) = 2, (D(x))(0) = 1, y(0) = 0, (D(y))(0) = 1;

sol1 := dsolve({DE1, DE2, ICs}, {x(t), y(t)}, numeric, event_initial = true, events = [[[0, x(t) < 2], halt]]);

plots:-odeplot(sol1, [t, x(t)], t = 0 .. 3);

The animation in this worksheet seems realistic but I wonder if the physics which is implemented is correct since there is no explicit mention of rotational energy.


this is my code:



q:=0; h(x,y):=Heaviside(g(x,y)-q);

(x, y) -> Heaviside(g(x, y))

p := int(int(h(x, y)*exp((-x^2-y^2)*(1/2))/(2*Pi), y = -infinity .. infinity), x = -infinity .. infinity);

Error, (in unknown) too many levels of recursion



I keep having this error message, is this integral impossible to process or is there someting wrong with the code?





for x to 10 do Disp[x] := plot([x+cos(a), sin(a), a = 0 .. 2*Pi], colorscheme = ["Red", "Green"]) end do;

display(seq(Disp[i], i = 1 .. x-1), insequence = true, scaling = constrained);

What coding in the color or colorscheme parameter in the plot command above would rotate the circle's colors in conjunction with the circle's movement along the x axis to give the appearance of a rolling wheel? 

how can I get a positive output in the following example

solve((hv*hw+lv*lw)/(lv+hv)*(lv+hv) = dw*(lv+hv), hv)

instead of -lv*(dw-lw)/(dw-hw) I would like lv*(dw-lw)/(hw-dw)

How can type 



op(S[1]), op(S[2]), op(S[3]), op(S[4]), op(S[5])





op(S[1]), op(S[2]), op(S[3]), op(S[4]), op(S[5])








[op(S[1]), op(S[2],op(S[3],op(S[4])]

Please help me

When I try:


I get:

Error, (in Statistics:-Sample) invalid input: expected distribution or algebraic expression with random variables, but received X

What am I doing wrong?


in maple, i can't use print? Help me. Please.

The third execution of procedure TP returns infinity but it has a real value in the integral's plot. Why is this?

Good day everyone.

Please can help me with this code on the taylor series expansion involving Fractional Differential Equation (FDE)? Particularlly, the lines highlighted blue and green respectively in relation to FDE.

Thank you and kind regards

P:=sum((a[k]*x^(k))/GAMMA(k+1-alpha), k=0..3):
assume(alpha>0,alpha < 1):
e1:=simplify(eval(P, x=q+h))=y[n+1]:
var:=seq(a[i], i=0..3):

Cc:=eval(<var>, solve(eval({M}),{var}) ):
for i from 1 to 4 do
end do:
E:=collect(Cf, [y[n+1], f[n], f[n+1],f[n+2]], recursive):
s2:=y[n+2]=simplify(eval(Cf, x=q+2*h)):
collect(%, [y[n+1], f[n], f[n+1],f[n+2]], recursive):

s1:=y[n]=simplify(eval(Cf, x=q)):
collect(%, [y[n+1], f[n], f[n+1],f[n+2]], recursive):

F[n]:=convert(taylor((D(y))(x), h = 0,12), polynom):
F[n+1]:=convert(taylor((D(y))(x+h), h = 0,12), polynom):
F[n+2]:=convert(taylor((D(y))(x+2*h), h = 0,12), polynom):




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