Hey Guys, 

I have plotted the following graph, its the position of an elementary volume before a supernovae takes place.

x₀=r₀sinθ0cosΦ₀

y0=r₀sinθ₀sinΦ₀

z₀=r₀cosθ₀

Created a mesh of 10 mesh points, divided the ranges  

0<r<1

0<θ<∏

Can somebody show me how to input this into maple please

f_i+1,j := (a_j)(f_i,j-1) + (b_j)(f_i,j) + (c_j)(f_i,j+1)

where 

a_j := [1/2 (r-q) j Delta t - 1/2 sigma^2 j^2 Delta t] 

b_j := [1 + sigma^2 j^2 Delta t + r Delta t] 

c_j := [-1/2 (r-q) j Delta t - 1/2 sigma^2 j^2 Delta t] 

with r:= 0.05, Delta t:= T/N , Delta S:=S/M, T:=0.4, q:=0, K:=1.1, S:=2, M=N:=10

Hey,

suppose I have got two quadratics varying with variable Q1 and two quadratics varying with Q2,

eg. P1=a+b*Q1^2 and P11=c+d*Q1^2

     P2=e+f*Q2^2 and P22=g-h*Q2^2

where a,b,c,d,e,f,g and h are finite constants such that curves P1 and P11 intersect at two values of Q1, one being negative and the other positive. The same is true for curves P2 and P22. Now, I can plot curves P1 and P2 on a 3D plot (x-axis = Q1 and y-axis=Q2 and z-axis=P1=P2...

Dear all,

I have a pair of equations,

-xt[1,1]+a*tanh(xt[1,1])+a[12]*tanh(xt[2,1])=0,-xt[2,1]+a[21]*tanh(xt[1,1])+a*tanh(xt[2,1])=0

and I want to calculate xt[1,1] and xt[2,1];

when i use

solve([-xt[1,1]+a*tanh(xt[1,1])+a[12]*tanh(xt[2,1])=0,-xt[2,1]+a[21]*tanh(xt[1,1])+a*tanh(xt[2,1])=0],[xt[1,1],xt[2,1]]);

warning comes out," Warning, solutions may have been lost", and the solution [0,0] I am looking for is really lost.

I am to solve this differential equation (nonlinear):

C'(t) = G[0]+0.1*sqrt(C(t))-C(t)*n+C[ude]*n

G[0] = 1.7*1.2*83

n = 6

C[ude] = 417

i need to use programming since the DE is nonlinear.

Does anybody have the solution to that?

Thanks in advance.

Best regards,

Marie

how we can convert mathematica5.2 worksheet into maple13 worksheet?

f:=x->x^3:
y[1]:=0.5:
y[3]:=subs(x=y[1],f(x)):
i := 2:
y[i] := y[i-1]-(y[1]-y[2+1])*(1/2):
i:=1:
g[i] := (y[i+1]-y[i+2])*x/(y[i]-y[i+1])+(y[i+2]*y[i]-y[i+1]*y[i+1])/(y[i]-y[i+1]):
g[i] := unapply(g[i], x):
h[i] := unapply(solve(y = g[i](x), x), y):
for i from 2 by 1 to 7 do
g[i]:=(f@h[i-1])(x);
g[i]:=unapply(g[i],x);
y[i+2]:=subs(x=y[i+1],g[i](x));
s:=solve(z=g[i](x),x);
h[i]:=unapply(s[1],z);
end do:

for i from -1 by -1 to -8 do

Is it possible to write a Maple program which prints itself and only itself?

Hey guys I'd be grateful if you can help me plot the following graph. It's the position of an elementary volume before a supernova takes place.

x₀=r₀sinθ₀cosΦ₀

y₀=r₀sinθ₀sinΦ₀

z₀=r₀cosθ₀

create a mesh of 10 mesh points, divide the ranges

0<r<any positive value

Usually Maple gives solutions in terms of radicals only up to degree=4
(for example using RootOf + allvalues).

Using 'irreduc' (to test first) and 'galois' (to check for the Galois
group, if degree <= 9) I have cases, where the result implies, that by
theory the roots can be given through radicals:

Though the results are difficult to read in my case they are 'C(6)' or
'C(8)', the cyclic groups of that order - thus abelian and IIRC those
groups are solvable.

Hi all,

I have maple 13 and have been trying to set up some plot of a single point on a plane, to teach plotting in two dimensional plane.The figure at the following page may help you to get an insight of what I am trying to do in maple:

http://en.wikipedia.org/wiki/Cartesian_coordinate_system

Any advice?

Regards!

Sohail

Hi,

Basically I have two curves, P1 and P2, which are dependent on Q1 and Q2 in a square manner. For example, P1:=1+Q1^2 and P2:=5-Q2^2. I now want to plot P1 and P2 on the same graph with the values of P1 and P2 on the same z-axis and Q1 and Q2 along the x- and y-axes. Consequently I would get two orthogonal parabolic curves. plot3d([P1,P2],Q1=-1..1,Q2=-1..1) doesn't work because it obviously draws two surfaces rather than two curves for P1 with Q2 always equal to...

Can anybody please help :

a_j := [-1/2 (r-q) j Delta t + 1/2 sigma^2 j^2 Delta t] / 1 + r Delta t

b_j := [1 - sigma^2 j^2 Delta t] / 1 + r Delta t

c_j := [1/2 (r-q) j Delta t + 1/2 sigma^2 j^2 Delta t] / 1 + r Delta t

with r:= 0.05, Delta t:= T/N , Delta S:=S/M, T:=0.4, q:=0, K:=1.1, S:=2, M=N:=10

Then I defined i:=N; such that f_{i,j}:=max(K - j*Delta S, 0);

And next I defined g_{k,j}:=f_{i,j}...

Could somebody help me I am having trouble with an if statement.

I have defined three functions a.b and c in terms of one variable j.

I have correctly defined g_{k,j} as well as other functions and now I have run into trouble with my if statement.

p:=(k,j)->if k=0 then max(K-j*S,0) else g_{k,j} end if

I think its clear what i am trying to achieve in the if statement as I want to be able to enter p(1,4), p(4,3) etc. . .  after the if statement...

Please help me with the following.

I have this function

y(theta,t):=1/4*Pi + A*sin(m*theta) + A*sin(m*theta)* (exp(-1/4 *alpha^2*beta^2)-2)/(4*Pi)*t;