Dear Users!

Hoped everyone fine with everything! I want to define a square matrix P whose elements are

p[i,j]=<φ[i],ψψ[j]>

printlevel := 2; for i from 0 while i <= 2^k-1 do for j from 0 while j <= M-1 do varphi[i, j] := t^j end do end do;
printlevel := 2; for i from 0 while i <= 2^k-1 do for j from 0 while j <= M-1 do phi[M*i+j+1] := varphi[i, j] end do end do;
printlevel := 2; for i from 0 while i <= 2^k-1 do for j from 0 while j <= M-1 do psi[i, j] := 2^((1/2)*k)*sqrt(2*j+1)*(sum((-1)^(j+i1)*factorial(j+i1)*(2*t-i)^i1/(factorial(j-i1)*factorial(i1)^2), i1 = 0 .. j)) end do end do;
printlevel := 2; for i from 0 while i <= 2^k-1 do for j from 0 while j <= M-1 do `&psi;&psi;`[M*i+j+1] := psi[i, j] end do end do;
 

take k=2,M=3

Hello Guys, I hope you are all fine. I have been struggling with creating an animation of the points (x,y) in maple. I have tried this example 
L := [[1, 1], [3, 2], [3.4, 6], [5, 3, 7], [3, 7, 9, 1], [2, 6, 8, 4, 5]];
animate(PointPlot, [L[trunc(t)]], t = 1 .. 6, frames = 150)
but in my case it shows two points at different location means it takes x and y seperate value and showed it on 1 and 2 on x axis but i want to animate it as the location of point. Please help me. 
Thank you in anticipation.

Help me about Tabulate and append text below.

file help: text_and_tabulate.mw

Dear Users!

Hope everyone is fine. I am want use the command of isolate only red color term. Then I need help to simplify the expression term by term like simplify the rational expressions, combine powers, simplify radicals etc..

x*c*(diff(f(eta), eta))*epsilon/(-epsilon*t+1)^2+(1/2)*x*c^2*(diff(f(eta), eta, eta))*y*epsilon/((-epsilon*t+1)^3*sqrt(c/(nu*(-epsilon*t+1)))*nu)+x*c^2*(diff(f(eta), eta))^2/(-epsilon*t+1)^2-sqrt(c*nu/(-epsilon*t+1))*f(eta)*x*c*(diff(f(eta), eta, eta))*sqrt(c/(nu*(-epsilon*t+1)))/(-epsilon*t+1) = a*x*epsilon/(-epsilon*t+1)^2+a^2*x/(-epsilon*t+1)^2+x*c^2*(diff(f(eta), eta, eta, eta))/(-epsilon*t+1)^2+sqrt(2)*GAMMA*x^2*c^3*(diff(f(eta), eta, eta))*sqrt(c/(nu*(-epsilon*t+1)))*(diff(f(eta), eta, eta, eta))/(-epsilon*t+1)^3+A*g*beta[T]*theta(eta)*T[w]-A*g*beta[T]*theta(eta)*T[infinity]+A*g*beta[C]*phi(eta)*C[w]-A*g*beta[C]*phi(eta)*C[infinity]-sigma*B^2*x*c*(diff(f(eta), eta))/(rho*(-epsilon*t+1))+sigma*B^2*a*x/(rho*(-epsilon*t+1))

I am waiting your postive answer on it.

please help me out with maple code/script for vander pol limit cycle using its differential equation..

How do I plot the following with Maple 2015?

Dear users!

Hope everyone should be fine. I have some query regarding to set title of the following problem in 2D

beta := 1; alpha[1] := .5

f := beta+alpha[1]+sin(x*beta);

plot(f, x = 0 .. Pi, axes = boxed, size = [300, 270], labels = [x, "f ' (x)"], labelfont = ["Arial", 10, Bold], labeldirections = [horizontal, vertical], title = "&beta; = 1, &alpha;[1] = 0.5\n", titlefont = ["Arial", 10, Bold]);

I got the following graph

but I need the following one

Want to shift title on right. Want to write alpha[1] and f like in figure 2. 

I am waiting your response. Thanks

Dear users! 

Hope everyone fine with everything. I want to chage the format of x-axis in the attached file. I need the x-axis in degree like (-180^o, -120^o, -60^o, 0 , 60^o, 120^o, 180^o). Also I want the lable on y-axis as:

