Maple Questions and Posts

These are Posts and Questions associated with the product, Maple

 with small step rather than two points?

I had writen my questions into this code hw2_finished.mw And I think there's maybe something wrong with A=[0..130]

I think you can understand what I want by reading this code but if you are confused please let me know.

Dear all

I have query related to generating commutator table from given Lie algebra spanning symmetries of PDE, I have ten vectors that I need to commute with each other to constitute commutator table. I have calculated all the commutations manually using "symmetrycommutator" command.

My queries is, can I convert this operation into short loop type program which automatically calculate cummutation of vectors from 1 to 10.

The Maple file is attaches with this query for reference.

Commutator_for_Infinitesimal_Generators_for_variable_coefficient_Zakharove_Kuznetsov_Equation.mw

Hello All,

 

I am doing research for my master's thesis and I am trying to type a PDE into maple and have been struggling for about 2 weeks now. I am new to Maple but have had a lot of Matlab experience. I've seen a couple of ways to do variable changes and things but I'm still.

 

These are the substitutions I'm trying to do.

into this equation, but and some of these terms are 0 (x*d/dt are 0, P, and M_e, and q_ye).

 So it is basically a Euler-Bernoulli beam with a free vibration.

I tried to use dchange to do the substitutions but I am having trouble defining ...

 or doing both

 

 

 I see why there is an error but I don't know how to fix it.

 

I've also tried

Which is close but still doesn't look quite right. because d/dt(omega*t/2) -> omega/2

 

Any thoughts?

 

Thanks in advance

Dear all,

I am using Maple to perform numerical integrations. When the final index in the loop is set to 5, the computation is fast and the results are quickly delivered. When I set a number higher than that, even 6, the program gets really slow and often crashes.

I herewith attach the script I use to generate the results. I guess that there should be a problem of memory management and I tried to use gc() as suggested in some forms but without success. I would appreciate it if someone here could explain the reason behind the problem.

Thank you,
question.mw

restart; Ts := 1.; sigma := 1.; C := 1.; B := 2./(1+C); with(inttrans); beta := B*Ts*omega; assume(Tb > 0); assume(u >= 0); FzzS := -(3/2)*u^3*((2*u+I*beta)*(exp(2*sigma*u)+sigma^2*exp(2*u))-4*sigma*u*((1+sigma)*u-1))/((2*u+I*beta)^2*exp((2*(1+sigma))*u)-4*u^2*((1+sigma)*u-1)^2); InvFzzS := simplify(invfourier(FzzS, omega, t)); logTimeMin := -2; logTimeMax := -1; NumSteps := 6; logTimeStep := evalf(1.0*(logTimeMax-logTimeMin)/NumSteps); curdirectory(); A := matrix(NumSteps, 2); T1 := Array(1 .. NumSteps); AF := Array(1 .. NumSteps); for i to NumSteps do logTime := evalf(i*logTimeStep+logTimeMin); curTime := evalf(10^logTime); A[i, 1] := curTime; A[i, 2] := evalf(Int(eval(InvFzzS, t = curTime), u = 0 .. infinity, epsilon = 10^(-5))); T1[i] := A[i, 1]; AF[i] := A[i, 2] end do

 

 

``

 

Download question.mw

hi,

here a comlicated formula,how i simplify

thanks  a lot.

``

f := (kappa*omega^2+omega^3)*(Y+(-sqrt(N)*omega^(3/2)*sin(theta[2])*cos(varphi[2])*lambda__b+sqrt(N)*omega^(3/2)*sin(theta[1])*cos(varphi[1])*lambda__a)/(2*(kappa*omega^2+omega^3)))^2/(2*omega)+(-kappa*omega^2+omega^3)*(X+(sqrt(N)*omega^(3/2)*sin(theta[2])*cos(varphi[2])*lambda__b+sqrt(N)*omega^(3/2)*sin(theta[1])*cos(varphi[1])*lambda__a)/(2*(-kappa*omega^2+omega^3)))^2/(2*omega)+(Omega*N*cos(theta[2])*omega+Omega*N*cos(theta[1])*omega-P__X^2*kappa+P__X^2*omega+P__Y^2*kappa+P__Y^2*omega)/(2*omega)-(sqrt(N)*omega^(3/2)*sin(theta[2])*cos(varphi[2])*lambda__b+sqrt(N)*omega^(3/2)*sin(theta[1])*cos(varphi[1])*lambda__a)^2/(8*omega*(-kappa*omega^2+omega^3))-(-sqrt(N)*omega^(3/2)*sin(theta[2])*cos(varphi[2])*lambda__b+sqrt(N)*omega^(3/2)*sin(theta[1])*cos(varphi[1])*lambda__a)^2/(8*omega*(kappa*omega^2+omega^3))

