that is ok I sloved the problem,thank you very much

it is ok after restart.

w := proc (x, theta, t) options operator, arrow; sin(m*Pi*x/L)*cos(n*theta)*exp(I*omega*t-I*g*x) end proc;l[1] := (D[`$`(1, 2)](w))(1, (1/3)*Pi, 2)+(1-sigma)*(D[`$`(2, 2)](w))(1, (1/3)*Pi, 2)/(2*R^2)-K*rho[C]*h*(D[`$`(3, 2)](w))(1, (1/3)*Pi, 2); L[1] := simplify(l[1]/exp((2*I)*omega-(1.326880917*10^14*I)*gg*kk))

above is ok

but later I want to continue l[2] := (1+sigma)*(D[1, 2](w))(1, (1/3)*Pi, 2)/(2*R); L[2] := simplify(l[2]/exp((2*I)*omega-(1.326880917*10^14*I)*gg*kk)); l[3] := -sigma*(D[1](w))(1, (1/3)*Pi, 2)/R; L[3] := simplify(l[3]/exp((2*I)*omega-(1.326880917*10^14*I)*gg*kk)), it still shows the error '  Error, (in simplify/table) too many levels of recursion
 

that is why?

Thank you

I have changed the semicolon into colon, That [animate, ..., tubeplot] dissapeared

but, I still can not get the fig of the relationship of kk and w,and there is no error shown by maple

can you help me ?

I am sorry ,I made a mistake, I donot know why the system upload another program, I am very sorry about that.

The question is that

I used the order implicitplot,but it disappeared that

[animate, animate3d, animatecurve, arrow, changecoords, complexplot, complexplot3d, conformal, conformal3d, contourplot, contourplot3d, coordplot, coordplot3d, densityplot, display, dualaxisplot, fieldplot, fieldplot3d, gradplot, gradplot3d, graphplot3d, implicitplot, implicitplot3d, inequal, interactive, interactiveparams, intersectplot, listcontplot, listcontplot3d, listdensityplot, listplot, listplot3d, loglogplot, logplot, matrixplot, multiple, odeplot, pareto, plotcompare, pointplot, pointplot3d, polarplot, polygonplot, polygonplot3d, polyhedra_supported, polyhedraplot, rootlocus, semilogplot, setcolors, setoptions, setoptions3d, spacecurve, sparsematrixplot, surfdata, textplot, textplot3d, tubeplot]

I really donot know how to correct it ,can you help me ?

Thank you very much

the program is in addition

View 13857_2009.12.31QQ.mw on MapleNet or Download 13857_2009.12.31QQ.mw
View file details

there is minimize(expr, opt1, opt2, ..., optn),maybe it is useful
 

I used the order implicitplot,but it disappeared that

[animate, animate3d, animatecurve, arrow, changecoords, complexplot, complexplot3d, conformal, conformal3d, contourplot, contourplot3d, coordplot, coordplot3d, densityplot, display, dualaxisplot, fieldplot, fieldplot3d, gradplot, gradplot3d, graphplot3d, implicitplot, implicitplot3d, inequal, interactive, interactiveparams, intersectplot, listcontplot, listcontplot3d, listdensityplot, listplot, listplot3d, loglogplot, logplot, matrixplot, multiple, odeplot, pareto, plotcompare, pointplot, pointplot3d, polarplot, polygonplot, polygonplot3d, polyhedra_supported, polyhedraplot, rootlocus, semilogplot, setcolors, setoptions, setoptions3d, spacecurve, sparsematrixplot, surfdata, textplot, textplot3d, tubeplot]

I really donot know how to correct it ,can you help me ?

Thank you very much

the program is in additionView 13857_Q1.mw on MapleNet or Download 13857_Q1.mw
View file details

what I do is in addition

maybe it is uesful for youView 13857_Q1.mw on MapleNet or Download 13857_Q1.mw
View file details

I have fix bb:=1;

there  is one more question, it is that when I want to get the relationship of w & kk,I used eliminate({L[A]+L[DD], S[A]+S[DD]}, omega),but I can not get a result,maybe it is over the capability maple,I am not sure?

