> restart;
> Digits := 10;
> m := 11;
> P := 100;
> alpha := 1;
> F[0] := 0;
> F[1] := epsilon;
> epsilon := 0;
> F[2] := A;
> T[0] := -T[1]/alpha-1;
> T[1] := B;
> for k from 0 to m do F[k+3] := (-(sum(F[k-r+2]*F[r]*(k-r+2)*(k-r+1), r = 0 .. k))-1+sum(F[r+1]*F[k-r+1]*(r+1)*(k-r+1), r = 0 .. k))*factorial(k)/factorial(k+3); T[k+2] := -P*(sum(F[r]*T[k-r+1]*(k-r+1), r = 0 .. k))*factorial(k)/factorial(k+2) end do;
> f := 0;
> t := 0;
>
> for k from 0 to m do f := f+F[k]*eta^k; t := t+T[k]*eta^k end do;
> print(f);
> print(t)

> with(numapprox);
> pade(f, eta, [4, 4]);

>pade(t, eta, [4, 4])

>solve({limit(pade(f, eta, [4, 4]), eta = infinity) = 0., limit(pade(t, eta, [4, 4]), eta = infinity) = 0.}, [A, B])

Hi!

Say, I got an expression that depends on two variables, x and y. How can I tell Maple, that y is actually just a (real) constant, so y does not depend on x?

Because when I apply a differentiation with the "D" - command, it would always also write out expressions, where y is differentiated w.r.t. x.

Thanks!

Hi There,

Can any one correct me the mistake in the following differentiation:

Above command gives the following error:

Error, (in simpl/abs) abs is not differentiable at non-real arguments

The differential equation I'm solving for is:

Hello fellow Maple users. 
I have been useing Maple for 5 month now and i think its a great program. But sometimes i find it a bit difficult to use. 
When i plot a function i would like to see the function on the line in the plot that represents the funtion. So if have several lines representing different functions, i would Maple to highlight the function when i hower over the function 
with the mouse .. Is this posible ? I know its posible in Inspires CAS. 

It must be posible for Maple to label the plots with the function ??? If not Maple programmers got a job...

I have a problem in excuting this differential equation in maple it takes a long time but yet no result.

> restart;

> Delta:= epsilon[2]-epsilon[1];

> epsilon[y] := epsilon[2]-(1/4)*Delta*(1-tanh(a*y))^2;

 > z:= tanh(a*y) ;

 > ODE[4]:= diff(Y(y),y,y)- ( a/2* Delta *(1-z)*(z^2-1))/(epsilon[2]- Delta*(1-z)/4)* diff(Y(y),y)-( beta^2+ mu[0]*epsilon[y]*omega^2)*Y(y) = 0;

> dsolve(ODE[4],Y(y));

does this always occur or i do have problem with my version of maple 15, 7 and 16.

Thank you, looking forward for your answers.

I have created a simple proc to return true or false depending on the value in a list of values.  now i would like to filter out the falses and return the actual value that is true from the list.

sol := dsolve([ode2, theta(0) = (1/180)*Pi*10, (D(theta))(0) = 0], numeric)

sol(0)

[t = 0., theta(t) = .174532925199433, diff(theta(t), t) = 0.]

given a numerical solution in the above form, how can I integrate a function

f(theta(t)) numerically?

I tried to isolate theta(t) first:

x:=t->rhs(op(2,sol(t)))

But I can't integrate

phi:=t->int(x(t),t=0..t)

I am having some difficulty animating the function shown in the attached file.  I am going to create an animation which will show the curve as a function of t.  My first question is that there is no way to compute K_n because the initial conditions I have are only given as arbritary functions F(z),G(z).  So I am not really sure how to proceed here.

My second question is that I also want to plot the Z dependent part of y as a function of z/b.  I have tried to incorporate this into Maple, however, all I get back is that there are 'unexpected variables present'

Thanks.

Is it possible to assign a matrix given values/expressions?

Instead of assigning each element separately like in the example below:

with(combinat):
list1 := permute([a, b, a, b, a, b], 3);
list1a := subs(b=1,subs(a=0, list1));
list1a := permute([seq(seq(k,k=0..1),k2=1..3)], 3);
list2 := permute([a, b, c, d, e, f, g, h, a, b, c, d, e, f, g, h, a, b, c, d, e, f, g, h], 3);
list3 := subs(h=18,subs(g=17,subs(f=16,subs(e=15,subs(d=14,subs(c=13,subs(b=12,subs(a=11,list2))))))));
list3 := permute([seq(seq(k,k=11..18),k2=1..3)], 3);
list5 := Matrix(nops(list1a)*nops(list3), 1);
count := 1;
for n from 1 to nops(list3) do
temp1 := subs(1=list1a[1],list3[n]);
for k from 11 to nops(list1a)+10 do
temp1 := subs(k=list1a[k-10],temp1);
od;
list5[count] := temp1;
count := count + 1;
od;
list6 := permute(list5, 2);

Error, (in combinat:-permute) 1st argument must be a list, set or a non-negative integer

I have a bit of an issue.

I'm trying to solve this:

eq1:= exp(-lambda)*exp(k)/factorial(k) > .95

lambda:=2

solve(eq1,k)

Warning, solutions may have been lost

Good afternoon sir.

I request your kind support to the above cited query.

With thanks & Regards

M.Anand

Assistant Professor in Mathematics

SR International Institute of Technology,

Hyderabad, Andhra Pradesh, INDIA.

Consider the differential equation zZ'' + Z' + a²Z = 0,  where Z = Z(z).  Using the change of variables x = \sqrt{z/b}with b a constant,  obtain the differential equation Z'' + (1/x)Z' + c²Z = 0, where Z = Z(x) and c = 2a \sqrt{b}.

I tried Maple help and it offers the dchange command, and what I have tried is shown below;

with(PDEtools):

DE:= ...

tr:= {z = x²b}

dchange(tr, DE)

This did not return anything however.  I am thinking I need to specify that b is a constant, however, I am a little unsure on how to do this. Is the above the correct way to proceed?  I don't see how I have specified anywhere that in the final PDE, I require Z=Z(x).

Thanks for any help.  This is my first post here, so apologies for the typesetting. If there is inbuilt latex, I will use it next time.

Hi;

Trying to find the maximum values of the implicit function.. I tree the following commands.

l:=1:alpha:=1:b:=100:k:=20:

(alpha+(l+alpha)*u+alpha*k*u^2)*a=u*(alpha+(l+alpha)*u+alpha*k*u^2)*(1+l*alpha*b/((alpha+(l+alpha)*u+alpha*k*u^2)));

with(plots):

implicitplot((alpha+(l+alpha)*u+alpha*k*u^2)*a=u*(alpha+(l+alpha)*u+alpha*k*u^2)*(1+l*alpha*b/((alpha+(l+alpha)*u+alpha*k*u^2))),a=0..100,u=0..20,numpoints=9000000,color=black,axes=boxed,font=[1,1,18],thickness=2,tickmarks=[3,3],view=[0..12,0..12],labels=[a,u]);

maximize((alpha+(l+alpha)*u+alpha*k*u^2)*a=u*(alpha+(l+alpha)*u+alpha*k*u^2)*(1+l*alpha*b/((alpha+(l+alpha)*u+alpha*k*u^2))), N=0.5..1, location);

Any help I'll appriciated