(1/2)*(kappa*omega^2+omega^3)*(Y+(-N^(1/2)*omega^(3/2)*sin(theta[2])*cos(varphi[2])*lambda__b+N^(1/2)*omega^(3/2)*sin(theta[1])*cos(varphi[1])*lambda__a)/(2*kappa*omega^2+2*omega^3))^2/omega+(1/2)*(-kappa*omega^2+omega^3)*(X+(N^(1/2)*omega^(3/2)*sin(theta[2])*cos(varphi[2])*lambda__b+N^(1/2)*omega^(3/2)*sin(theta[1])*cos(varphi[1])*lambda__a)/(-2*kappa*omega^2+2*omega^3))^2/omega+(1/2)*(Omega*N*cos(theta[2])*omega+Omega*N*cos(theta[1])*omega-P__X^2*kappa+P__X^2*omega+P__Y^2*kappa+P__Y^2*omega)/omega-(1/8)*(N^(1/2)*omega^(3/2)*sin(theta[2])*cos(varphi[2])*lambda__b+N^(1/2)*omega^(3/2)*sin(theta[1])*cos(varphi[1])*lambda__a)^2/(omega*(-kappa*omega^2+omega^3))-(1/8)*(-N^(1/2)*omega^(3/2)*sin(theta[2])*cos(varphi[2])*lambda__b+N^(1/2)*omega^(3/2)*sin(theta[1])*cos(varphi[1])*lambda__a)^2/(omega*(kappa*omega^2+omega^3))

(1)

``

(1/2)*(kappa*omega^2+omega^3)*(Y+(-N^(1/2)*omega^(3/2)*sin(theta[2])*cos(varphi[2])*lambda__b+N^(1/2)*omega^(3/2)*sin(theta[1])*cos(varphi[1])*lambda__a)/(2*kappa*omega^2+2*omega^3))^2/omega

(2)

``

    f is a complicated function,i want to make it more simplify,but i want to keep square style,

 let coefficients of X and Y keep one unit,and simplify terms  containd special symbol of omega

 

Download Q1119.mw

it what i wanted.

How using Maple one can obtain commutator and adjoint table for given set of infinitesimal generators spanning Lie symmetries for PDE?

 

 

The plot command outputs a graph of symbols where the line is made up of "H"s, vertical lines on the axes are "+"s, horizontal lines on the axes are "-"s, and intercepts are "*"s.

 

How do I fix this problem so that the graph displays a graph normally?

Hello,

I need to be able to find the difference between the ith+1 and the ith element of a sequence and store that difference in another list. My code so far doesn't seem to be able to do that. My code is below. Thank you in advance to anyone taking their time looking at this, any help is greatly appreciated!



Worksheet_Assignment_4_Ex_2_(v).mw

 

Kind regards,

Gambia Man

The Maple 2015.2 update is now available for Windows and Linux users. It provides official support for Windows 10, connectivity to the latest MATLAB release, support for MapleSim 2015.2, updates to Physics, and a variety of small improvements throughout the product.  It is available through the Check for Updates system, and is also available from our website on the Maple 2015.2 download page.

The Mac update will be released as soon as we have finished fixing the problems on Mac OS X 10.11. We’re working on it.

eithne

Hello I am trying to plot this differential equation,  have not had any success can any one help me.


plot(q(t)= Ce^-1/2 + sin2t +cos4t)

not sure what i am doing wrong and i just used the plot command

Thanks

Hi everyone. Please excuse me in advance, as I am new on this website and to Maple.

So I am using Maple 16 and I want it to give me the derivative of a function that has constants in it, for example f(x)=(x-a)(x-b). I wrote the (probably stupid) commands on Maple:

> f(x):=(x-a)(x-b);

> diff(f(x), x);

but the result is not the expected (2x-a-b), but rather (D(x))(x-b)-(D(a))(x-b). What would be the right things to write to get what I need?

 

Thank you in advance,

David

Hello,

I'm trying to improve the fvisual aspect of a worksheet.

For some inputs or texts, I would like to reuse some formats.

Is a way to duplicate formatting such as for software like word or others ?

In others words, is there in Maple this kind of option   Duplicate formatting ?

Thanks a lot for your help.

 

Has anyone tried to run the following in Maple command-line mode (i.e. in terminal window, type "maple" to start it without the graphic interface),

"

expr1:=t1+t2+t3+t4+t5+t6+t7+t8+t9+t10+t11+t12+t13+t14+t15+t16+t17+t18+t19+t20+t21+t22+t0-t0+t23;
expr2:=t1+t2+t3+t4+t5+t6+t7+t8+t9+t10+t11+t12+t13+t14+t15+t16+t17+t18+t19+t20+t21+t22+t0-t0+t23;
print(expr1-expr2);

"

Surprisingly, I didn't get "0" with my Maple 17 (under Linux platform) or 18 (under Mac OSX platform). Can anyone help me confirm this?

Hello everyone,

I'm working on a simulation for standing wave to prove that the combination of 2 waves in opposite direction can create standing wave. So I use these:

> restart;
> with(plots):
> W1:=A*cos(omega*t-k*x);

> W2:=A*cos(omega*t+k*x);

> W:=W1+W2;

> SW:=(A,omega,k)->animate(plot,[{W1,W2,W},x=-4..4,y=-4..4,color=[red,green,blue],scaling=constrained],t=0..5,frames=10);

> display(SW(2,2*Pi,5),insequence);

It did work if SW is a function with one variable, now I need 3 variables (A,omega,k);

It said: "Plotting error, empty plot"

Please show me my mistake or an another method. Thank you

First 1180 1181 1182 1183 1184 1185 1186 Last Page 1182 of 2224