Maple Questions and Posts

These are Posts and Questions associated with the product, Maple

Hello,

In this procedure I get the error, 

Error, `:=` unexpected

I know what it means I just can't seem to resolve it. Any help would be greatly appreciated! Thank you in advance for looking at this code!

Kind regards,

Gambia Man

HamilMat := proc (K::integer) local ni, mi, nj, mj, N, Hamil, Eigenvec, i, j, res; option remember; global Vij, U, L; N := K^2; ni := Vector(N); mi := Vector(N); nj := Vector[row](N); mj := Vector[row](N); for i to N do for j to K do res := (i+K-j)/K; if type(res, integer) = true then ni[i] := j; nj[i] := j; mi[i] := res; mj[i] := res end if end do end do; Hamil := Matrix(N, shape = symmetric); for i to N do for j from i to N do if i <> j then Hamil(i, j) := Vij(ni[i], mi[i], nj[j], mj[j]) elif i = j then Hamil(i, j) := Vij(ni[i], mi[i], nj[j], mj[j])+(1/2)*(ni[i]^2+mi[i]^2)*Pi^2/L^2 end if end do end do; return Eigenvec := Eigenvectors(Hamil, output = ['values', 'vectors']), Hamil end proc

Error, `:=` unexpected

 

 

Download small_error.mw

Hi all,

 

I tried to find the real solution of the unlinear integral equation:

 

exp(-h^2/T)*(Int(exp(-x^2/T)*BesselI(0, h*x/T)*x, x = 0 .. 1))/T

 

but I got the warning and an complex solution:

 

 solve(subs(T = 1, eq)-.99 = 0, h)

Warning, solutions may have been lost

-1.232350242*10^(-32)-1.130417828*I

 

Can anyone help me to find a real solution for this issue (if possible)...?

I would like to thank you in advance.

 

Hi. I have a commutation table, that I've got as output of Greg Reid algorithm.

Can I use to define Lie algebra for components DifferentialGeometry & Lie Algebras?

Worksheet is attached...

FCT.mw

 

In the following problem though b and c are same (except the way denominator 2 is hanfled), command ' a-b ' readily answers zero, but a-c not so. Why? Only on condition of assumption real it gives zero!

a := (1/2)*(kappa*omega^2+omega^3)*(Y+(1/2)*(-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)/(kappa*omega^2+omega^3))^2/omega:

b := (1/2)*(kappa*omega^2+omega^3)*(Y+(1/2)*(-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)/(kappa*omega^2+omega^3))^2/omega:

a-b;

0

(1)

c := (1/2)*(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/omega:

a-c;

(1/2)*(kappa*omega^2+omega^3)*(Y+(1/2)*(-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)/(kappa*omega^2+omega^3))^2/omega-(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)

"(->)"

0

(3)

``

Why the answer is not given as zero?

``

``

 

Download what_is_the_difference_between_b_and_c.mw

What difference therms b and c make for Maple? Are they not same?

Ramakrishnan V

rukmini_ramki@hotmail.com

are these quations in wave form for harmonic ocillator?

i can not plot them

plot(exp(I*2*t) - exp(-I*2*t), t=0..10);

plot(cos(t) - I*sin(t), t=0..10);

 

why same solution but their plot graph are different?

plot(cos(t) + sin(t), t=0..10);
complexplot(exp(I*t), t=0..10);

 

Hello all, first time posting here.  Anyways, I am trying to make a velocity vector plot from known u and v velocity components at known x and y locations i.e. at (0.1, 0.1), u=-0.5 m/s and v=0.43 m/s.  My calculations have been performed with a FORTRAN code, and the results output in matrix format.  What I am left with are 4 vectors: x location, y location, u velocity component, and v velocity component.  I know that MATLAB has a specific function, called "quiver," for creating a velocity vector plot from known data points (through the use of vectors containing those points).  I have trouble believing that there is not an equivalent function, however convoluted it may be, in Maple to create these plots.  I have tried to manipulate the fieldplot, vectorplot, and arrow functions to no avail.

Any ideas on how to create a velocity vector plot from known data points?  I am using Maple 18.

Here is a link to an example of what I am looking for:

https://people.sc.fsu.edu/~jburkardt/py_src/spiral_data/spiral_matplotlib.png

Thanks, Brandon.

I have used data and the curve fitting procedure to find a curve and I now have to put error bars on the graph to complete my task? I know that there is an errorplot function but I cant seem to find a way to do it with the curvefitting thing?

 

I have attached the whole worksheet, I hope someone can help! Thanks Detectors_Tutorial.mwDetectors_Tutorial.mw

 

Even if there was some method of exporting the graph to another medium that would also be very helpful

 

 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?

First 1210 1211 1212 1213 1214 1215 1216 Last Page 1212 of 2254