A wealth of knowledge is on display in MaplePrimes as our contributors share their expertise and step up to answer others’ queries. This post picks out one such response and further elucidates the answers to the posted question. I hope these explanations appeal to those of our readers who might not be familiar with the techniques embedded in the original responses.

The Question: How to sketch a line in space?

vahid65 wanted to know how to sketch the line with equation (x-2)/3 = (y-1)/4 = (z-3)/3 in a three-dimensional space.

This question was answered using two different approaches that we will discuss.

The first approach, given by Preben Alsholm, suggested using these commands: solve, subs, and plots:-spacecurve.

Preben provided the following lines of code:

1. {(x-2)/3 , (y-1)/4 , (z-3)/3} =~t;
2. solve(%,{x,y,z});
3. L:=subs(%,[x,y,z]);
4. plots:-spacecurve(L,t=-5..5,thickness=3,color=red);

The first line sets the three expressions equal to t using the element wise operator =~.  This distributes the operation of equality over the elements of the set of expressions, forming a set of equations.

Result-1: {(x-2)/3=t, (y-1)/4=t, (z-3)/3=t}

The second line invokes the solve command. The solve command, solve(equations, variables), solves one or more equations or inequalities for their unknowns. So in this line, this command was used to solve each expression for its corresponding unknown.