I am waiting your positive respons. 

Help_new.mw

I wanted to make a 3d plot with lines linking points, but I am struggling with it
(I couldn't find anything in the plot3d documentation, so it might not be supported, which would be odd as its such a normal thing to want to do)

I have tried with no success:

plot3d([1, 2, 3], [4, 5, 6]);
plot3d(Matrix(2, 3, {(1, 1) = 1, (1, 2) = 2, (1, 3) = 3, (2, 1) = 4, (2, 2) = 5, (2, 3) = 6}))

Any thoughts?

hi, i am trying to get a list of equations using the coeff function. I have the following equation:

restart;
Lambda:=-(1/8)*(D[4](a))(x, y, t, u(x, y, t), diff(u(x, y, t), x), diff(u(x, y, t), y), diff(u(x, y, t), t), diff(u(x, y, t), x, x), diff(u(x, y, t), y, x), diff(u(x, y, t), x, t), diff(u(x, y, t), y, t), diff(u(x, y, t), t, t))*(diff(u(x, y, t), x))-(1/8)*(D[6](a))(x, y, t, u(x, y, t), diff(u(x, y, t), x), diff(u(x, y, t), y), diff(u(x, y, t), t), diff(u(x, y, t), x, x), diff(u(x, y, t), y, x), diff(u(x, y, t), x, t), diff(u(x, y, t), y, t), diff(u(x, y, t), t, t))*(diff(u(x, y, t), y, x))-(1/8)*(D[7](a))(x, y, t, u(x, y, t), diff(u(x, y, t), x), diff(u(x, y, t), y), diff(u(x, y, t), t), diff(u(x, y, t), x, x), diff(u(x, y, t), y, x), diff(u(x, y, t), x, t), diff(u(x, y, t), y, t), diff(u(x, y, t), t, t))*(diff(u(x, y, t), x, t))-(1/8)*(D[5](a))(x, y, t, u(x, y, t), diff(u(x, y, t), x), diff(u(x, y, t), y), diff(u(x, y, t), t), diff(u(x, y, t), x, x), diff(u(x, y, t), y, x), diff(u(x, y, t), x, t), diff(u(x, y, t), y, t), diff(u(x, y, t), t, t))*(diff(u(x, y, t), x, x))-(1/8)*(D[12](a))(x, y, t, u(x, y, t), diff(u(x, y, t), x), diff(u(x, y, t), y), diff(u(x, y, t), t), diff(u(x, y, t), x, x), diff(u(x, y, t), y, x), diff(u(x, y, t), x, t), diff(u(x, y, t), y, t), diff(u(x, y, t), t, t))*(diff(u(x, y, t), x, t, t))-(1/8)*(D[10](a))(x, y, t, u(x, y, t), diff(u(x, y, t), x), diff(u(x, y, t), y), diff(u(x, y, t), t), diff(u(x, y, t), x, x), diff(u(x, y, t), y, x), diff(u(x, y, t), x, t), diff(u(x, y, t), y, t), diff(u(x, y, t), t, t))*(diff(u(x, y, t), x, x, t))-(1/8)*(D[11](a))(x, y, t, u(x, y, t), diff(u(x, y, t), x), diff(u(x, y, t), y), diff(u(x, y, t), t), diff(u(x, y, t), x, x), diff(u(x, y, t), y, x), diff(u(x, y, t), x, t), diff(u(x, y, t), y, t), diff(u(x, y, t), t, t))*(diff(u(x, y, t), y, x, t))-(1/8)*(D[9](a))(x, y, t, u(x, y, t), diff(u(x, y, t), x), diff(u(x, y, t), y), diff(u(x, y, t), t), diff(u(x, y, t), x, x), diff(u(x, y, t), y, x), diff(u(x, y, t), x, t), diff(u(x, y, t), y, t), diff(u(x, y, t), t, t))*(diff(u(x, y, t), y, x, x))+(1/8*(-3*(D[8](a))(x, y, t, u(x, y, t), diff(u(x, y, t), x), diff(u(x, y, t), y), diff(u(x, y, t), t), diff(u(x, y, t), x, x), diff(u(x, y, t), y, x), diff(u(x, y, t), x, t), diff(u(x, y, t), y, t), diff(u(x, y, t), t, t))-4))*(diff(u(x, y, t), x, x, x))-(1/8)*(D[1](a))(x, y, t, u(x, y, t), diff(u(x, y, t), x), diff(u(x, y, t), y), diff(u(x, y, t), t), diff(u(x, y, t), x, x), diff(u(x, y, t), y, x), diff(u(x, y, t), x, t), diff(u(x, y, t), y, t), diff(u(x, y, t), t, t))+(1/4)*(D[8](a))(x, y, t, u(x, y, t), diff(u(x, y, t), x), diff(u(x, y, t), y), diff(u(x, y, t), t), diff(u(x, y, t), x, x), diff(u(x, y, t), y, x), diff(u(x, y, t), x, t), diff(u(x, y, t), y, t), diff(u(x, y, t), t, t))*(diff(u(x, y, t), x, x, x));

Now I am trying to get the coefficients of the derivatives of u(x,y,t). So by hand I have done the calculation of taking the coeffient of diff(u(x,y,t),x,x,x) and setting this equal to zero, and then the coefficient of diff(u(x,y,t),x,x,t) and so on. However when I use the coeff function, it will only allow me to compute the coeff of diff(u(x,y,t),x,x,x) and others of third order. However it will not let me use the deriviateves of u of second and first order.

coeff(Lambda,diff(u(x,y,t),x));

this will not work but 

coeff(Lambda,diff(u(x,y,t),x,x,x));

will work.

Any help would be great thanks.

I have a relatively complicated ODE that i am plotting. One of the variables in particular (B[2211], purple line in the graph) should under go an exponential-like decay to zero, but instead flies off into negative territory  (see graph below) despite having a very simple equation:

problem variable rate of change= -problem variable*constant +linear combination of variables that  are always posotive.

My intuition is that this is because i have somehow used the odeplot tool wrong - possibly due to a problem with stepsizes.

Here is a MWE i've made:

aa_problem_MWE.mw

Hi everybody,

Please take the example given in the help pages of DocumentTools[Tabulate] (the one where a cardinal sine is plotted)
Change plot(sin(x)/x) by plot(sin(x),/x, legend=sinc(x))

The legend doesn't appear if the list (named "A" in the help page) is displayed through DocumentTools:-tabulate.

Is it possible to circumvent this problem ?

TIA

Hi all,

I don't know why this file doesn't work anymore! It was ok last week. I even reinstall the Maple on my laptop, but it went wrong again...

Please help me.

Plot_Diff_Temperature.mw

"There were problems during the loading process. Your worksheet may be incomplete."

Pls fix this problem for me. I tried various things. The text in the worksheet is in Danish and contains danish letters (æøå) which may have caused the problem.

Maple Worksheet - Error

Failed to load the worksheet /maplenet/convert/1._Opgaver.mw .
 

Download 1._Opgaver.mw

Download 1._Opgaver.mw

Hi everybody,

I have written a module (let's say MyModule) that  I use as a package in a worksheet (with(MyModule)).

At some point in my worksheet I call the procedure MyProc which is part of MyModule.
I find it not to work as expected. So I modify it within MyModule and generate again the archive which contains this module.

Because my worksheet does a lot of things before the call to MyProc, I would like to test quickly the above modifications.
The idea is to do :
unwith(MyModule):  # packages() no longer exhibits its name
with(MyModule);      # to load the corrected one

Unfortunately, contrary to what happens with a "native" package,  the command unwith(MyModule) is ineffective:

• once done showstat(MyProc) still displays the content of the procedure, and running  the command MyProc(...) proves that it still "exists" in the worksheet
   
• forcing a reload of MyModule ( with(MyModule) ) and acanning again MyProc ( showstat(MyProc) ) reveals the code MyProc had when MyModule has been loaded for the first time.

Is it possible to "free" a user package through the "unwith" command ?

Hope to read you soon, TIA


PS : to be clearer

A worksheet contains  the definition of N procedures, plus the one of MyModule, and ends with the commands to generate an archive file named MyModule.mla.
MyModule is defined that way

MyModule := module()
option package

export Proc1 := eval(:-Proc1),
           ......
           MyProc := eval(:-MyProc),
           ......
           ProcN := eval(:-ProcN):
end module