quartic_12.mw

Hellow,

 any one know how to solve the second order differential equation using finite element method (Using Linear shape functions)

Hi all,

I've just started learning Maple. I can't figure it out how to transform my problem into maple syntax.

My problem includes solving a simple algebraic equation, that is,

Y[k+1] = -10*Y[k]

I would like to solve this equation for different values of k, for example, k = 0,1,2,3.

and express my results in terms of Y[0], i.e,.

Y[1] = -10*Y[0];

Y[2] = -10*Y[1] = 100*Y[0]

Y[3] = -1000*Y[0]

I don't know...

After typing

with(LinearAlgebra):
avector:=Vector(3,symbol=a):
bvector:=Vector(3,symbol=b):
'a x b'=CrossProduct(avector,bvector)

Maple writes out

relativistic_64.mwrelativistic_64.mw

relativistic_64.mwrelativistic_64.mw

A beam of neutrons bombards a reactor's wall. Consider the motion of the neutrons as a random walk on the (x,y) plane.

The neutrons then have to comply with the following conditions:

1. Only four directions of motion are possible (left, right, up or down).

relativistic_63.mwrelativistic_63.mwrelativistic_63.mwrelativistic_63.mwrelativistic_63.mw

I have written the maple code to calculate time domain value of the laplace domain equation. It contain slowly converge double series so that I used epsilon algorithm to accelerate convergence. one summation contain 7 terms and another summation terms reduce to 850 terms after the use of epsilon alorithms. I also used deHoog algorithm with 9 terms to obtain time domain data but maple give the following error:

Error: (in evalf/exp/general) object too large

Do...

Dear sir,

i want to change the horizontal and vertical axis in maple plot. on horizontal axis  i  want velocity  u(y) values and on vertical axis y values i am attaching the file also My codes are 

U:=(y^2-2*h*y)*(-(5/2)*c1/(2*h^2*a^2+5));

 c1:=1:

 h:=1:

  a0:=subs(a=1,U):

 a1:=subs(a=1.5,U):

 a2:=subs(a=2.5,U):

 a3:=subs(a=3.5,U):

 a4:=subs(a=4.5,U):

 a5:=subs(a=5,U):

  plot([a0,a1, a2,a3,a4,a5], y=0..1, title = " Velocity profile for c1=1 and h=1", legend = [alpha = 1, alpha= 1.5,alpha= 2.5,alpha= 2.5,alpha= 4.5,alpha= 5], titlefont = ["ROMAN", 15], labels = ["y ", "u(y)"], labeldirections = ["horizontal", "vertical"], labelfont = ["HELVETICA", 10], style = [POINT],symbol = [BOX,CROSS,CIRCLE,DIAGONALCROSS,DIAMOND,SOLIDBOX],color = [black, black,black,black,black,black], axesfont = ["HELVETICA", "ROMAN", 8], legendstyle = [font = ["HELVETICA", 9], location = bottom],axes=boxed);

If I have the parameters-(v,k,lambda) of a combinatorical design, I wonder whether I can get the incidence matrix of the design or not with maple.

Remark:v is the number points of the design; k is the number of points on each line; lambda is the number of points common on every two lines.

Thanks!

hello,everyone,
I have a problem to make a figure, commands are below:

restart;
1+(2*(S1*e2/(S2*e1)+S2*e1/(S1*e2)))*tanh(S2*k*d);
S1 = sqrt(beta^2-e1), S2 = sqrt(beta^2-e2);
e1 = 4, e2 = -19-.53*I;
subs(`%%`, %, k = 2*Pi/(0.633e-6), `%%%`);

I want to plot the relation between d and beta(real part named betar and imaginary part betai),but meet some trouble.
I got some points from the article is when d=400e-10, betar in the range (2.1,2.2...

amu^4+bmu^3+cmu^2+dmu+e=0

The unknown is mu

I intend to calculate the eigenvalues and eigenvectors, but for the same matrix when i executed the Eigenvectors command in mapple twice the eigenvalues generated are same both the times but the values of eigenvector are infinite at one instance whereas they are in numeric domain in the second instance.

Can anyone please tell me why this is occuring?