Result-2: {(x= 2+3t, y=1+4t, z= 3+3t}

You may have noticed that the % symbol is used within the command. This symbol is referring to the equation label that was created chronologically last. 

The third line uses the subs command. The subs command, subs(x=a,expr), substitutes “a” for “x” in the expression expr. In this case, “a” is the reference to the equation label referenced by (%) and expr is the sequence of variables x,y,z. The square brackets around expr forms a list, and the replacement imposed by the subs command replaces each of the three names x, y, and z with their equivalents in the set of equations returned by the solve command.

Result-3: [2+3t,1+4t, 3+3t ]

Finally the last line uses the plots:-spacecurve function.  This function, spacecurve(sc,r,opts),  graphs a parametrically defined curve in three-dimensional Cartesian space. 

In this example the first parameter, sc, is replaced by L since L has the properties of a list. The parameter  r is an equation containing the parameter name and the parameter range; here it is the equation t=-5..5. Last but not least, the opts parameter (which is optional) gives the user the opportunity to specify the optional properties of the graph, properties such as thickness and color of the space curve. 

Result 4: 

Another contributer, Carl Love, suggested that these commands could be combined using the zip function, zip( f, u, v), which is an unusual application of it. The zip command applies elementwise to corresponding members of the two lists u and v, the binary function f, creating a new list r, whose members are f(u_k,v_k).

Carl provided the following code:

plots:-spacecurve(  zip(solve, [(x-2)/3 , (y-1)/4 , (z-3)/3] =~ t, [x,y,z]),  t= -5..5, thickness= 3, color= red, axes= normal);

In this case zip is applying solve to u, a list of equations, and v, a list of variables. The equations in u are those in Result-1. The zip function returns the solution of each equation for the indicated variable, that is, the new list r is the list in Result-3. This list is the first argument of plots:-spacecurve. Finally, the parameter range and opts are up to the user’s preferences.

A second approach was suggested by Kitonum. The difference between this approach and the first one is that here we are not using the solve command. Instead, we are treating the line as an intersection of two planes.

Kitonum's code is as follows:

1. L := [(x-2)*(1/3), (y-1)*(1/4), (z-3)*(1/3)]:
2. plots[intersectplot](L[1]-L[2], L[2]-L[3], x = -4 .. 4, y = -4 .. 4, z = -4 .. 4, linestyle = 1, thickness = 3, axes = normal, orientation = [60, 75], view = [-1 .. 3, -4 .. 2, -1 .. 4]);

The first statement initializes a list L that contains all the desired expressions in the given symmetric form of the line.

In the second statement Kionum used the plots[intersectplot] command. This command graphs, in three-dimensional Cartesian space, the intersection of a pair of surfaces. This example uses the calling sequence: intersectplot(expr1, expr2, x=a..b, y=c..d, z= e..f, options). 

What are expr1 and expr2 used in Kitonum’s solution? The pairwise differences of L[1] and L[2], and L[2] and L[3] are planes, the intersection of which is the desired line in space.

expr1 = L[1] – L[2] = (x-2)*(1/3)- (y-1)*(1/4)

expr2= L[2] - L[3] = (y-1)*(1/4)-(z-3)*(1/3)

The variables x,y, and z are the coordinate names and define the axes, whereas a, b, c, d, e, f define the ranges. In this case they were taken as -4 or 4.

In the options parameter Kitonum used linestyle, thickness, axes, orientation, and  view. The help page “plot3d,option” details these and other options that apply.

This blog was written by our intern Pia under the supervision of Dr. Robert Lopez. We both hope that you find this useful. If there is a particular question on MaplePrimes that you would like further explained, please let us know. 

Hi,

say I have defined some (very simple, acting on a finite set M) functions f,g,h. I want to evaluate certain permutations, so I want to make a for loop that goes through possible combinations of functions in my array A=Array([f,g,h]). (In reality my array is a little larger.) No I don't understand how I can use the elements of this array as actual functions. Even if f is defined and I can evaluate, say, f(m), I cannot use A(1)(m) to evaluate m. The output is just "f(m)". How can I make maple recognize that we are having a function here?

Thanks it advance.

Hi,

I'm not sure that I mean datatable component corectly.

I also consider that I was done somthing wrong

Thank you for advanced for any help.

wzel

Download datatable_problem.mw

Is there a way in Maple to convert a list into a string of numbers?

So [1,2,3,4,5] -> 12345

I have tried Maple help, but there are too many convert commands to try them all.
And the ones I have tried didn't work.

Thanks.

Hello,

In the creation of a list, I would like to use the assume function twice.

Here it is a print screen of my issue.

The second part of the list is not taken into account.

Do you have ideas so that my list takes into account the second term ?

Thanks a lot for your feedback

Hello! Hope every is fine. I want to expand all expression of exp of the attached file like this

exp(c[1]*t+d[1]*n-d) = exp(c[1]*t+d[1]*n)*exp(-d)

waiting your kind response.

Help.mw

Mob #: 0086-13001903838

Hello! Hope everything fine with you. Please share the command to find the max and min values of the attached function, I had tried but I was only for one variable. I am waiting your quick response.

Examples.pdf

Mob #: 0086-13001903838

Hello,

I would like to solve this equation :

tan(gamma0(t)) = tan(gamma[1](t)-theta[1](t)-psi[1](t)

I would like to select the solution of arctan by assuming the domain of variation of some variables. In my case, if -Pi/2 < gamma[1](t) - theta[1](t) - psi[1](t) < Pi/2, the arctan gives only one solution.

Can you help me to solve this equation ? 

eq:=tan(gamma0(t)) = tan(gamma[1](t)-theta[1](t)-psi[1](t))
solve(eq,gamma0(t)) assuming -Pi/2 < gamma[1](t) - theta[1](t) - psi[1](t) < Pi/2

I guess that I don't use properly the "assume" function.

The result that I would like to obain is quiet simple :

gamma0(t)) = gamma[1](t)-theta[1](t)-psi[1](t)

Thanks a lot for your help

eulermac(1/(n*ln(n)^2),n=2..N,1);  #Error
Error, (in SumTools:-DefiniteSum:-ClosedForm) summand is singular in the interval of summation


eulermac(1/(n*ln(n)^2+1),n=2..N,1);  #nonsense

Here is my Maple 16 code:

 I expected to get outuput

a [a,b,c]

a [a,c,b]

But I get no output.

Why?

2*t*exp(t) + 2*t

would like to apply op recurively 

and output list

[2,t,exp(t),2,t]

not the same ordering every time of monomials after determinant and map sign positive and op in maple 15

sometimes i need to use Reverse or Rotate List to adjust.

why ordering is different in list of monomials?

is it caused by virus?

Hi!

I wonder how is it possible to numerically evaluate two-dimensional sum, something like this:

Hi my dear friends, I am haunted by a problem of how to convert a very very very long equation to Latex or Word properly?

Thank you for taking a look at the following 'long and boring' equation and sharing your brilliant idea.

"1/12*R^2*h^3*rho*diff(diff(alpha(t),t),t)*cos(beta(t))^2*Pi+1/2*R^4*h*rho*sin(gamma(t))*diff(diff(phi(t),t),t)*Pi-1/4*R^4*h*rho*Pi*cos(beta(t))^2*diff(diff(alpha(t),t),t)-1/4*R^4*h*rho*cos(alpha(t))*diff(beta(t),t)*diff(gamma(t),t)*Pi-1/12*R^2*h^3*rho*cos(alpha(t))*diff(beta(t),t)*diff(gamma(t),t)*Pi+1/2*R^4*h*rho*diff(gamma(t),t)*cos(gamma(t))*diff(phi(t),t)*Pi+1/12*R^2*h^3*rho*sin(gamma(t))*cos(beta(t))^2*diff(diff(phi(t),t),t)*Pi+1/2*R^4*h*rho*cos(gamma(t))*diff(diff(theta(t),t),t)*cos(phi(t))*Pi-1/4*R^4*h*rho*Pi*sin(gamma(t))*cos(beta(t))^2*diff(diff(phi(t),t),t)+1/2*R^4*h*rho*sin(gamma(t))*diff(diff(psi(t),t),t)*sin(theta(t))*Pi-1/12*R^2*h^3*rho*cos(beta(t))^2*cos(gamma(t))*diff(diff(psi(t),t),t)*cos(theta(t))*sin(phi(t))*Pi-1/4*R^4*h*rho*cos(beta(t))^2*sin(alpha(t))*diff(beta(t),t)*cos(gamma(t))*diff(phi(t),t)*Pi+1/4*R^4*h*rho*cos(alpha(t))*cos(beta(t))^2*diff(beta(t),t)*diff(theta(t),t)*sin(phi(t))*Pi-1/12*R^2*h^3*rho*cos(alpha(t))*cos(beta(t))^2*diff(beta(t),t)*diff(theta(t),t)*sin(phi(t))*Pi+1/4*R^4*h*rho*cos(alpha(t))*cos(beta(t))*sin(beta(t))*diff(diff(theta(t),t),t)*sin(phi(t))*Pi-1/4*R^4*h*rho*cos(alpha(t))*diff(beta(t),t)*sin(beta(t))^2*diff(theta(t),t)*sin(phi(t))*Pi-1/12*R^2*h^3*rho*cos(alpha(t))*cos(beta(t))*sin(beta(t))*diff(diff(theta(t),t),t)*sin(phi(t))*Pi+1/12*R^2*h^3*rho*cos(alpha(t))*diff(beta(t),t)*sin(beta(t))^2*diff(theta(t),t)*sin(phi(t))*Pi-1/6*R^2*h^3*rho*sin(gamma(t))*cos(beta(t))*diff(phi(t),t)*Pi*diff(beta(t),t)*sin(beta(t))+1/2*R^4*h*rho*Pi*sin(gamma(t))*cos(beta(t))*diff(phi(t),t)*diff(beta(t),t)*sin(beta(t))-1/12*R^2*h^3*rho*cos(beta(t))^2*diff(gamma(t),t)*sin(gamma(t))*diff(theta(t),t)*cos(phi(t))*Pi-1/4*R^4*h*rho*cos(beta(t))*sin(alpha(t))*sin(beta(t))*cos(gamma(t))*diff(diff(phi(t),t),t)*Pi+1/4*R^4*h*rho*diff(beta(t),t)*sin(beta(t))^2*sin(alpha(t))*cos(gamma(t))*diff(phi(t),t)*Pi+1/12*R^2*h^3*rho*cos(beta(t))*sin(alpha(t))*sin(beta(t))*cos(gamma(t))*diff(diff(phi(t),t),t)*Pi-1/12*R^2*h^3*rho*diff(beta(t),t)*sin(beta(t))^2*sin(alpha(t))*cos(gamma(t))*diff(phi(t),t)*Pi-1/12*R^2*h^3*rho*cos(alpha(t))^2*cos(beta(t))^2*diff(gamma(t),t)*cos(gamma(t))*diff(psi(t),t)*sin(theta(t))*Pi-1/12*R^2*h^3*rho*cos(alpha(t))^2*cos(beta(t))^2*cos(gamma(t))*diff(theta(t),t)*diff(phi(t),t)*sin(phi(t))*Pi-1/12*R^2*h^3*rho*cos(alpha(t))*cos(beta(t))^2*diff(beta(t),t)*diff(psi(t),t)*cos(phi(t))*cos(theta(t))*Pi+1/4*R^4*h*rho*sin(gamma(t))*cos(beta(t))^2*sin(alpha(t))*diff(beta(t),t)*diff(theta(t),t)*cos(phi(t))*Pi-1/4*R^4*h*rho*cos(beta(t))^2*sin(alpha(t))*diff(beta(t),t)*cos(gamma(t))*diff(psi(t),t)*sin(theta(t))*Pi+1/12*R^2*h^3*rho*cos(beta(t))^2*diff(gamma(t),t)*sin(gamma(t))*diff(psi(t),t)*cos(theta(t))*sin(phi(t))*Pi+1/12*R^2*h^3*rho*cos(beta(t))^2*cos(gamma(t))*diff(psi(t),t)*diff(theta(t),t)*sin(theta(t))*sin(phi(t))*Pi-1/12*R^2*h^3*rho*cos(beta(t))^2*cos(gamma(t))*diff(psi(t),t)*cos(theta(t))*diff(phi(t),t)*cos(phi(t))*Pi-1/6*R^2*h^3*rho*cos(beta(t))*cos(gamma(t))*diff(theta(t),t)*cos(phi(t))*Pi*diff(beta(t),t)*sin(beta(t))-1/12*R^2*h^3*rho*sin(gamma(t))*cos(beta(t))*cos(alpha(t))*sin(beta(t))*cos(gamma(t))*diff(phi(t),t)^2*Pi+1/4*R^4*h*rho*diff(theta(t),t)*Pi*sin(gamma(t))*cos(beta(t))^2*cos(alpha(t))^2*diff(psi(t),t)*cos(theta(t))+1/4*R^4*h*rho*sin(gamma(t))*sin(alpha(t))*diff(beta(t),t)*diff(psi(t),t)*cos(theta(t))*sin(phi(t))*Pi-1/4*R^4*h*rho*diff(theta(t),t)*Pi*sin(gamma(t))*cos(beta(t))^2*sin(alpha(t))^2*diff(psi(t),t)*cos(theta(t))+1/12*R^2*h^3*rho*sin(alpha(t))^2*cos(beta(t))^2*cos(gamma(t))*diff(theta(t),t)*diff(phi(t),t)*sin(phi(t))*Pi+1/12*R^2*h^3*rho*sin(gamma(t))*sin(alpha(t))^2*cos(beta(t))^2*diff(psi(t),t)*diff(theta(t),t)*cos(theta(t))*Pi+1/12*R^2*h^3*rho*sin(alpha(t))^2*cos(beta(t))^2*diff(gamma(t),t)*cos(gamma(t))*diff(psi(t),t)*sin(theta(t))*Pi-1/12*R^2*h^3*rho*sin(gamma(t))*sin(alpha(t))^2*cos(beta(t))^2*diff(gamma(t),t)*diff(theta(t),t)*cos(phi(t))*Pi-1/4*R^4*h*rho*diff(theta(t),t)*Pi*cos(beta(t))^2*sin(alpha(t))^2*sin(phi(t))*diff(phi(t),t)*cos(gamma(t))+1/12*R^2*h^3*rho*sin(gamma(t))*cos(alpha(t))^2*cos(beta(t))^2*diff(theta(t),t)^2*cos(phi(t))*sin(phi(t))*Pi+1/4*R^4*h*rho*diff(theta(t),t)*Pi*cos(beta(t))^2*cos(alpha(t))^2*sin(phi(t))*diff(phi(t),t)*cos(gamma(t))+1/12*R^2*h^3*rho*sin(gamma(t))*cos(alpha(t))^2*cos(beta(t))^2*diff(gamma(t),t)*diff(theta(t),t)*cos(phi(t))*Pi+1/4*R^4*h*rho*sin(gamma(t))*cos(beta(t))*cos(alpha(t))*sin(beta(t))*cos(gamma(t))*diff(phi(t),t)^2*Pi+1/12*R^2*h^3*rho*sin(gamma(t))*sin(alpha(t))*diff(beta(t),t)*diff(psi(t),t)*cos(theta(t))*sin(phi(t))*Pi-1/12*R^2*h^3*rho*sin(gamma(t))*cos(alpha(t))^2*cos(beta(t))^2*diff(psi(t),t)*diff(theta(t),t)*cos(theta(t))*Pi+1/6*R^2*h^3*rho*cos(alpha(t))*cos(beta(t))*sin(beta(t))*diff(psi(t),t)*cos(phi(t))*diff(theta(t),t)*sin(theta(t))*Pi+1/2*R^4*h*rho*diff(gamma(t),t)*cos(gamma(t))*cos(beta(t))*sin(alpha(t))*sin(beta(t))*diff(theta(t),t)*cos(phi(t))*Pi+1/2*R^4*h*rho*cos(beta(t))*sin(alpha(t))*sin(beta(t))*diff(gamma(t),t)*sin(gamma(t))*diff(psi(t),t)*sin(theta(t))*Pi-1/2*R^4*h*rho*cos(beta(t))*sin(alpha(t))*sin(beta(t))*cos(gamma(t))*diff(psi(t),t)*diff(theta(t),t)*cos(theta(t))*Pi-1/6*R^2*h^3*rho*cos(beta(t))*cos(alpha(t))*sin(beta(t))*cos(gamma(t))^2*diff(theta(t),t)*cos(phi(t))*diff(phi(t),t)*Pi-1/2*R^4*h*rho*sin(gamma(t))*cos(beta(t))*cos(alpha(t))*sin(beta(t))*cos(gamma(t))*diff(psi(t),t)^2*cos(theta(t))^2*Pi+1/4*R^4*h*rho*sin(alpha(t))*cos(beta(t))*sin(beta(t))*cos(gamma(t))*diff(theta(t),t)^2*cos(phi(t))*sin(phi(t))*Pi-1/2*R^4*h*rho*diff(theta(t),t)*Pi*sin(gamma(t))*cos(beta(t))^2*cos(alpha(t))^2*diff(psi(t),t)*cos(phi(t))^2*cos(theta(t))-1/12*R^2*h^3*rho*cos(alpha(t))^2*cos(beta(t))^2*cos(gamma(t))*diff(psi(t),t)*diff(theta(t),t)*sin(theta(t))*sin(phi(t))*Pi+1/6*R^2*h^3*rho*sin(gamma(t))*cos(alpha(t))^2*cos(beta(t))^2*diff(psi(t),t)*diff(theta(t),t)*cos(phi(t))^2*cos(theta(t))*Pi-1/12*R^2*h^3*rho*cos(alpha(t))^2*cos(beta(t))^2*cos(gamma(t))*diff(psi(t),t)*cos(phi(t))*cos(theta(t))*diff(phi(t),t)*Pi-1/12*R^2*h^3*rho*cos(alpha(t))^2*cos(beta(t))^2*cos(gamma(t))*diff(psi(t),t)^2*sin(theta(t))*cos(phi(t))*cos(theta(t))*Pi-1/12*R^2*h^3*rho*sin(gamma(t))*cos(alpha(t))^2*cos(beta(t))^2*diff(gamma(t),t)*diff(psi(t),t)*cos(theta(t))*sin(phi(t))*Pi+1/12*R^2*h^3*rho*diff(theta(t),t)^2*sin(beta(t))*Pi*sin(gamma(t))*cos(beta(t))*cos(alpha(t))*cos(phi(t))^2*cos(gamma(t))+1/6*R^2*h^3*rho*sin(gamma(t))*cos(beta(t))*cos(alpha(t))*sin(beta(t))*cos(gamma(t))*diff(psi(t),t)^2*cos(theta(t))^2*Pi+1/2*R^4*h*rho*sin(gamma(t))*cos(beta(t))*sin(alpha(t))*sin(beta(t))*diff(psi(t),t)*diff(theta(t),t)*sin(theta(t))*sin(phi(t))*Pi+1/6*R^2*h^3*rho*diff(gamma(t),t)*cos(gamma(t))*cos(beta(t))*sin(alpha(t))*sin(beta(t))*diff(psi(t),t)*cos(theta(t))*sin(phi(t))*Pi-1/6*R^2*h^3*rho*sin(gamma(t))*cos(beta(t))*sin(alpha(t))*sin(beta(t))*diff(psi(t),t)*diff(theta(t),t)*sin(theta(t))*sin(phi(t))*Pi+1/12*R^2*h^3*rho*sin(beta(t))*Pi*cos(beta(t))*diff(psi(t),t)^2*sin(phi(t))*cos(phi(t))*sin(alpha(t))*cos(theta(t))^2*cos(gamma(t))+1/2*R^4*h*rho*diff(theta(t),t)*Pi*sin(gamma(t))*cos(beta(t))^2*cos(phi(t))*diff(phi(t),t)*cos(alpha(t))*cos(gamma(t))*sin(alpha(t))-1/2*R^4*h*rho*cos(beta(t))*cos(alpha(t))*sin(beta(t))*cos(gamma(t))^2*diff(psi(t),t)^2*sin(theta(t))*cos(theta(t))*sin(phi(t))*Pi+1/2*R^4*h*rho*sin(alpha(t))*cos(beta(t))*sin(beta(t))*cos(gamma(t))*diff(psi(t),t)*diff(theta(t),t)*cos(phi(t))^2*cos(theta(t))*Pi+1/2*R^4*h*rho*sin(gamma(t))*cos(beta(t))*cos(alpha(t))*sin(beta(t))*cos(gamma(t))*diff(psi(t),t)*sin(theta(t))*diff(phi(t),t)*Pi+1/4*R^4*h*rho*sin(gamma(t))*sin(alpha(t))*cos(beta(t))*sin(beta(t))*diff(psi(t),t)^2*sin(theta(t))*cos(phi(t))*cos(theta(t))*Pi-1/6*R^2*h^3*rho*diff(theta(t),t)*Pi*sin(gamma(t))*cos(beta(t))^2*cos(phi(t))*diff(phi(t),t)*cos(alpha(t))*cos(gamma(t))*sin(alpha(t))+1/2*R^4*h*rho*cos(beta(t))*cos(alpha(t))*sin(beta(t))*cos(gamma(t))^2*diff(psi(t),t)*diff(theta(t),t)*sin(theta(t))*cos(phi(t))*Pi-1/2*R^4*h*rho*cos(beta(t))*cos(alpha(t))*sin(beta(t))*cos(gamma(t))^2*diff(psi(t),t)*cos(theta(t))*diff(phi(t),t)*sin(phi(t))*Pi-1/6*R^2*h^3*rho*sin(alpha(t))*cos(beta(t))*sin(beta(t))*cos(gamma(t))*diff(psi(t),t)*diff(theta(t),t)*cos(phi(t))^2*cos(theta(t))*Pi-1/12*R^2*h^3*rho*sin(gamma(t))*sin(alpha(t))*cos(beta(t))*sin(beta(t))*diff(psi(t),t)^2*sin(theta(t))*cos(phi(t))*cos(theta(t))*Pi-1/6*R^2*h^3*rho*cos(beta(t))*cos(alpha(t))*sin(beta(t))*cos(gamma(t))^2*diff(psi(t),t)*diff(theta(t),t)*sin(theta(t))*cos(phi(t))*Pi+1/6*R^2*h^3*rho*cos(beta(t))*cos(alpha(t))*sin(beta(t))*cos(gamma(t))^2*diff(psi(t),t)*cos(theta(t))*diff(phi(t),t)*sin(phi(t))*Pi-1/12*R^2*h^3*rho*sin(beta(t))*Pi*sin(gamma(t))*cos(beta(t))*cos(alpha(t))*diff(psi(t),t)^2*cos(phi(t))^2*cos(theta(t))^2*cos(gamma(t))-1/6*R^2*h^3*rho*sin(gamma(t))*cos(beta(t))*cos(alpha(t))*sin(beta(t))*cos(gamma(t))*diff(psi(t),t)*sin(theta(t))*diff(phi(t),t)*Pi-1/2*R^4*h*rho*diff(gamma(t),t)*cos(gamma(t))*cos(beta(t))*sin(alpha(t))*sin(beta(t))*diff(psi(t),t)*cos(theta(t))*sin(phi(t))*Pi+1/2*R^4*h*rho*diff(diff(alpha(t),t),t)*Pi-1/2*R^4*h*rho*cos(gamma(t))*diff(diff(psi(t),t),t)*cos(theta(t))*sin(phi(t))*Pi-1/4*R^4*h*rho*Pi*sin(gamma(t))*cos(beta(t))^2*sin(theta(t))*diff(diff(psi(t),t),t)-1/4*R^4*h*rho*diff(diff(theta(t),t),t)*Pi*cos(beta(t))^2*cos(phi(t))*cos(gamma(t))+1/12*R^2*h^3*rho*sin(gamma(t))*cos(beta(t))^2*diff(diff(psi(t),t),t)*sin(theta(t))*Pi-1/12*R^2*h^3*rho*cos(alpha(t))*cos(beta(t))^2*diff(beta(t),t)*diff(gamma(t),t)*Pi+1/4*R^4*h*rho*cos(alpha(t))*cos(beta(t))^2*diff(beta(t),t)*diff(gamma(t),t)*Pi+1/4*R^4*h*rho*Pi*cos(beta(t))^2*diff(gamma(t),t)^2*cos(alpha(t))*sin(alpha(t))+1/4*R^4*h*rho*diff(theta(t),t)^2*Pi*cos(beta(t))^2*cos(alpha(t))*sin(alpha(t))-1/12*R^2*h^3*rho*diff(theta(t),t)^2*Pi*cos(beta(t))^2*cos(alpha(t))*sin(alpha(t))-1/12*R^2*h^3*rho*cos(alpha(t))*diff(beta(t),t)*diff(theta(t),t)*sin(phi(t))*Pi+1/12*R^2*h^3*rho*sin(alpha(t))*diff(beta(t),t)*cos(gamma(t))*diff(phi(t),t)*Pi-1/4*R^4*h*rho*cos(alpha(t))*diff(beta(t),t)*diff(theta(t),t)*sin(phi(t))*Pi+1/4*R^4*h*rho*sin(alpha(t))*diff(beta(t),t)*cos(gamma(t))*diff(phi(t),t)*Pi-1/12*R^2*h^3*rho*Pi*cos(beta(t))^2*diff(gamma(t),t)^2*cos(alpha(t))*sin(alpha(t))+1/2*R^4*h*rho*Pi*cos(beta(t))*diff(alpha(t),t)*diff(beta(t),t)*sin(beta(t))-1/6*R^2*h^3*rho*diff(alpha(t),t)*cos(beta(t))*Pi*diff(beta(t),t)*sin(beta(t))+1/12*R^2*h^3*rho*diff(gamma(t),t)*cos(gamma(t))*cos(beta(t))^2*diff(phi(t),t)*Pi-1/2*R^4*h*rho*diff(gamma(t),t)*sin(gamma(t))*diff(theta(t),t)*cos(phi(t))*Pi-1/2*R^4*h*rho*cos(gamma(t))*diff(theta(t),t)*diff(phi(t),t)*sin(phi(t))*Pi-1/4*R^4*h*rho*Pi*diff(gamma(t),t)*cos(gamma(t))*cos(beta(t))^2*diff(phi(t),t)+1/2*R^4*h*rho*diff(gamma(t),t)*cos(gamma(t))*diff(psi(t),t)*sin(theta(t))*Pi+1/2*R^4*h*rho*sin(gamma(t))*diff(psi(t),t)*diff(theta(t),t)*cos(theta(t))*Pi+1/12*R^2*h^3*rho*cos(beta(t))^2*cos(gamma(t))*diff(diff(theta(t),t),t)*cos(phi(t))*Pi-1/12*R^2*h^3*rho*cos(alpha(t))*cos(beta(t))*sin(beta(t))*diff(diff(gamma(t),t),t)*Pi+1/12*R^2*h^3*rho*cos(alpha(t))*diff(beta(t),t)*sin(beta(t))^2*diff(gamma(t),t)*Pi+1/4*R^4*h*rho*cos(alpha(t))*cos(beta(t))*sin(beta(t))*diff(diff(gamma(t),t),t)*Pi-1/4*R^4*h*rho*cos(alpha(t))*diff(beta(t),t)*sin(beta(t))^2*diff(gamma(t),t)*Pi+1/2*R^4*h*rho*cos(beta(t))*cos(alpha(t))*sin(beta(t))*cos(gamma(t))^2*diff(theta(t),t)*cos(phi(t))*diff(phi(t),t)*Pi+1/4*R^4*h*rho*Pi*sin(gamma(t))*cos(beta(t))^2*cos(alpha(t))^2*diff(psi(t),t)^2*sin(phi(t))*cos(phi(t))*cos(theta(t))^2+1/4*R^4*h*rho*diff(theta(t),t)*Pi*cos(beta(t))^2*sin(theta(t))*cos(alpha(t))^2*diff(psi(t),t)*sin(phi(t))*cos(gamma(t))+1/4*R^4*h*rho*cos(beta(t))*cos(alpha(t))*sin(beta(t))*diff(psi(t),t)^2*sin(theta(t))*cos(theta(t))*sin(phi(t))*Pi-1/4*R^4*h*rho*diff(theta(t),t)*Pi*cos(beta(t))^2*sin(theta(t))*sin(alpha(t))^2*diff(psi(t),t)*sin(phi(t))*cos(gamma(t))-1/4*R^4*h*rho*Pi*sin(gamma(t))*cos(beta(t))^2*sin(alpha(t))^2*diff(psi(t),t)^2*sin(phi(t))*cos(phi(t))*cos(theta(t))^2-1/4*R^4*h*rho*Pi*sin(gamma(t))*cos(beta(t))^2*diff(gamma(t),t)*sin(alpha(t))^2*diff(psi(t),t)*sin(phi(t))*cos(theta(t))-1/4*R^4*h*rho*Pi*cos(beta(t))^2*sin(theta(t))*sin(alpha(t))^2*diff(psi(t),t)^2*cos(phi(t))*cos(theta(t))*cos(gamma(t))-1/4*R^4*h*rho*Pi*cos(beta(t))^2*sin(alpha(t))^2*diff(psi(t),t)*cos(phi(t))*diff(phi(t),t)*cos(theta(t))*cos(gamma(t))-1/12*R^2*h^3*rho*sin(alpha(t))*cos(beta(t))*sin(beta(t))*cos(gamma(t))*diff(theta(t),t)^2*cos(phi(t))*sin(phi(t))*Pi+1/4*R^4*h*rho*Pi*cos(beta(t))^2*cos(alpha(t))^2*diff(psi(t),t)*cos(phi(t))*diff(phi(t),t)*cos(theta(t))*cos(gamma(t))+1/4*R^4*h*rho*Pi*cos(beta(t))^2*sin(theta(t))*cos(alpha(t))^2*diff(psi(t),t)^2*cos(phi(t))*cos(theta(t))*cos(gamma(t))+1/4*R^4*h*rho*Pi*sin(gamma(t))*cos(beta(t))^2*diff(gamma(t),t)*cos(alpha(t))^2*diff(psi(t),t)*sin(phi(t))*cos(theta(t))-1/12*R^2*h^3*rho*Pi*sin(gamma(t))*cos(beta(t))^2*cos(alpha(t))^2*diff(psi(t),t)^2*sin(phi(t))*cos(phi(t))*cos(theta(t))^2+1/6*R^2*h^3*rho*Pi*cos(beta(t))^2*sin(theta(t))*diff(psi(t),t)*diff(phi(t),t)*cos(alpha(t))*cos(gamma(t))^2*sin(alpha(t))-1/2*R^4*h*rho*Pi*cos(beta(t))^2*sin(theta(t))*diff(psi(t),t)*diff(phi(t),t)*cos(alpha(t))*cos(gamma(t))^2*sin(alpha(t))+1/12*R^2*h^3*rho*Pi*cos(beta(t))^2*diff(psi(t),t)^2*cos(phi(t))^2*cos(alpha(t))*cos(theta(t))^2*cos(gamma(t))^2*sin(alpha(t))-1/6*R^2*h^3*rho*Pi*cos(beta(t))^2*diff(gamma(t),t)*diff(psi(t),t)*cos(phi(t))*cos(alpha(t))*cos(theta(t))*sin(alpha(t))-1/4*R^4*h*rho*Pi*cos(beta(t))^2*diff(psi(t),t)^2*cos(phi(t))^2*cos(alpha(t))*cos(theta(t))^2*cos(gamma(t))^2*sin(alpha(t))-1/2*R^4*h*rho*cos(alpha(t))*cos(beta(t))*sin(beta(t))*diff(psi(t),t)*cos(phi(t))*diff(theta(t),t)*sin(theta(t))*Pi-1/6*R^2*h^3*rho*cos(beta(t))*sin(alpha(t))*sin(beta(t))*diff(gamma(t),t)*sin(gamma(t))*diff(psi(t),t)*sin(theta(t))*Pi+1/6*R^2*h^3*rho*cos(beta(t))*sin(alpha(t))*sin(beta(t))*cos(gamma(t))*diff(psi(t),t)*diff(theta(t),t)*cos(theta(t))*Pi-1/6*R^2*h^3*rho*diff(gamma(t),t)*cos(gamma(t))*cos(beta(t))*sin(alpha(t))*sin(beta(t))*diff(theta(t),t)*cos(phi(t))*Pi-1/4*R^4*h*rho*diff(theta(t),t)^2*sin(beta(t))*Pi*sin(gamma(t))*cos(beta(t))*cos(alpha(t))*cos(phi(t))^2*cos(gamma(t))-1/12*R^2*h^3*rho*cos(beta(t))*cos(alpha(t))*sin(beta(t))*diff(psi(t),t)^2*sin(theta(t))*cos(theta(t))*sin(phi(t))*Pi+1/2*R^4*h*rho*Pi*cos(beta(t))^2*diff(gamma(t),t)*diff(psi(t),t)*cos(phi(t))*cos(alpha(t))*cos(theta(t))*sin(alpha(t))+1/12*R^2*h^3*rho*Pi*sin(gamma(t))*cos(beta(t))^2*sin(alpha(t))^2*diff(psi(t),t)^2*sin(phi(t))*cos(phi(t))*cos(theta(t))^2+1/12*R^2*h^3*rho*sin(gamma(t))*sin(alpha(t))^2*cos(beta(t))^2*diff(gamma(t),t)*diff(psi(t),t)*cos(theta(t))*sin(phi(t))*Pi+1/12*R^2*h^3*rho*sin(alpha(t))^2*cos(beta(t))^2*cos(gamma(t))*diff(psi(t),t)^2*sin(theta(t))*cos(phi(t))*cos(theta(t))*Pi+1/12*R^2*h^3*rho*sin(alpha(t))^2*cos(beta(t))^2*cos(gamma(t))*diff(psi(t),t)*cos(phi(t))*cos(theta(t))*diff(phi(t),t)*Pi+1/12*R^2*h^3*rho*sin(alpha(t))^2*cos(beta(t))^2*cos(gamma(t))*diff(psi(t),t)*diff(theta(t),t)*sin(theta(t))*sin(phi(t))*Pi-1/6*R^2*h^3*rho*sin(gamma(t))*sin(alpha(t))^2*cos(beta(t))^2*diff(psi(t),t)*diff(theta(t),t)*cos(phi(t))^2*cos(theta(t))*Pi+1/2*R^4*h*rho*diff(theta(t),t)*Pi*sin(gamma(t))*cos(beta(t))^2*sin(alpha(t))^2*diff(psi(t),t)*cos(phi(t))^2*cos(theta(t))-1/4*R^4*h*rho*sin(gamma(t))*cos(beta(t))^2*sin(alpha(t))*diff(beta(t),t)*diff(psi(t),t)*cos(theta(t))*sin(phi(t))*Pi+1/12*R^2*h^3*rho*sin(gamma(t))*cos(beta(t))^2*sin(alpha(t))*diff(beta(t),t)*diff(psi(t),t)*cos(theta(t))*sin(phi(t))*Pi+1/4*R^4*h*rho*sin(gamma(t))*diff(beta(t),t)*sin(beta(t))^2*sin(alpha(t))*diff(psi(t),t)*cos(theta(t))*sin(phi(t))*Pi+1/12*R^2*h^3*rho*sin(gamma(t))*cos(beta(t))*sin(alpha(t))*sin(beta(t))*diff(diff(psi(t),t),t)*cos(theta(t))*sin(phi(t))*Pi-1/12*R^2*h^3*rho*sin(gamma(t))*diff(beta(t),t)*sin(beta(t))^2*sin(alpha(t))*diff(psi(t),t)*cos(theta(t))*sin(phi(t))*Pi+1/6*R^2*h^3*rho*cos(beta(t))*cos(gamma(t))*diff(psi(t),t)*cos(theta(t))*sin(phi(t))*Pi*diff(beta(t),t)*sin(beta(t))-1/4*R^4*h*rho*sin(gamma(t))*cos(beta(t))*sin(alpha(t))*sin(beta(t))*diff(diff(psi(t),t),t)*cos(theta(t))*sin(phi(t))*Pi-1/2*R^4*h*rho*Pi*cos(beta(t))*diff(psi(t),t)*sin(phi(t))*cos(theta(t))*cos(gamma(t))*diff(beta(t),t)*sin(beta(t))-1/6*R^2*h^3*rho*Pi*cos(beta(t))^2*diff(psi(t),t)^2*cos(phi(t))^2*cos(alpha(t))*cos(theta(t))^2*sin(alpha(t))-1/6*R^2*h^3*rho*diff(theta(t),t)*Pi*cos(beta(t))^2*diff(gamma(t),t)*sin(phi(t))*cos(alpha(t))*sin(alpha(t))+1/12*R^2*h^3*rho*cos(beta(t))^2*sin(alpha(t))*diff(beta(t),t)*cos(gamma(t))*diff(psi(t),t)*sin(theta(t))*Pi-1/12*R^2*h^3*rho*sin(gamma(t))*cos(beta(t))^2*sin(alpha(t))*diff(beta(t),t)*diff(theta(t),t)*cos(phi(t))*Pi+1/2*R^4*h*rho*Pi*cos(beta(t))^2*diff(psi(t),t)^2*cos(alpha(t))*cos(theta(t))^2*cos(gamma(t))^2*sin(alpha(t))-1/12*R^2*h^3*rho*diff(theta(t),t)^2*Pi*cos(beta(t))^2*cos(phi(t))^2*cos(alpha(t))*cos(gamma(t))^2*sin(alpha(t))-1/6*R^2*h^3*rho*Pi*cos(beta(t))^2*diff(psi(t),t)^2*cos(alpha(t))*cos(theta(t))^2*cos(gamma(t))^2*sin(alpha(t))+1/4*R^4*h*rho*Pi*cos(beta(t))^2*diff(gamma(t),t)*sin(theta(t))*cos(alpha(t))^2*diff(psi(t),t)*cos(gamma(t))+1/4*R^4*h*rho*sin(gamma(t))*cos(beta(t))*cos(alpha(t))*sin(beta(t))*cos(gamma(t))*diff(psi(t),t)^2*Pi-1/12*R^2*h^3*rho*sin(gamma(t))*sin(alpha(t))^2*cos(beta(t))^2*diff(theta(t),t)^2*cos(phi(t))*sin(phi(t))*Pi+1/4*R^4*h*rho*diff(theta(t),t)*Pi*sin(gamma(t))*cos(beta(t))^2*diff(gamma(t),t)*sin(alpha(t))^2*cos(phi(t))+1/4*R^4*h*rho*diff(theta(t),t)^2*Pi*sin(gamma(t))*cos(beta(t))^2*sin(alpha(t))^2*sin(phi(t))*cos(phi(t))-1/4*R^4*h*rho*Pi*cos(beta(t))^2*diff(gamma(t),t)*sin(theta(t))*sin(alpha(t))^2*diff(psi(t),t)*cos(gamma(t))+1/2*R^4*h*rho*Pi*cos(beta(t))^2*diff(psi(t),t)^2*cos(phi(t))^2*cos(alpha(t))*cos(theta(t))^2*sin(alpha(t))+1/2*R^4*h*rho*diff(theta(t),t)*Pi*cos(beta(t))^2*diff(gamma(t),t)*sin(phi(t))*cos(alpha(t))*sin(alpha(t))+1/4*R^4*h*rho*diff(theta(t),t)^2*Pi*cos(beta(t))^2*cos(phi(t))^2*cos(alpha(t))*cos(gamma(t))^2*sin(alpha(t))-1/4*R^4*h*rho*diff(theta(t),t)^2*Pi*sin(gamma(t))*cos(beta(t))^2*cos(alpha(t))^2*sin(phi(t))*cos(phi(t))-1/4*R^4*h*rho*diff(theta(t),t)*Pi*sin(gamma(t))*cos(beta(t))^2*diff(gamma(t),t)*cos(alpha(t))^2*cos(phi(t))-1/12*R^2*h^3*rho*sin(gamma(t))*cos(beta(t))*cos(alpha(t))*sin(beta(t))*cos(gamma(t))*diff(psi(t),t)^2*Pi+1/4*R^4*h*rho*Pi*cos(beta(t))^2*diff(psi(t),t)*diff(phi(t),t)*cos(phi(t))*cos(theta(t))*cos(gamma(t))-1/4*R^4*h*rho*Pi*cos(beta(t))^2*diff(psi(t),t)*sin(phi(t))*diff(theta(t),t)*sin(theta(t))*cos(gamma(t))-1/4*R^4*h*rho*Pi*cos(beta(t))^2*diff(psi(t),t)*sin(phi(t))*cos(theta(t))*diff(gamma(t),t)*sin(gamma(t))+1/2*R^4*h*rho*cos(beta(t))*sin(alpha(t))*sin(beta(t))*diff(gamma(t),t)*sin(gamma(t))*diff(phi(t),t)*Pi-1/6*R^2*h^3*rho*cos(beta(t))*sin(alpha(t))*sin(beta(t))*diff(gamma(t),t)*sin(gamma(t))*diff(phi(t),t)*Pi-1/4*R^4*h*rho*cos(beta(t))*sin(alpha(t))*sin(beta(t))*cos(gamma(t))*diff(diff(psi(t),t),t)*sin(theta(t))*Pi+1/4*R^4*h*rho*diff(beta(t),t)*sin(beta(t))^2*sin(alpha(t))*cos(gamma(t))*diff(psi(t),t)*sin(theta(t))*Pi-1/12*R^2*h^3*rho*cos(alpha(t))*cos(beta(t))*sin(beta(t))*diff(diff(psi(t),t),t)*cos(phi(t))*cos(theta(t))*Pi+1/12*R^2*h^3*rho*cos(alpha(t))*diff(beta(t),t)*sin(beta(t))^2*diff(psi(t),t)*cos(phi(t))*cos(theta(t))*Pi+1/4*R^4*h*rho*sin(gamma(t))*cos(beta(t))*sin(alpha(t))*sin(beta(t))*diff(diff(theta(t),t),t)*cos(phi(t))*Pi-1/4*R^4*h*rho*sin(gamma(t))*diff(beta(t),t)*sin(beta(t))^2*sin(alpha(t))*diff(theta(t),t)*cos(phi(t))*Pi-1/4*R^4*h*rho*cos(alpha(t))*diff(beta(t),t)*sin(beta(t))^2*diff(psi(t),t)*cos(phi(t))*cos(theta(t))*Pi+1/12*R^2*h^3*rho*cos(beta(t))*sin(alpha(t))*sin(beta(t))*cos(gamma(t))*diff(diff(psi(t),t),t)*sin(theta(t))*Pi-1/12*R^2*h^3*rho*diff(beta(t),t)*sin(beta(t))^2*sin(alpha(t))*cos(gamma(t))*diff(psi(t),t)*sin(theta(t))*Pi-1/12*R^2*h^3*rho*sin(gamma(t))*cos(beta(t))*sin(alpha(t))*sin(beta(t))*diff(diff(theta(t),t),t)*cos(phi(t))*Pi+1/12*R^2*h^3*rho*sin(gamma(t))*diff(beta(t),t)*sin(beta(t))^2*sin(alpha(t))*diff(theta(t),t)*cos(phi(t))*Pi+1/2*R^4*h*rho*Pi*sin(gamma(t))*cos(beta(t))*sin(theta(t))*diff(psi(t),t)*diff(beta(t),t)*sin(beta(t))+1/2*R^4*h*rho*diff(theta(t),t)*Pi*cos(beta(t))*cos(phi(t))*cos(gamma(t))*diff(beta(t),t)*sin(beta(t))-1/6*R^2*h^3*rho*sin(gamma(t))*cos(beta(t))*diff(psi(t),t)*sin(theta(t))*Pi*diff(beta(t),t)*sin(beta(t))+1/4*R^4*h*rho*cos(alpha(t))*cos(beta(t))*sin(beta(t))*diff(diff(psi(t),t),t)*cos(phi(t))*cos(theta(t))*Pi+1/4*R^4*h*rho*cos(alpha(t))*cos(beta(t))^2*diff(beta(t),t)*diff(psi(t),t)*cos(phi(t))*cos(theta(t))*Pi-1/6*R^2*h^3*rho*sin(gamma(t))*cos(beta(t))*cos(alpha(t))*sin(beta(t))*cos(gamma(t))*diff(psi(t),t)*diff(theta(t),t)*cos(phi(t))*cos(theta(t))*sin(phi(t))*Pi+1/2*R^4*h*rho*sin(gamma(t))*cos(beta(t))*cos(alpha(t))*sin(beta(t))*cos(gamma(t))*diff(psi(t),t)*diff(theta(t),t)*cos(phi(t))*cos(theta(t))*sin(phi(t))*Pi+1/12*R^2*h^3*rho*cos(beta(t))^2*sin(alpha(t))*diff(beta(t),t)*cos(gamma(t))*diff(phi(t),t)*Pi+1/4*R^4*h*rho*diff(theta(t),t)*Pi*cos(beta(t))^2*diff(phi(t),t)*sin(phi(t))*cos(gamma(t))+1/4*R^4*h*rho*diff(theta(t),t)*Pi*cos(beta(t))^2*cos(phi(t))*diff(gamma(t),t)*sin(gamma(t))+1/12*R^2*h^3*rho*diff(gamma(t),t)*cos(gamma(t))*cos(beta(t))^2*diff(psi(t),t)*sin(theta(t))*Pi+1/12*R^2*h^3*rho*sin(gamma(t))*cos(beta(t))^2*diff(psi(t),t)*diff(theta(t),t)*cos(theta(t))*Pi-1/12*R^2*h^3*rho*cos(beta(t))^2*cos(gamma(t))*diff(theta(t),t)*diff(phi(t),t)*sin(phi(t))*Pi+1/2*R^4*h*rho*diff(gamma(t),t)*sin(gamma(t))*diff(psi(t),t)*cos(theta(t))*sin(phi(t))*Pi+1/2*R^4*h*rho*cos(gamma(t))*diff(psi(t),t)*diff(theta(t),t)*sin(theta(t))*sin(phi(t))*Pi-1/2*R^4*h*rho*cos(gamma(t))*diff(psi(t),t)*cos(theta(t))*diff(phi(t),t)*cos(phi(t))*Pi-1/4*R^4*h*rho*Pi*diff(gamma(t),t)*cos(gamma(t))*cos(beta(t))^2*sin(theta(t))*diff(psi(t),t)-1/4*R^4*h*rho*Pi*sin(gamma(t))*cos(beta(t))^2*diff(theta(t),t)*cos(theta(t))*diff(psi(t),t)-1/12*R^2*h^3*rho*sin(gamma(t))*sin(alpha(t))*diff(beta(t),t)*diff(theta(t),t)*cos(phi(t))*Pi-1/4*R^4*h*rho*Pi*cos(beta(t))^2*diff(phi(t),t)^2*cos(alpha(t))*cos(gamma(t))^2*sin(alpha(t))-1/2*R^4*h*rho*diff(theta(t),t)^2*Pi*cos(beta(t))^2*cos(phi(t))^2*cos(alpha(t))*sin(alpha(t))+1/12*R^2*h^3*rho*Pi*cos(beta(t))^2*diff(psi(t),t)^2*cos(alpha(t))*cos(theta(t))^2*sin(alpha(t))+1/12*R^2*h^3*rho*Pi*cos(beta(t))^2*diff(psi(t),t)^2*cos(alpha(t))*cos(gamma(t))^2*sin(alpha(t))-1/4*R^4*h*rho*Pi*cos(beta(t))^2*diff(psi(t),t)^2*cos(alpha(t))*cos(theta(t))^2*sin(alpha(t))-1/4*R^4*h*rho*Pi*cos(beta(t))^2*diff(psi(t),t)^2*cos(alpha(t))*cos(gamma(t))^2*sin(alpha(t))+1/12*R^2*h^3*rho*Pi*cos(beta(t))^2*diff(phi(t),t)^2*cos(alpha(t))*cos(gamma(t))^2*sin(alpha(t))+1/6*R^2*h^3*rho*cos(alpha(t))*cos(beta(t))^2*diff(theta(t),t)^2*cos(phi(t))^2*Pi*sin(alpha(t))+1/4*R^4*h*rho*Pi*cos(beta(t))^2*diff(gamma(t),t)*cos(alpha(t))^2*diff(phi(t),t)*cos(gamma(t))-1/12*R^2*h^3*rho*cos(alpha(t))*diff(beta(t),t)*diff(psi(t),t)*cos(phi(t))*cos(theta(t))*Pi+1/12*R^2*h^3*rho*sin(alpha(t))^2*cos(beta(t))^2*diff(gamma(t),t)*cos(gamma(t))*diff(phi(t),t)*Pi-1/4*R^4*h*rho*Pi*cos(beta(t))^2*diff(gamma(t),t)*sin(alpha(t))^2*diff(phi(t),t)*cos(gamma(t))-1/12*R^2*h^3*rho*cos(alpha(t))^2*cos(beta(t))^2*diff(gamma(t),t)*cos(gamma(t))*diff(phi(t),t)*Pi-1/4*R^4*h*rho*sin(gamma(t))*sin(alpha(t))*diff(beta(t),t)*diff(theta(t),t)*cos(phi(t))*Pi+1/4*R^4*h*rho*sin(alpha(t))*diff(beta(t),t)*cos(gamma(t))*diff(psi(t),t)*sin(theta(t))*Pi-1/4*R^4*h*rho*cos(alpha(t))*diff(beta(t),t)*diff(psi(t),t)*cos(phi(t))*cos(theta(t))*Pi+1/12*R^2*h^3*rho*sin(alpha(t))*diff(beta(t),t)*cos(gamma(t))*diff(psi(t),t)*sin(theta(t))*Pi+1/4*R^4*h*rho*Pi*cos(beta(t))^2*diff(diff(psi(t),t),t)*sin(phi(t))*cos(theta(t))*cos(gamma(t))-1/3*R^2*h^3*rho*cos(alpha(t))*cos(beta(t))^2*diff(psi(t),t)*diff(theta(t),t)*cos(phi(t))*cos(theta(t))*sin(phi(t))*Pi*sin(alpha(t))+1/4*R^4*h*rho*sin(beta(t))*Pi*sin(gamma(t))*cos(beta(t))*cos(alpha(t))*diff(psi(t),t)^2*cos(phi(t))^2*cos(theta(t))^2*cos(gamma(t))-1/4*R^4*h*rho*sin(beta(t))*Pi*cos(beta(t))*diff(psi(t),t)^2*sin(phi(t))*cos(phi(t))*sin(alpha(t))*cos(theta(t))^2*cos(gamma(t))+1/6*R^2*h^3*rho*cos(beta(t))*cos(alpha(t))*sin(beta(t))*cos(gamma(t))^2*diff(psi(t),t)^2*sin(theta(t))*cos(theta(t))*sin(phi(t))*Pi+R^4*h*rho*diff(theta(t),t)*Pi*cos(beta(t))^2*diff(psi(t),t)*sin(phi(t))*cos(phi(t))*cos(alpha(t))*cos(theta(t))*sin(alpha(t))+1/6*R^2*h^3*rho*Pi*sin(gamma(t))*cos(beta(t))^2*sin(theta(t))*diff(psi(t),t)^2*sin(phi(t))*cos(alpha(t))*cos(theta(t))*cos(gamma(t))*sin(alpha(t))-1/2*R^4*h*rho*diff(theta(t),t)*Pi*cos(beta(t))^2*diff(psi(t),t)*sin(phi(t))*cos(phi(t))*cos(alpha(t))*cos(theta(t))*cos(gamma(t))^2*sin(alpha(t))-1/2*R^4*h*rho*Pi*sin(gamma(t))*cos(beta(t))^2*sin(theta(t))*diff(psi(t),t)^2*sin(phi(t))*cos(alpha(t))*cos(theta(t))*cos(gamma(t))*sin(alpha(t))-1/2*R^4*h*rho*Pi*sin(gamma(t))*cos(beta(t))^2*diff(psi(t),t)*sin(phi(t))*diff(phi(t),t)*cos(alpha(t))*cos(theta(t))*cos(gamma(t))*sin(alpha(t))+1/2*R^4*h*rho*diff(theta(t),t)*Pi*sin(gamma(t))*cos(beta(t))^2*sin(theta(t))*diff(psi(t),t)*cos(phi(t))*cos(alpha(t))*cos(gamma(t))*sin(alpha(t))-1/6*R^2*h^3*rho*diff(theta(t),t)*Pi*sin(gamma(t))*cos(beta(t))^2*sin(theta(t))*diff(psi(t),t)*cos(phi(t))*cos(alpha(t))*cos(gamma(t))*sin(alpha(t))+1/6*R^2*h^3*rho*Pi*sin(gamma(t))*cos(beta(t))^2*diff(psi(t),t)*sin(phi(t))*diff(phi(t),t)*cos(alpha(t))*cos(theta(t))*cos(gamma(t))*sin(alpha(t))+1/6*R^2*h^3*rho*diff(theta(t),t)*Pi*cos(beta(t))^2*diff(psi(t),t)*sin(phi(t))*cos(phi(t))*cos(alpha(t))*cos(theta(t))*cos(gamma(t))^2*sin(alpha(t)) = 0"