How can I solve this problem?

maybe only after I do this one, I can plot the relation.

thank you

I show my order in additon, in the program I want to get the relationship of w & kk,but in the end the soft show that it has been over max area, what should I do,can you help me ?

Thank youView 13857_kk-w.mw on MapleNet or Download 13857_kk-w.mw
View file details

thank you

yes, the eqations is for a torus, I want to get is a tube in force

in your expressionplot3d([x, r*cos(theta),r*sin(theta)], x = 0..L, theta = 0..2*Pi, 'scaling=constrained'); I get a cylinder

now I want to make the cylinder in force,maybe

your  x                  my mu(x, theta)

your  r*cos(theta),                 my r*cos(theta)-upsilon(x, theta)*sin(theta)+omega(x, theta)*cos(theta)

your  r*sin(theta),                 r*sin(theta)+upsilon(x, theta)*cos(theta)+omega(x, theta)*sin(theta)

is it right?

I  only make a cylinder in force,

thank you very much

yes, the eqations is for a torus, I want to get is a tube in force

thank you,I solved Q1,

for the Q2, after I use 'scaling=constrained' ,I found the pic of tube is not what I expect,

I suppose it is a long tube, the length of the axial direction is L,the radius is r;

I want define

axial direction:mu := proc (x, theta) options operator, arrow; A*cos(m*Pi*x/L-n*theta) end proc

ciral direction:upsilon := proc (x, theta) options operator, arrow; B*cos(m*Pi*x/L-n*theta) end proc

radius direction:omega := proc (x, theta) options operator, arrow; C*sin(m*Pi*x/L-n*theta) end proc

and after change

fx:r*cos(theta)-upsilon(x, theta)*sin(theta)+omega(x, theta)*cos(theta)

fy:r*sin(theta)+upsilon(x, theta)*cos(theta)+omega(x, theta)*sin(theta)

fz:mu(x, theta)

the whole order as follows:

m := 3; n := 1; A := .1; B := 0; C := .1; r := 1; L := 8; mu := proc (x, theta) options operator, arrow; A*cos(m*Pi*x/L-n*theta) end proc; upsilon := proc (x, theta) options operator, arrow; B*cos(m*Pi*x/L-n*theta) end proc; omega := proc (x, theta) options operator, arrow; C*sin(m*Pi*x/L-n*theta) end proc; with(plots); plot3d([r*cos(theta)-upsilon(x, theta)*sin(theta)+omega(x, theta)*cos(theta), r*sin(theta)+upsilon(x, theta)*cos(theta)+omega(x, theta)*sin(theta), mu(x, theta)], x = 0 .. L, theta = 0 .. 2*Pi, scaling = constrained)

but I got a ring,in fact it is a tube,maybe there are some mistake,can you help me?

thank you

I do an order as follows,

m := 3; n := 1; A := .1; B := 0; C := .1; r := 1; L := 8; mu := proc (x, theta) options operator, arrow; A*cos(m*Pi*x/L-n*theta) end proc; upsilon := proc (x, theta) options operator, arrow; B*cos(m*Pi*x/L-n*theta) end proc; omega := proc (x, theta) options operator, arrow; C*sin(m*Pi*x/L-n*theta) end proc; with(plots); plot3d([r*cos(theta)-upsilon(x, theta)*sin(theta)+omega(x, theta)*cos(theta), r*sin(theta)+upsilon(x, theta)*cos(theta)+omega(x, theta)*sin(theta), mu(x, theta)], x = 0 .. L, theta = 0 .. 2*Pi)

In the pic I got,there are 2 Qs,

1. why the r/L is so large ,and if I change the r/L there is no differences between them?

2.the tube in fz direction, if I only draw the pic of with(plots); plot3d(mu(x, theta), theta = 0 .. 2*Pi, x = 0 .. L),that is what I expected pic of tube in fz,but in the first pic, the tube in fz doesnot change,although I change m

Thank you very much.