Hello,

I have created a test document to check an equation and I would like to feed a range of k values (ranging from 0 to say 512) in increments of 1) to the main function, and then plot the result as a graph of Zuncomp (y-axis) vs k(x-axis).

I can verify the values of the function by typing say Zuncomp(0) or Zuncomp(10) etc but how do you run this so that you do Zuncomp(k) to read all values of k and generate the output to plot?

Using Maple 13 (Windows 7...

restart;

assume(alph>0);

assume(alph,real);

f_exp:=exp(-alph*r^2);

ff_deriv:=simplify(r^2*f_exp^2);

ff:=simplify(int(ff_deriv,r=0..infinity));

It seems to understand what i'm trying to integrate but when I try and find ff it comes up with an awful expression full of erf functions and with a lim r--> infinity at the front. However on wolfram alpha I get the answer that I want: sqrt(pi/2)/(8*alpha^(3/2)).

I am trying to solve two differential equations numericaly:

I recently get a new macbook pro with retina display. Unfortunately Maple 16 seems to be not retina ready. Searching the net I found a tool called retinizer. Using it Maple 16 seems to be displaying a little better, but still the outputs of calculations looks blured. however, whenever certain texts are selected, the whole Maple window display becomes very sharp, so I guess Maple is somehow retina ready but this functionality is disabled by default. Checking the "Get Info" window...

Hi 

In trying to solve:

 a := diff(u(t), t$2)+(p^2-I+t^2)*u(t) = 0;

I get the following solutions,              

sol1 := u(t) = _C1*WhittakerM(-1/4-(1/4*I)*p^2, 1/4, I*t^2)/sqrt(t)+_C2*WhittakerW(-1/4-(1/4*I)*p^2, 1/4, I*t^2)/sqrt(t)        

It should be possible to expand these into parabolic cylinder functions and but im not sure how, i would appreciate any help.

thanks

Hi I have three differential equations: 

u := diff(P(t), t) = -7*10^(-8)*P(t)*t/(P(t)*t+R(t))^(1/2),

diff(R(t), t) = 7*10^(-8)*t^2*P(t)/(P(t)*t+R(t))^(1/2)+600*(Z(t)^2-10^5*t^3*(1/(1.15*10^12))^(2/3)*e^(-1.15*10^12))/(t*(P(t)*t+R(t))^(1/2)),

diff(Z(t), t) = -4*10^5*(Z(t)^2-10^5*t^3*(1/(1.15*10^12))^(2/3)*e^(-1.15*10^12))/(t^2*(P(t)*t+R(t))^(1/2))

and i want to solve them with initial conditions:

initial := R(0) = 0, Z(0) = 0, P(0) = P;

Hello and thank you in advance,

I am used to Mathematica and am looking to switch, conditionals in Mathematica seem easy, like say

How do you replicate this in an easy to read way without lots of nested if statements?

Hi,

r:=[3,3,3,3,4,4,4,3,3,3,3,3,3,3,3,2,4,3,1];

What's the shortest way to create a frequency table from it?

Ie, tells me elements are 1,2,3,4

the corresponding counts are: 1,1,13 and 4

and then if it can return the number with highest counts: 3

Many thanks!

Casper

Hellow, can any body give the suggestion regarding multiple plot in array

i have 

V:=(2*x+3*gamma)*y^2

d0:=subs(x=0.1,gamma1=0,V):
d1:=subs(x=0.1,gamma1=2,V):
d2:=subs(x=0.1,gamma1=4,V):

plot([d0,d1, d2], y=0..0.1,  legend = ["gamma = 0", "gamma= 2","gamma= 4"], titlefont = ["ROMAN", 15], labels = ["y ", "v"], labeldirections = ["horizontal", "vertical"], labelfont = ["HELVETICA", 10], linestyle = [solid, longdash,dashdot],color = [black, red,green], axesfont = ["HELVETICA", "ROMAN", 8], legendstyle = [font = ["HELVETICA", 9], location = bottom],axes=boxed);

and 

e0:=subs(x=0.5,gamma1=0,V):
e1:=subs(x=0.5,gamma1=2,V):
e2:=subs(x=0.5,gamma1=4,V):

plot([e0,e1, e2], y=0..0.1,  legend = ["gamma = 0", "gamma= 2","gamma= 4"], titlefont = ["ROMAN", 15], labels = ["y ", "v"], labeldirections = ["horizontal", "vertical"], labelfont = ["HELVETICA", 10], linestyle = [solid, longdash,dashdot],color = [black, red,green], axesfont = ["HELVETICA", "ROMAN", 8], legendstyle = [font = ["HELVETICA", 9], location = bottom],axes=boxed);

I want to display  this side by side  i,e set {do,d1,d2} and {eo,e1,e2} 

 thanks in advance

 exa.doc

Hi, I am trying to plot a solution curve on a vector field with an initial condition of y(0)=0 and I keep getting error messages. This is what I have so far:

with(DEtools); dfieldplot(diff(y(x), x) = x^2+y(x)^2-1, y(x), x = 0 .. 5, y = -1 .. 5, arrows = line, title = 'Slope*Field');

Thank you. 

Greetings,

I have a finite Fourier series of the form, Asin(x)cos(z) + Bsin(2x)cos(z) + Csin(3x)cos(2z) +... I would like to be able to extract the coefficients, A,B,C etc., that correspond to a specific mode. I have tried using the coeff and coeffs functions. They work for a one-dimensional Fourier series (i.e. if S= Asin(x) + Bsin(2x) + Csin(3x); is my series, then coeff(S,sin(2x) returns B). I cannot however get this to work for the 2-D case. Any suggestions?

Hello, I'm currently trying to approximate the solution of a PDE with Maple, but I'm having troubles with the boundarie conditions .
For instance, when I set : (dT/dr)(r=3,t) = 0 , everything works fine .

With (dT/dr)(r=3,t) = -T(r=3,t), it's still okay .
But if I set (dT/dr)(r=3) = -T(r=3,t)^4, Maple just ignores the second member and works as if it was 0 . There is no error notification, but I don't get a satisfying solution ..

I have occasionally a need to map a function over the components of an expression; but not just at the top level---which map() does---but to an arbitrary level. One application is e.g. to throw simplify at subexpressions of a complicated object. Sooner or later I will want to specify the deepest level also... I looked through the docs and Help but could not find anything like that. I had expected map[n]() to do something like this, but it does not.

Does Maple have a...

Dear Maple users

We are using Maple at our school. Usually activation with the Activation code works well if the computer has a good connection to the Internet. Our data instructor have however experienced a problem with two computers he cannot seem to solve. He is not able to activate Maple 16 on those computers. After a failure in activation he even reset the license.dat file. Maybe someone have an idea to how we can get those computers activated and what the problem...

hello every body

how can i solve this equation with respect to z

a, b are constants

3a(1+z)^2=2(a(1+z)^3+b)^(1/2)