Hello guys

I have a problem in runge kutta method for coupled differential equations .For example consider these equations:

a*Φ''''(x)+b*Φ''(x)+c*Φ(x)+d*Ψ''(x)+e*Ψ(x):=0

d*Φ''(x)+e*Φ(x)+j*Ψ''(x)+h*Ψ(x):=0

The boundry conditions are:

Φ(0.83)=1;Φ'(0.83)=0;Φ(-0.83)=1;Φ'(-0.83)=0;Ψ(0.83)=1;Ψ(-0.83)=1

Φ(x) and Ψ(x) are unknown functions and a,b,c,d,e,j,h are constants . I want to get answers by runge kutta numerical method .

i used this form of algorithm but i cant do it :

please help me to solve these types of equations.

restart:

FF := Q-1+(1/5)*K*dp^3*h^5+(1/3)*dp*h^3+h+h1*h;

DDP:=[solve(FF,dp)]:

h:=1+phi*cos(2*Pi*x):

h1:=2*Pi*alpha*beta*phi*cos(2*Pi*x);

beta:=1:alpha:=0:
phi:=0.5:

dpdx:=evalf(DDP[1]):

dpp:=Int((dpdx),x=0..1):

plot([subs(K=-0.1,Re(evalf(dpp)))],Q=-1..1,axes=box,color=[blue],adaptive= false);

I am unable to get any type of output. 

Cheers!

Hello guys,

I think that the title explains the question very well. Is there any function in MAPLE that allows me to generate N random numbers considering a mean value, standard deviation and a percentile?

Thank you,

hi

I need to export or save my plots during the run out of maple as jpeg or gif but I can not find any command to do it. the only way i found is right click on the plot and then export it but i need program do it by itself during the run.

thanks

I'm trying to use Rodrigues' Rotation Formula to graph a circle--by rotating a point, naturally--and I'm not sure how to implement that in Maple. For something like a function, allowing a variable--say, x--to vary and produce a curve is natural, but the evalm() function doesn't handle that very well: I can use the formula to *evaluate* the rotations, but not to use them to graph the circle.

I looked for something like solids of revolution--after all, I'm doing something pretty similar--but didn't find anything helpful; apparently there's a way to make solids of revolution, but it's pretty specific and would be difficult to adapt?

Hello Those who are in Mapleprimes.

I want to know the meaning of [].

That is,

a:=[1,2,3];

b:=a[];

What were returned was 

b:=1,2,3

This result is the same as that obtained from

b:=op(a);

What kind of meaning does a[] have?

I want to know the logic of this expression.

Can I ask you about this?

Thank you in advance.

taro yamada 

Hi every body:

how can i solve this equations(without numerical method):

eq1 := (D[1, 1](eta11))(t, a*t, a^2*t)+1.326096634*10^8*Pi^2*eta11(t, a*t, a^2*t)-3.315241586*10^7*Pi*eta21(t, a*t, a^2*t) = 0

eq2 := 2.054901810*10^13*eta21(t, a*t, a^2*t)+(D[1, 1](eta21))(t, a*t, a^2*t)-8.219607239*10^13*Pi*eta11(t, a*t, a^2*t)+4.137421500*10^8*Pi^2*eta21(t, a*t, a^2*t) = 0

eq3 := (D[1, 1](eta31))(t, a*t, a^2*t)+4.137421500*10^8*Pi^2*eta31(t, a*t, a^2*t) = 0

I am making something for my math prof as a token of appreciation. It is a spiral slide rule designed to fit in a watch face, but it can be made bigger. I threw together a program to show what I mean:

https://www.khanacademy.org/cs/spiral-slide-rule/5831476045152256

I was told that maple can be used to create things like this, but since I am new to the software, I really don't know how to put something like this together. The best I could do was to make a spiral in polar coordinates, and plot my lists of points on top of it. However, this does not annotate the points and I can't make nice looking lines (the only point options in polarplot are stars, crosses, and other things that aren't what Im looking for).

Does anyone know how I could model this, and then export it as a format where I might be able to send to a printer to get it printed on specialty paper?

Thank you very much

EDIT: I uploaded my plot here, in case you want to see it:

http://randomwrites.wikispaces.com/file/view/mediumspiral1.bmp/508315780/mediumspiral1.bmp

Please see the link below of the screen shot from Maple

https://www.dropbox.com/s/r7xn2uqnn4qfbp7/Screenshot%202014-05-12%2015.21.07.png

AOA... Dear when i expand

sum(sum(binomial(n-1, i)*x^(n-i-alpha)*(-a*n)^i*c[n]*GAMMA(n-i+1)/GAMMA(n-i-alpha+1), i = 0 .. n-1), n = ceil(alpha) .. M)

for M=2 and alpha=1/2 its answer is 

-sqrt(x)*c[1]*sqrt(-(1-x)/x)*(2*x-1)/(sqrt(Pi)*(1-x))-(1/4)*c[1]*hypergeom([3/2, 2], [3], 1/x)/(x*sqrt(x*Pi))-(4/3)*x^(3/2)*c[2]*(-(2-x)/x)^(3/2)*(2*x-1)/(sqrt(Pi)*(2-x))+(2/3)*c[2]*hypergeom([3/2, 2], [4], 2/x)/(x*sqrt(x*Pi))

which is very difficulty i want its answer in Gamma form i.e.

2*sqrt(x)*c[1]/sqrt(Pi)+(8/3)*x^(3/2)*c[2]/sqrt(Pi)-4*sqrt(x)*c[2]/sqrt(Pi)

Pl help me

Hi all,

I ploted the cylinder like this:

with(plottools);
with(plots);
display(cylinder([1, 1, 1], 1, 3), orientation = [45, 70], scaling = constrained, grid = [2, 2, 2]);

Is it possible to export nodal points of the plot?

Thank a lot.

Is there a book or a place where to start learning programming in Maple from the very beginnings?

I have some experience with RealBasic (now Xojo) and C, so I am not a complete newbe but I am at loss with Maple and I have hard times following the help system to learn.

Thank you

Martina 

AOA.. I want to generate a matrix for arbitrary value of n

I want to generate a random matrix of size m x n as follows: for each i between 1 and m, take a random p_i number according to some fixed distribution (for example, uniform distribution in [0,10]). These numbers p_i are independent. Then I'll form a random matrix whose (i,j)-entry is chosen according to a probability distribution that depends on the parameter p_i (for example, it could be Poisson(p_i)). So entries in the same row are independent.

How can I do that?

Hi everyone!

Would really appreciate some help here.

So I'm going to have my exam next week with all these math problems. We're aloud to check our answers with Maple once.

So I wonder how I would go about this with the Taylorpolynomials? I've tried getting it to work with Maple now but I can't figure it out.

Lets say I got a question like this:

"Given the intial value problem y'(x)=x(y(x))^2, y(1)=2. Calculate the 2.order Taylorpolynomial with the "development point" 1.

I know the answer is suppose to be this:2+4(t-1)+10(t-1)^2 

But I can't get that answer in Maple...any ideas?