Hello , Happy weekend ... I need a help on plot of a wavefunction give below , i have made the plot using plot3d in cylindrical coordinate. But what i need now is the Polar plot of it. And it always comes with the code ensure calling sequence is correct.

psi(r,phi):= C[0, .1]*BesselI(4, 2*lambda*exp((1/2*I)*phi))*BesselJ(2, 1/r)+C[1, .1]*BesselI(6, 2*lambda*exp((1/2*I)*phi))*BesselJ(3, 1/r)+ C[2, .1]*BesselI(8, 2*lambda*exp((1/2*I)*phi))*BesselJ(4, 1/r)+ C[3, .1]*BesselI(10, 2*lambda*exp((1/2*I)*phi))*BesselJ(5, 1/r)+ C[4, .1]*BesselI(12,2*lambda*exp((1/2*I)*phi))*BesselJ(6, 1/r)+ C[5, .1]*BesselI(14, 2*lambda*exp((1/2*I)*phi))*BesselJ(7, 1/r)+C[6, .1]*BesselI(16, 2*lambda*exp((1/2*I)*phi))*BesselJ(8, 1/r)+ C[7, .1]*BesselI(18, 2*lambda*exp((1/2*I)*phi))*BesselJ(9, 1/r);

C[0, .1]:= 2.222222222*10^11;  C[1, .1]:=8.333333333*10^18 ; C[2, .1]:=6.250000000*10^26 ; C[3, .1]:= 2.127659574*10^32; C[4, .1]:= 3.125000000*10^43 ; C[5, .1]:=1.639344262*10^52 ; C[6, .1]:=1.408450704*10^61; C[7, .1]:=1.886792453*10^70;

My 1st case:
lambda:=0.1; plot3d(abs(psi(r,phi))^2,phi=0..2*Pi,r=0..5, coords=cylindrical);

2nd case: What i want.

lambda:=0.1; plot(abs(psi(r,phi))^2,r=-0..5, phi= 0..2*Pi, coords = polar);

ex := taylor(exp(x), x = 0, 5);

it can not use with subs(x=1,ex);

how to remove O(x^5) part

How can I use maple and separation of variables ,solve the wave equation u_tt=u_xx for a string of lengh 4 with u(x,0)=sin (Pi/2 x) e^-x . u_t(x,0)=0.u(0,t)=0 and u(4,t)=0.

Is there a way to solve the following integral completely?

int(int(x^2+y^2, y = -sqrt(2*a*x-x^2) .. sqrt(2*a*x-x^2)), x = 0 .. 2*a)

When I try solving that all I get is the following:

All is well when I calculate it using polar coordinates though:

int(int(r^3, r = 0 .. 2*a*cos(t)), t = -(1/2)*Pi .. (1/2)*Pi)

=         

Why doesn't maple solve the first one completely and is there a way around this?

Thanks

Calor

Hello everybody,

While doing a document about quantum mechanics, I stumble on this strange behaviour.

>with(ScientificConstants);
>GetConstant();
>GetValue()

Error message:

Error, (in ScientificConstants:-GetValue) `a[0]` is not a scientific constant object

But in the Assistant Scientific Constant, the value is there: 5.291772e-11m

Why is it not directly available like I did the first time?

Thank you in advance for your help.

--------------------------------------
Mario Lemelin

Maple 18 Ubuntu 13.10 - 64 bits
Maple 18 Win 7 -  64 bits
messagerie : mario.lemelin@cgocable.ca
téléphone :  (819) 376-0987

Hello,

I would like to ask if one can divide the maple worksheet (mw) to multiple files like series of command. It could be simillar to include in C++. I have found only save and read command but they can be use only for variables or they can be use in other way ? 

Thank you for any possible solution,

vidocq

want to write in maple code

to generate all commutative diagram 

with adjacency matrix

however, i only know a -> b, b-> c , a->d , d-> c

google no information about all commutative diagram, 

another problem is

would like to enrich theory , however, do not know how to connect property such as equations with diagram

Hello!

I am trying to write the following function in maple:

f(P(x)) = (x-1)*P'(x)-x*P(1) 

Where P(x) = 1 + x + x^2 

My approach thus far has yielded limited results; i am able to get some, but not all, of the results expected.

This is my code,( the relevant bits ):

p := x-> 1+x+x^2

f:= P -> (x-1)*P'(x)-x*P(1); 

if i then type f(1); in maple i get the expected result of -x

if however i type f(x) i get an equation that is unable to produce the expected result when simplified - seemingly, it doesn't parse for maple. The expected result would be -1 

Using unapply make's no difference. Either i have made a mistake typing in my functions - in which case i would like an ëxplanation as to what i did wrong, or seemingly there is a different approach here? Perhaps solving this problem with maple is simply not doable? - anyway my line of reasoning leads me to finally ask how i write a function within a function in this nice program :) 

Any help in this matter is greatly appreciated - i only just recently started using this program. 

I have the following d.e.:

I need to change the s variable into a different one, where the new variable is defined by

(the old s shows up in the limit of the integral)

I tried dchange, but it chokes on this as I don't have an explicit representation of s in terms of Theta.

(I know the overall solution as other people smarter than me have solved this a long time ago, but I 'd like to have the derivation to understand it).

Mac Dude

Let A and B be matrices of dimensions 3,2 and 2,3 corresondingly. It is
known Multiply(A,B)=Matrix([[8, 2, -2], [2, 5, 4], [-2, 4, 5]])

 What is Multiply(B,A)?

limit(t*(int(exp(-t*tan(x)), x = 0 .. (1/2)*Pi)), t = infinity)?

Is it possible to find it in Maple? The command

MultiSeries:-limit(t*(int(exp(-t*tan(x)), x = 0 .. (1/2)*Pi)), t = infinity);
outputs
.

Dear All,

I have problems connected with a double integration in Maple.

The simplified task with its description is in the uploaded file.

Any help will be appreciated.

Thank You! 

mp_qn_double_integra.mw

Dear all,

I would like to solve the Fredholm Integral equation, using numerical method.
This is my code.

there is a problem with subs, does not working here.
# Then, we obtain from (9) the coeficient A[n] and B[n].

Then I woulk like to recompute (2), and then compute (1).
# Puting x=m*h, in (1), how can we generate a linear Matrix from (1).

Fred.mw

Thanks

f:=C->sum(x[cat(new,c)],c=1..C);

f(4) # does not give me what i want , which is

x[cat(new,1)]+x[cat(new,2)]+x[cat(new,3)]+x[cat(new,4)];

seq(x[cat(new,c)],c=1..4);
convert([%],`+`);

This works. But it's just a simple example to illustrate my problem.

Is there a way to fix this 'cat' problem in the procedure 'f'?

Thanks,

casper

Hi everybody,

I am using maple to solve differential equations.  I have been using maple for some time and came across the infolevel[dsolve] command or option.  I tried finding the answer to my question online but cant seem to come across a straight answer.

My question is what does the integer do that you have to equal the infolevel[dsolve] command?  

For example I am using a book that states

"An even more important diagnostic tool is the infolevel[dsolve] command, which will give information on what methods are used in attempting to solve the ODE when dsolve is applied, even if unsuccessful.  An integer between 1 and 5 must be specified, with generally more detailed information being provided as the number is increased. On applying the dsolve command to ode2, the method of attack is summarized in the following output and, in this case, the general solution y(x) given with two arbitrary coefficients C1 and C2.

What does the integer between 1 and 5 mean or do? In the example above the book uses 5, why 5?.

Thanks,

Matt