Hi,

I have the following problem with plots in Maple 2020: I wanted to create several plots which should have the exact same size. I used the command size = [400,400] and Maple is creating a plot with that size, but Maple applies a white frame to my plots which does not have the same size in all cases (see the red mark of different length in the pictures below), so that the effective size differs, which is extremely ugly if you want to arrange several pictures in a document. So my question is

1.) How to remove this white frame, so that the efftictive picture size is actually 400x400?

2.) If 1.) is not possible, how can I adjust the frame so that it always has the same size?

My Code:

with(plots);
with(ColorTools);
with(plottools);
plot1 := inequal(0 <= y^3 + x^2, x = -5 .. 5, y = -5 .. 5, filledregions, color = blue, background = "Gainsboro", size = [400, 400]);
plot2 := plot(x = -5 .. 5, y = -5 .. 5, background = "Gainsboro");
ll := line([-3, 0], [3, 0]), color = blue, thickness = 5;
l := line([-3, 2], [5, 5]), color = blue, thickness = 2;
display(plot2, l, size = [400, 400]);

Thanks for your help.

I can't figure how to make my own type, which is rational and greater than one.

There are buildin types for postive and posint, and so on. But what if I want to make for rational and greater than some value, say 1?

This is easy to do using patmatch, using the conditional. But do not know how to do it to make my own type.

Here is a simple example. I want to check for sin(x)^n, where n is rational and must be >1.  Using pathmatch

restart;
patmatch(sin(x)^(1/2),conditional(sin(x)^a::rational,a>1),'la');la;

does it. Using structured type, the best I could do is this

restart;

TypeTools:-AddType('my_sin',specfunc(sin)^And(rational,positive));
type(sin(x)^(1/2),'my_sin');

But this does not check for >1, only positive.

Any suggestions? I know I could do this using other means, by direct parsing, using op, and so on. But I'd like to learn how to do it using structured type, just to learn the syntax if there is one.

is it possible to use conditional with structured types? But need  name to do that, like with patmatch, but there is no such syntax in structured types. 

ps. I think conditional does not make much sense with structured type. But I need to figure how to make my own type, which is rational and say >1, or integer and say >2 and so on. I just do not know how to do that yet. But I am sure there is a way. Will try to figure it out.

Hi please  help me in this problem in maple 18 

How do I solve the system K=B and find values 

x_{0},y_{0},z{0}

I posed the problem in the form pdf and mw

thank you 

problem.mw

problem.pdf

I can't get pdsolve to solve this pde. Here are my tries below. One of them works, but the analytical solution Maple gives is wrong. So I am not sure if it needs some additional hints or some other help to make it give the correct solution.

This is Laplace pde but with non-zero on RHS. On disk centered at origin and with radius 1, with given BC on edge of disk.

restart;
the_rhs := r^2*cos(theta)*sin(theta);
pde := VectorCalculus:-Laplacian(u(r,theta),'polar'[r,theta])=the_rhs;
bc := u(1, theta) = cos(theta)*sin(theta);
sol:=pdsolve([pde, bc], u(r, theta));

Maple gives

But this does not look right. There is a complex exponential on its own there. The solution should be real. 

I tried to solve it using pdsolve numerically, but could not make it work. It kepts complaining about missing BC. I am not good with numerical solvers in Maple. May be someone could try to do it.

evalf(eval(rhs(sol),[r= sqrt(1/4 + 1/4),theta=Pi/4]))

gives 0.5951891582 which is wrong,. It should be  0.23958

It tried to give it the symmetry conditions on theta, but now it did not solve it

restart;
the_rhs := r^2*cos(theta)*sin(theta);
pde := VectorCalculus:-Laplacian(u(r,theta),'polar'[r,theta])=the_rhs;
bc := u(1, theta) = cos(theta)*sin(theta),u(r,0) = u(r,2*Pi),(D[2](u))(r, 0) = (D[2](u))(r, 2*Pi);
sol:=pdsolve([pde, bc], u(r, theta));

No solution. I tried adding HINT = boundedseries(r=0) but that also did not help. No solution.

Could someone solve this PDE numerically in Maple and check the solution at the above location? It should be 0.23958 which shows the analytical solution given is not correct. Can Maple solve this numerically?

Maple 2021.1

ps. Analytical solution was obtained by Mathematica. I did not solve this by hand.

Hello,

I have lists of points of different ranges that I want to plot on the same graph. I would not want to manually cut the list for the range I want to see (there are many lists) so I thought about the view option. But there seems to be a strange behaviour of pointplot combined with view:

restart;
with(plots);
pointplot([[0.1, 0.5], [0.7, 0.7]], view = [0 .. 1, 0 .. 0.1]);

Hello Everyone,
I have written a procedure which takes one input and gives one output. Now, if I want to find the root of that procedure (as in when my procedure gives a '0' output), Is there any method to do it directly?
One idea that I got is to use maple's inbuilt Secant method (with(Student[NumericalAnalysis]). But this method can be done only for functions. Hence, my question.
How can we convert a maple procedure which takes one input and gives one output into a maple function? Also, is there a method to apply secant method on such function?
Please consider the following example:

As can be seen, it seems not possible to apply secant method on procedures.
Can someone please help me in doing this or any alternate method to get the job done?


 

I would like to assign the values of basis vectors which are calculated by PlanePlot. Is there anyone who knows how to use them or how to imitate the command?

PlanePlot(z+tan(-Pi/2+alph||29)*(x*sin(-Pi/2+gam||29)-y*cos(-Pi/2+gam||29))=0, [x,y,z],showbasis):
normal vector: <-.2497, .7343e-1, 1.>
equation of plane: -.2497*x+.7343e-1*y+1.*z = 0.
point on plane nearest origin: <-0., 0., 0.>
basis vectors: <.7106e-1, .9959, -.5539e-1>, <.9678, -.5539e-1, .2457>

It looks like as below:

I was expanding an expressions, and found that Maple expands trig also. Then found there is a way to turn that off. Which is great. I tested it in my worksheet and it works. But when I put the same code inside my module, I get an error.  It seems like some scoping issue which I do not understand.

Here is a MWE

restart;
expr:=(1+x)*sin(3*x+k);
forget(expand);
expand(expandoff());
expandoff(cos,sin,exp);
F:=expand(expr);
forget(expand);

which gives 

Which does what I want. But inside a proc, the code fails

foo:=proc(expr)
local F;

forget(expand);
expand(expandoff());
expandoff(cos,sin,exp);
F:=expand(expr);
forget(expand);

end proc;

And now when calling it as 

expr:=x*sin(3*x+k);
foo(expr)

What Am I doing wrong? And how to make it work? I need to prevent expand() from expanding these functions. Help says that expand has memory table, and I want to clear that before and after doing this, so it does not affect later code.

Maple 2021.1

Update

I just found out, if I add a restart before defining my function, then the error goes away

restart;
expr:=(1+x)*sin(3*x+k);
forget(expand);
expand(expandoff());
expandoff(cos,sin,exp);
F:=expand(expr);
forget(expand);

restart;  # had to add this. But why??
foo:=proc(expr)
local F;

forget(expand);
expand(expandoff());
expandoff(cos,sin,exp);
F:=expand(expr);
forget(expand);
end proc;

#
expr:=(1+x)*sin(3*x+k);
foo(expr)

#no error. now it works

I need this to work repeatedly inside a module, so I can't do restart each time ofcourse like the above.

Even if I do restart before defining the function, next time I call it, it will fail:

restart; 
foo:=proc(expr)
local F;

forget(expand);
expand(expandoff());
expandoff(cos,sin,exp);
F:=expand(expr);
forget(expand);
end proc;

#
expr:=(1+x)*sin(3*x+k);
foo(expr); #OK
foo(expr); #FAILED

Very strange. 

We consider a fixed circle (C) tangent to a fixed line Δ at a given point O of this line. 
Circles Γ tangent to circles C in M and to the right Δ in N are studied.
 Show that the MN line passes through a fixed point I. Infer that the circles Γ remain orthogonal to a fixed circle.
My code is :
 restart; with(geometry);with(plots);
_EnvHorizontalName := 'x';_EnvVerticalName := 'y';
dist := proc (M, N) sqrt(Vdot(expand(M-N), expand(M-N))) end proc;
point(oo, 0, 3); p := 6;
point(N, 5, 0);
line(Delta, y = 0, [x, y]);
para := x^2 = 2*p*y;
solve(subs(x = 5, para), y); point(varpi, 5, 25/12);
line(alpha, [oo, varpi]); k := 3/(25/12);
point(M, (0+5*k)/(1+k), (3+25*k*(1/12))/(1+k)); 
circle(C, x^2+(y-3)^2 = 9, [x, y]);cir := implicitplot(x^2+(y-3)^2 = 9, x = -5 .. 5, y = -5 .. 7, color = blue);
Para := implicitplot(para, x = -40 .. 40, y = 0 .. 40, linestyle = 3, color = coral);
homothety(J, N, -k, M); coordinates(J);
circle(C1, (y-25/12)^2+(x-5)^2 = (25/12)^2, [x, y]);line(lNJ, [N, J]);
triangle(T1, [J, oo, M]); triangle(T2, [N, varpi, M]);
C1 := implicitplot((y-25/12)^2+(x-5)^2 = (25/12)^2, x = 2 .. 8, y = 0 .. 5, color = magenta);dr1 := draw([oo, Delta, varpi, N, M, J], printtext = true); dr2 := draw([alpha(color = black), lNJ(color = black), T1(color = green, filled = true), T2(color = green, filled = true)]);
inversion(M, M, C);
inversion(N, M, C);
Fig := proc (xOm)
local cir, c2, C2, C1, c3, C3, k, M, N, J, sol, dr, varpi;
global p, para, Para;
sol := solve(subs(x = xOm, para), y);
cir := (y-sol)^2+(x-xOm)^2 = sol^2; c2 := x^2+(y-3)^2 = 9;
geometry:-point(N, xOm, 0); sol := solve(subs(x = xOm, para), y);
geometry:-point(varpi, xOm, sol); k := 3/sol;
geometry:-point(M, xOm*k/(k+1), (3+k*sol)/(k+1));
geometry:-homothety(J, N, -k, M);
c3 := (x-(1/2)*xOm)^2+(y-3)^2 = (1/4)*dist(N, J)^2;
C1 := plots:-implicitplot(cir, x = -xOm .. 3*xOm, y = 0 .. 3*xOm, color = magenta);
C2 := plots:-implicitplot(c2, x = -xOm .. 2*xOm, y = 0 .. 2*xOm, color = blue);
C3 := plots:-implicitplot(c3, x = -xOm .. 2*xOm, y = 0 .. 2*xOm, color = blue);
dr := geometry:-draw([varpi, M, J]);
plots:-display([Para, C2, C1, C3, dr], view = [-xOm .. 3*xOm, -1 .. 3*xOm], axes = normal, scaling = constrained) end proc;

Fig(8);
display([seq(Fig(4+.8*i), i = 4 .. 15)]);
display({C1, Para, cir, dr1, dr2}, view = [-8 .. 8, -1 .. 8], axes = normal, scaling = constrained, size = [500, 500]);
I don't know what is that orthogonal circle to each tangent circles. Thank you to help me.

.

I can't find how to use Maple structured type, to check for say sin(x) and say sin(x)^2 without having to duplicate the code and make a structured type for each of the two cases. An example will explain.

I need to check if expression is of this form   m*sin(anything)^n*cos(anything)^r  Where in this example below, m,n,r can be integer or rational.     

The problem is that 'specfunc'(sin)^Or(integer,rational) will not match sin(x) but will only match if sin(x) is raised to an actual power different from default 1.

So I have to duplicate the check below, by writing Or('specfunc'(sin),'specfunc'(sin)^Or(integer,rational) and the same for cos(x). I am trying to avoid all this duplication, because as the structured type becomes more involved, it will hard to read and work with like this.

Is there a way to eliminate this? In Mathematica for example, this is done using the pattern   x_. where the little dot at the end there, allows for zero of more. So for the above, it will be  Sin[x_]^n_. and this will match Sin[anything] and also match Sin[anything]^2 

I know it is possible to use patmatch in Maple to do this, but is it possible without using pathmatch and just using structured types? as I am trying to do everything without having to use patmatch now.

restart;
my_type:=''`*`'( { Or('specfunc'(sin),'specfunc'(sin)^Or(integer,rational)),
                   Or('specfunc'(cos),'specfunc'(cos)^Or(integer,rational)),
                   Or(integer,rational)})';
type(3*sin(x)^2*cos(x)^3,my_type);
type(sin(x)^2*99*cos(x),my_type);
type(sin(x)*cos(x),my_type);
type(cos(x)*sin(x)^(1/2),my_type);

gives

btw, the help page for structured type mentiones zero or more occurances of but do not know how to do this for the above example

$-\Delta u=1$ in $D$

$\frac{\partial u}{\partial\nu}+(x_1^2+1)u=0$ on $\partial D$

where $D=\{(x_1,x_2) : x_1^2+x_2^2<1\}$.

restart;
line := x/100 - 1/2;
wave := cos(x / 5) * sin(x / 2); ## -1 <= wave <= 1
eq:= line - wave;
r1 := solve(line < -1); ## RealRange(-infinity,Open(-50))
r2 := solve(line > 1);  ## RealRange(Open(150),infinity)
## There are 27 solutions in RealRange(-50,150)

I searched Questions.  Found many solutions.  The only one that seems to work is Student[Calculus1]:-Roots()

Is there another way to do this?

I am trying to solve a system of polynomial equations (with rational number coefficients). For this I am computing its Gröbner basis using the F4 algorithm implemented in Maple. (`Groebner[Basis]`).

As far as I understand it, the algorithms solves the problem modulo one or more prime number and later reconstructs the full (rational number) solution. Usually it only takes a handfull of primes. But in my case it by now solved the problem successfully modulo about ~300 different primes, each time doing exactly the same computation (according to the log files running with `infolevel[GroebnerBasis]:=5`).

Could somebody enlighten me what this means? In particular, can I interpret this as indication that the system of equations does have or does not have a solution?

• My system of equations is rather large. Each individual solve modulo a prime takes half an hour on a powerful workstation pc. Other computeralgebra systems without Faugère's algorithms cant solve it at all.
• (part of) the output can be found here: https://gist.github.com/krox/484252f075eb19edd0ac865099a564ba . The curious thing to me is that each solve behaves exactly the same. So why is maple repeating it with different primes over and over again?

Hi there, I just wanna solve a nonlinear diff equation for a 2 variables function phi.

Boundary conditions :

1. at the wall over NACA 0012 : the derivatives wrt x & y are null
2. for y at infinity : the derivatives wrt x & y are null

How to deal with curved body unlike the flat plate ?

Thanks for your help in advance!

I want to compute the Lie derivative of a large expression in maple.
The expression is attached, there are 20 vars in total. But Maple cannot do that (after 7hrs of working!), but this is doable with another math tool which first turns the expression into subexpressions.
However, I want to do this in Maple.
I have already tried to simplify/combine/factor/use the Lie algebra package, none of them helped.
Any suggestion is highly appreciated

(6*m__1*(2*x[13]*x[12]*((x[12]*(cos(x[7])^2 - 2)*(x[17] - g__1)*cos(x[8])^2 - 2*x[12]*(x[17] - g__1)*cos(x[7])^2 + x[18]*cos(x[7])*sin(x[7])*sin(x[8]) + x[12]*(x[17] - g__1))*cos(x[6])^2 - 2*cos(x[8])*(x[12]*sin(x[8])*(x[17] - g__1)*sin(x[7]) + x[18]*cos(x[7])/2)*sin(x[6])*cos(x[6]) + x[12]*cos(x[8])^2*(x[17] - g__1))*cos(x[4])^5 + ((-x[12]*x[13]*x[18]*sin(x[4])*(cos(x[7])^2 - 2)*cos(x[8])^2 + (((x[17] - g__1)*x[12]^2 + (2*x[17] - 2*g__1)*x[13]^2 - x[19])*x[12]*sin(x[8])*cos(x[7])^2 - 3*(x[12]^2*x[18] - (2*x[20])/9)*sin(x[7])*cos(x[7])/2 - 2*((x[17] - g__1)*x[12]^2 + (2*x[17] - 2*g__1)*x[13]^2 - x[19])*x[12]*sin(x[8]))*cos(x[8]) + 4*(x[18]*cos(x[7])^2/2 + x[12]*sin(x[7])*sin(x[8])*(x[17] - g__1)*cos(x[7]) - x[18]/4)*sin(x[4])*x[13]*x[12])*cos(x[6])^2 + (2*((x[17] - g__1)*x[12]^2 + (2*x[17] - 2*g__1)*x[13]^2 - x[19])*sin(x[7])*x[12]*sin(x[6])*cos(x[8])^2 - 4*(x[12]*x[13]*sin(x[4])*sin(x[6])*(x[17] - g__1)*cos(x[7]) - x[13]*x[18]*sin(x[4])*sin(x[6])*sin(x[7])*sin(x[8])/2 - x[12]*x[15]*(x[17] - g__1)/24)*x[12]*cos(x[8]) - (3*(x[12]^2*x[18] - (2*x[20])/9)*sin(x[8])*cos(x[7])/2 + ((x[17] - g__1)*x[12]^2 + (2*x[17] - 2*g__1)*x[13]^2 - x[19])*sin(x[7])*x[12])*sin(x[6]))*cos(x[6]) + x[12]*(-x[13]*x[18]*cos(x[8])^2*sin(x[4]) + ((x[17] - g__1)*x[12]^2 + (2*x[17] - 2*g__1)*x[13]^2 - x[19])*sin(x[8])*cos(x[8]) + x[12]*sin(x[8])*(x[17] - g__1)*(sin(x[6])*sin(x[7])*x[15] + x[14]*cos(x[7]))/6))*cos(x[4])^4 + ((-5*x[12]^2*x[13]*(cos(x[7])^2 - 2)*(x[17] - g__1)*cos(x[8])^2 - sin(x[4])*(3*(x[12]^2*x[18] - (2*x[20])/9)*sin(x[8])*cos(x[7])^2/2 + ((x[17] - g__1)*x[12]^2 + (2*x[17] - 2*g__1)*x[13]^2 - x[19])*sin(x[7])*x[12]*cos(x[7]) - 3*(x[12]^2*x[18] - (2*x[20])/9)*sin(x[8]))*cos(x[8]) + 7*(x[12]*(x[17] - g__1)*cos(x[7])^2 - (2*x[18]*cos(x[7])*sin(x[7])*sin(x[8]))/7 - (5*x[12]*(x[17] - g__1))/7)*x[13]*x[12])*cos(x[6])^2 - (3*sin(x[4])*(x[12]^2*x[18] - (2*x[20])/9)*sin(x[7])*cos(x[8])^2 - 10*(x[18]*cos(x[7])/5 + x[12]*sin(x[8])*(x[17] - g__1)*sin(x[7]))*x[13]*x[12]*cos(x[8]) + sin(x[4])*(((x[17] - g__1)*x[12]^2 + (2*x[17] - 2*g__1)*x[13]^2 - x[19])*x[12]*sin(x[8])*cos(x[7]) - 3*(x[12]^2*x[18] - (2*x[20])/9)*sin(x[7])/2))*sin(x[6])*cos(x[6]) - 5*x[12]^2*x[13]*(x[17] - g__1)*cos(x[8])^2 - (3*sin(x[4])*(x[12]^2*x[18] - (2*x[20])/9)*sin(x[8])*cos(x[8]))/2 + (x[17] - g__1)*(x[15]*cos(x[7])*sin(x[4])*sin(x[6])/6 - x[14]*sin(x[4])*sin(x[7])/6 + x[13])*x[12]^2)*cos(x[4])^3 - (7*x[12]*((-(4*x[13]*x[18]*sin(x[4])*(cos(x[7])^2 - 2)*cos(x[8])^2)/7 + (((x[17] - g__1)*x[12]^2 + ((4*x[17])/7 - (4*g__1)/7)*x[13]^2 - (2*x[19])/7)*sin(x[8])*cos(x[7])^2 - (5*x[12]*x[18]*cos(x[7])*sin(x[7]))/7 - 2*((x[17] - g__1)*x[12]^2 + ((4*x[17])/7 - (4*g__1)/7)*x[13]^2 - (2*x[19])/7)*sin(x[8]))*cos(x[8]) + 10*(x[18]*cos(x[7])^2/5 + x[12]*sin(x[7])*sin(x[8])*(x[17] - g__1)*cos(x[7]) - (2*x[18])/5)*sin(x[4])*x[13]/7)*cos(x[6])^2 + (2*((x[17] - g__1)*x[12]^2 + ((4*x[17])/7 - (4*g__1)/7)*x[13]^2 - (2*x[19])/7)*sin(x[7])*sin(x[6])*cos(x[8])^2 + (-(10*x[12]*x[13]*sin(x[4])*sin(x[6])*(x[17] - g__1)*cos(x[7]))/7 + (8*x[13]*x[18]*sin(x[4])*sin(x[6])*sin(x[7])*sin(x[8]))/7 + x[12]*x[15]*(x[17] - g__1)/21)*cos(x[8]) - ((5*x[12]*x[18]*cos(x[7])*sin(x[8]))/7 + ((x[17] - g__1)*x[12]^2 + ((4*x[17])/7 - (4*g__1)/7)*x[13]^2 - (2*x[19])/7)*sin(x[7]))*sin(x[6]))*cos(x[6]) - (4*x[13]*x[18]*cos(x[8])^2*sin(x[4]))/7 + ((x[17] - g__1)*x[12]^2 + ((4*x[17])/7 - (4*g__1)/7)*x[13]^2 - (2*x[19])/7)*sin(x[8])*cos(x[8]) + x[12]*x[14]*sin(x[8])*(x[17] - g__1)*cos(x[7])/21 + (2*x[13]*x[18]*sin(x[4]))/7 + x[12]*x[15]*sin(x[6])*sin(x[7])*sin(x[8])*(x[17] - g__1)/21)*cos(x[4])^2)/2 + 3*(((4*x[13]*(cos(x[7])^2 - 2)*(x[17] - g__1)*cos(x[8])^2)/3 + sin(x[4])*((2*x[18]*cos(x[7])^2*sin(x[8]))/3 + x[12]*sin(x[7])*(x[17] - g__1)*cos(x[7]) - (4*x[18]*sin(x[8]))/3)*cos(x[8]) - (2*x[13]*(cos(x[7])^2 - 2)*(x[17] - g__1))/3)*cos(x[6])^2 + ((4*x[18]*cos(x[8])^2*sin(x[4])*sin(x[7]))/3 - (8*x[13]*sin(x[7])*sin(x[8])*(x[17] - g__1)*cos(x[8]))/3 + (x[12]*sin(x[8])*(x[17] - g__1)*cos(x[7]) - (2*x[18]*sin(x[7]))/3)*sin(x[4]))*sin(x[6])*cos(x[6]) + (4*x[13]*(x[17] - g__1)*cos(x[8])^2)/3 + (2*x[18]*cos(x[8])*sin(x[4])*sin(x[8]))/3 - (2*x[13]*(x[17] - g__1))/3)*x[12]^2*cos(x[4]) + 2*(x[17] - g__1)*(cos(x[8])*sin(x[8])*(cos(x[7])^2 - 2)*cos(x[6])^2 + 2*(cos(x[8])^2 - 1/2)*sin(x[7])*sin(x[6])*cos(x[6]) + cos(x[8])*sin(x[8]))*x[12]^3)*L*cos(x[5])^6 + (12*m__1*sin(x[5])*((-((x[12]^2*x[18] - x[20]/6)*(cos(x[7])^2 - 2)*cos(x[8])^2)/2 + (x[12]^2*x[18] - x[20]/6)*cos(x[7])^2 + sin(x[7])*((x[17] - g__1)*x[12]^2 + (x[17] - g__1)*x[13]^2 - x[19]/2)*x[12]*sin(x[8])*cos(x[7]) - x[12]^2*x[18]/2 + x[20]/12)*cos(x[6])^2 - (((x[17] - g__1)*x[12]^2 + (x[17] - g__1)*x[13]^2 - x[19]/2)*x[12]*cos(x[7]) - (x[12]^2*x[18] - x[20]/6)*sin(x[7])*sin(x[8]))*cos(x[8])*sin(x[6])*cos(x[6]) - cos(x[8])^2*(x[12]^2*x[18] - x[20]/6)/2)*L*cos(x[4])^5 - 6*m__1*sin(x[5])*((sin(x[4])*(cos(x[7])^2 - 2)*((x[17] - g__1)*x[12]^2 + (x[17] - g__1)*x[13]^2 - x[19]/2)*x[12]*cos(x[8])^2 + 3*x[13]*((2*x[18]*cos(x[7])^2*sin(x[8]))/3 + x[12]*sin(x[7])*(x[17] - g__1)*cos(x[7]) - (4*x[18]*sin(x[8]))/3)*x[12]*cos(x[8]) - 2*sin(x[4])*(((x[17] - g__1)*x[12]^2 + (x[17] - g__1)*x[13]^2 - x[19]/2)*x[12]*cos(x[7])^2 - (x[12]^2*x[18] - x[20]/6)*sin(x[7])*sin(x[8])*cos(x[7]) - (((x[17] - g__1)*x[12]^2 + (x[17] - g__1)*x[13]^2 - x[19]/2)*x[12])/2))*cos(x[6])^2 - 2*(-2*x[12]*x[13]*x[18]*cos(x[8])^2*sin(x[7]) + sin(x[4])*((x[12]^2*x[18] - x[20]/6)*cos(x[7]) + sin(x[7])*((x[17] - g__1)*x[12]^2 + (x[17] - g__1)*x[13]^2 - x[19]/2)*x[12]*sin(x[8]))*cos(x[8]) - 3*(x[12]*sin(x[8])*(x[17] - g__1)*cos(x[7]) - (2*x[18]*sin(x[7]))/3)*x[13]*x[12]/2)*sin(x[6])*cos(x[6]) + cos(x[8])*(sin(x[4])*((x[17] - g__1)*x[12]^2 + (x[17] - g__1)*x[13]^2 - x[19]/2)*cos(x[8]) + 2*x[13]*x[18]*sin(x[8]))*x[12])*L*cos(x[4])^4 + ((15*m__1*sin(x[5])*(cos(x[7])^2 - 2)*(x[12]^2*x[18] - (2*x[20])/15)*L*cos(x[8])^2 - 18*x[12]*((x[17] - g__1)*(L*x[12]*x[13]*sin(x[4])*sin(x[5])*m__1 - 2*(m__1 + m__p)*(x[17] - g__1)/9)*sin(x[8])*cos(x[7])^2 - (2*L*x[13]*x[18]*cos(x[7])*sin(x[4])*sin(x[5])*sin(x[7])*m__1)/3 - 2*(x[17] - g__1)*(L*x[12]*x[13]*sin(x[4])*sin(x[5])*m__1 - 2*(m__1 + m__p)*(x[17] - g__1)/9)*sin(x[8]))*cos(x[8]) - 33*m__1*sin(x[5])*(((7*x[12]^2*x[18])/11 - x[20]/33)*cos(x[7])^2 + ((x[17] - g__1)*x[12]^2 + ((4*x[17])/11 - (4*g__1)/11)*x[13]^2 - (2*x[19])/11)*sin(x[7])*x[12]*sin(x[8])*cos(x[7]) - (5*x[12]^2*x[18])/11 + (2*x[20])/33)*L)*cos(x[6])^2 + 33*(-12*(x[17] - g__1)*(L*x[12]*x[13]*sin(x[4])*sin(x[5])*m__1 - 2*(m__1 + m__p)*(x[17] - g__1)/9)*sin(x[7])*x[12]*cos(x[8])^2/11 + m__1*sin(x[5])*(((x[17] - g__1)*x[12]^2 + ((4*x[17])/11 - (4*g__1)/11)*x[13]^2 - (2*x[19])/11)*x[12]*cos(x[7]) - 10*(x[12]^2*x[18] - (2*x[20])/15)*sin(x[7])*sin(x[8])/11)*L*cos(x[8]) + (6*x[12]*((2*L*x[13]*x[18]*cos(x[7])*sin(x[4])*sin(x[5])*sin(x[8])*m__1)/3 + (x[17] - g__1)*(L*x[12]*x[13]*sin(x[4])*sin(x[5])*m__1 - 2*(m__1 + m__p)*(x[17] - g__1)/9)*sin(x[7])))/11)*sin(x[6])*cos(x[6]) + 15*m__1*sin(x[5])*(x[12]^2*x[18] - (2*x[20])/15)*L*cos(x[8])^2 - 18*(x[17] - g__1)*(L*x[12]*x[13]*sin(x[4])*sin(x[5])*m__1 - 2*(m__1 + m__p)*(x[17] - g__1)/9)*x[12]*sin(x[8])*cos(x[8]) - 3*(x[12]^2*x[18] - x[20]/3)*m__1*sin(x[5])*L)*cos(x[4])^3 + 15*x[12]*((m__1*sin(x[4])*sin(x[5])*(cos(x[7])^2 - 2)*L*((x[17] - g__1)*x[12]^2 + ((4*x[17])/5 - (4*g__1)/5)*x[13]^2 - (2*x[19])/5)*cos(x[8])^2 + ((4*L*x[13]*x[18]*cos(x[7])^2*sin(x[5])*sin(x[8])*m__1)/5 + 2*(x[17] - g__1)*sin(x[7])*(L*x[12]*x[13]*sin(x[5])*m__1 - (2*sin(x[4])*(m__1 + m__p)*(x[17] - g__1))/15)*cos(x[7]) - (8*L*x[13]*x[18]*sin(x[5])*sin(x[8])*m__1)/5)*cos(x[8]) - (9*m__1*sin(x[4])*(((x[17] - g__1)*x[12]^2 + ((2*x[17])/9 - (2*g__1)/9)*x[13]^2 - x[19]/9)*cos(x[7])^2 - (5*x[12]*x[18]*cos(x[7])*sin(x[7])*sin(x[8]))/9 + (-(5*x[17])/9 + (5*g__1)/9)*x[12]^2 + (-(4*x[17])/9 + (4*g__1)/9)*x[13]^2 + (2*x[19])/9)*sin(x[5])*L)/5)*cos(x[6])^2 + ((8*L*x[13]*x[18]*cos(x[8])^2*sin(x[5])*sin(x[6])*sin(x[7])*m__1)/5 - 2*(x[12]*x[18]*cos(x[7])/2 + sin(x[7])*((x[17] - g__1)*x[12]^2 + ((4*x[17])/5 - (4*g__1)/5)*x[13]^2 - (2*x[19])/5)*sin(x[8]))*m__1*sin(x[4])*sin(x[5])*L*sin(x[6])*cos(x[8]) + 2*(x[17] - g__1)*(L*x[12]*x[13]*sin(x[5])*m__1 - (2*sin(x[4])*(m__1 + m__p)*(x[17] - g__1))/15)*sin(x[8])*sin(x[6])*cos(x[7]) - m__1*(x[12]*x[15]*sin(x[8])*(x[17] - g__1)*sin(x[4]) + 12*x[13]*x[18]*sin(x[6])*sin(x[7]))*sin(x[5])*L/15)*cos(x[6]) + m__1*sin(x[5])*(sin(x[4])*((x[17] - g__1)*x[12]^2 + ((4*x[17])/5 - (4*g__1)/5)*x[13]^2 - (2*x[19])/5)*cos(x[8])^2 + (x[12]*x[14]*sin(x[4])*(x[17] - g__1)*cos(x[7])/15 + x[12]*x[15]*sin(x[6])*sin(x[7])*(x[17] - g__1)*sin(x[4])/15 + (4*x[13]*x[18]*sin(x[8]))/5)*cos(x[8]) - sin(x[4])*((x[17] - g__1)*x[12]^2 + (2*x[17] - 2*g__1)*x[13]^2 - x[19])/5)*L)*cos(x[4])^2 + 18*x[12]*((-(2*L*x[12]*x[18]*sin(x[5])*m__1*(cos(x[7])^2 - 2)*cos(x[8])^2)/3 + 4*(x[17] - g__1)*(cos(x[7])^2 - 2)*sin(x[8])*(L*x[12]*x[13]*sin(x[4])*sin(x[5])*m__1 - ((m__1 + m__p)*(x[17] - g__1))/6)*cos(x[8])/3 + m__1*sin(x[5])*x[12]*L*(x[18]*cos(x[7])^2/3 + x[12]*sin(x[7])*sin(x[8])*(x[17] - g__1)*cos(x[7]) - (2*x[18])/3))*cos(x[6])^2 - (-8*(x[17] - g__1)*sin(x[7])*(L*x[12]*x[13]*sin(x[4])*sin(x[5])*m__1 - ((m__1 + m__p)*(x[17] - g__1))/6)*cos(x[8])^2/3 + m__1*sin(x[5])*(x[12]*cos(x[7])*(x[17] - g__1) - (4*x[18]*sin(x[7])*sin(x[8]))/3)*x[12]*L*cos(x[8]) + 4*(x[17] - g__1)*sin(x[7])*(L*x[12]*x[13]*sin(x[4])*sin(x[5])*m__1 - ((m__1 + m__p)*(x[17] - g__1))/6)/3)*sin(x[6])*cos(x[6]) - (2*L*x[12]*x[18]*cos(x[8])^2*sin(x[5])*m__1)/3 + 4*(x[17] - g__1)*sin(x[8])*(L*x[12]*x[13]*sin(x[4])*sin(x[5])*m__1 - ((m__1 + m__p)*(x[17] - g__1))/6)*cos(x[8])/3 + L*x[12]*x[18]*sin(x[5])*m__1/3)*cos(x[4]) - 12*(x[17] - g__1)*m__1*sin(x[4])*sin(x[5])*x[12]^3*L*((cos(x[8])^2 - 1/2)*(cos(x[7])^2 - 2)*cos(x[6])^2 - 2*cos(x[6])*sin(x[7])*cos(x[8])*sin(x[6])*sin(x[8]) + cos(x[8])^2 - 1/2))*cos(x[5])^5 + (-18*m__1*x[13]*L*(((cos(x[7])^2 - 2)*((x[17] - g__1)*x[12]^2 + (x[17]/3 - g__1/3)*x[13]^2 - x[19]/6)*cos(x[8])^2 + ((-2*x[17] + 2*g__1)*x[12]^2 + (-(2*x[17])/3 + (2*g__1)/3)*x[13]^2 + x[19]/3)*cos(x[7])^2 + (4*x[12]*x[18]*cos(x[7])*sin(x[7])*sin(x[8]))/3 + (x[17] - g__1)*x[12]^2 + (x[17]/3 - g__1/3)*x[13]^2 - x[19]/6)*cos(x[6])^2 - 2*((2*x[12]*x[18]*cos(x[7]))/3 + sin(x[7])*sin(x[8])*((x[17] - g__1)*x[12]^2 + (x[17]/3 - g__1/3)*x[13]^2 - x[19]/6))*cos(x[8])*sin(x[6])*cos(x[6]) + cos(x[8])^2*((x[17] - g__1)*x[12]^2 + (x[17]/3 - g__1/3)*x[13]^2 - x[19]/6))*cos(x[4])^5 + ((12*L*sin(x[4])*m__1*x[12]*x[13]*x[18]*(cos(x[7])^2 - 2)*cos(x[8])^2 - 12*m__1*L*(((x[17] - g__1)*x[12]^2 + ((9*x[17])/2 - (9*g__1)/2)*x[13]^2 - x[19]/2)*x[12]*sin(x[8])*cos(x[7])^2 - sin(x[7])*(x[12]^2*x[18] + 5/4*x[13]^2*x[18] - 1/6*x[20])*cos(x[7]) - 2*((x[17] - g__1)*x[12]^2 + ((9*x[17])/2 - (9*g__1)/2)*x[13]^2 - x[19]/2)*x[12]*sin(x[8]))*cos(x[8]) - 24*L*x[12]*x[13]*x[18]*cos(x[7])^2*sin(x[4])*m__1 - 36*sin(x[7])*(-x[12]*(x[17] - g__1)^2*(m__1 + m__p)*sin(x[5])/9 + m__1*sin(x[4])*x[13]*L*((x[17] - g__1)*x[12]^2 + (x[17]/3 - g__1/3)*x[13]^2 - x[19]/6))*sin(x[8])*cos(x[7]) + 12*L*sin(x[4])*m__1*x[12]*x[13]*x[18])*cos(x[6])^2 + (-24*m__1*((x[17] - g__1)*x[12]^2 + ((9*x[17])/2 - (9*g__1)/2)*x[13]^2 - x[19]/2)*sin(x[7])*x[12]*L*sin(x[6])*cos(x[8])^2 + (36*(-x[12]*(x[17] - g__1)^2*(m__1 + m__p)*sin(x[5])/9 + m__1*sin(x[4])*x[13]*L*((x[17] - g__1)*x[12]^2 + (x[17]/3 - g__1/3)*x[13]^2 - x[19]/6))*sin(x[6])*cos(x[7]) - 24*m__1*(x[12]*x[18]*sin(x[4])*sin(x[6])*sin(x[7])*sin(x[8]) - x[13]*x[15]*(x[17] - g__1)/24)*x[13]*L)*cos(x[8]) + 12*m__1*(sin(x[8])*(x[12]^2*x[18] + 5/4*x[13]^2*x[18] - 1/6*x[20])*cos(x[7]) + ((x[17] - g__1)*x[12]^2 + ((9*x[17])/2 - (9*g__1)/2)*x[13]^2 - x[19]/2)*sin(x[7])*x[12])*L*sin(x[6]))*cos(x[6]) - 12*m__1*(-x[12]*x[13]*x[18]*cos(x[8])^2*sin(x[4]) + ((x[17] - g__1)*x[12]^2 + ((9*x[17])/2 - (9*g__1)/2)*x[13]^2 - x[19]/2)*x[12]*sin(x[8])*cos(x[8]) - x[13]^2*sin(x[8])*(x[17] - g__1)*(sin(x[6])*sin(x[7])*x[15] + x[14]*cos(x[7]))/12)*L)*cos(x[4])^4 + ((63*(cos(x[7])^2 - 2)*(-(2*x[12]*sin(x[4])*(x[17] - g__1)^2*(m__1 + m__p)*sin(x[5]))/63 + m__1*x[13]*L*((x[17] - g__1)*x[12]^2 + ((2*x[17])/7 - (2*g__1)/7)*x[13]^2 - (2*x[19])/21))*cos(x[8])^2 + 12*m__1*sin(x[4])*(sin(x[8])*(x[12]^2*x[18] + 5/4*x[13]^2*x[18] - 1/6*x[20])*cos(x[7])^2 + ((x[17] - g__1)*x[12]^2 + ((9*x[17])/2 - (9*g__1)/2)*x[13]^2 - x[19]/2)*sin(x[7])*x[12]*cos(x[7]) - 2*sin(x[8])*(x[12]^2*x[18] + 5/4*x[13]^2*x[18] - 1/6*x[20]))*L*cos(x[8]) + (4*x[12]*sin(x[4])*(x[17] - g__1)^2*(m__1 + m__p)*sin(x[5]) - 63*m__1*((x[17] - g__1)*x[12]^2 - x[19]/21)*x[13]*L)*cos(x[7])^2 + 24*L*x[12]*x[13]*x[18]*cos(x[7])*sin(x[7])*sin(x[8])*m__1 - 2*x[12]*sin(x[4])*(x[17] - g__1)^2*(m__1 + m__p)*sin(x[5]) + 63*m__1*x[13]*L*((x[17] - g__1)*x[12]^2 + ((2*x[17])/7 - (2*g__1)/7)*x[13]^2 - (2*x[19])/21))*cos(x[6])^2 + 12*(2*m__1*sin(x[4])*sin(x[7])*L*(x[12]^2*x[18] + 5/4*x[13]^2*x[18] - 1/6*x[20])*cos(x[8])^2 + (-2*L*x[12]*x[13]*x[18]*cos(x[7])*m__1 - (21*sin(x[7])*(-(2*x[12]*sin(x[4])*(x[17] - g__1)^2*(m__1 + m__p)*sin(x[5]))/63 + m__1*x[13]*L*((x[17] - g__1)*x[12]^2 + ((2*x[17])/7 - (2*g__1)/7)*x[13]^2 - (2*x[19])/21))*sin(x[8]))/2)*cos(x[8]) + m__1*sin(x[4])*(((x[17] - g__1)*x[12]^2 + ((9*x[17])/2 - (9*g__1)/2)*x[13]^2 - x[19]/2)*x[12]*sin(x[8])*cos(x[7]) - sin(x[7])*(x[12]^2*x[18] + 5/4*x[13]^2*x[18] - 1/6*x[20]))*L)*sin(x[6])*cos(x[6]) + (-2*x[12]*sin(x[4])*(x[17] - g__1)^2*(m__1 + m__p)*sin(x[5]) + 63*m__1*x[13]*L*((x[17] - g__1)*x[12]^2 + ((2*x[17])/7 - (2*g__1)/7)*x[13]^2 - (2*x[19])/21))*cos(x[8])^2 + 12*m__1*sin(x[4])*L*(x[12]^2*x[18] + 5/4*x[13]^2*x[18] - 1/6*x[20])*sin(x[8])*cos(x[8]) - 27*m__1*x[13]*(-x[13]*x[15]*sin(x[4])*sin(x[6])*(x[17] - g__1)*cos(x[7])/27 + x[13]*x[14]*sin(x[7])*(x[17] - g__1)*sin(x[4])/27 + (x[17] - g__1)*x[12]^2 + ((4*x[17])/9 - (4*g__1)/9)*x[13]^2 - x[19]/9)*L)*cos(x[4])^3 + ((-30*L*sin(x[4])*m__1*x[12]*x[13]*x[18]*(cos(x[7])^2 - 2)*cos(x[8])^2 + 24*m__1*(((x[17] - g__1)*x[12]^2 + ((17*x[17])/8 - (17*g__1)/8)*x[13]^2 - x[19]/4)*x[12]*sin(x[8])*cos(x[7])^2 - (5*sin(x[7])*x[18]*(x[12]^2 - x[13]^2/5)*cos(x[7]))/8 - 2*((x[17] - g__1)*x[12]^2 + ((17*x[17])/8 - (17*g__1)/8)*x[13]^2 - x[19]/4)*x[12]*sin(x[8]))*L*cos(x[8]) + 12*L*x[12]*x[13]*x[18]*cos(x[7])^2*sin(x[4])*m__1 + 45*(x[17] - g__1)*sin(x[7])*(-(4*x[12]*(m__1 + m__p)*(x[17] - g__1)*sin(x[5]))/45 + L*x[13]*sin(x[4])*m__1*(x[12]^2 - (2*x[13]^2)/15))*sin(x[8])*cos(x[7]) - 30*L*sin(x[4])*m__1*x[12]*x[13]*x[18])*cos(x[6])^2 + (48*m__1*sin(x[7])*((x[17] - g__1)*x[12]^2 + ((17*x[17])/8 - (17*g__1)/8)*x[13]^2 - x[19]/4)*x[12]*L*sin(x[6])*cos(x[8])^2 + (-45*(x[17] - g__1)*(-(4*x[12]*(m__1 + m__p)*(x[17] - g__1)*sin(x[5]))/45 + L*x[13]*sin(x[4])*m__1*(x[12]^2 - (2*x[13]^2)/15))*sin(x[6])*cos(x[7]) + m__1*(60*x[12]*x[13]*x[18]*sin(x[4])*sin(x[6])*sin(x[7])*sin(x[8]) + x[15]*(x[12] - x[13])*(x[12] + x[13])*(x[17] - g__1))*L)*cos(x[8]) - 24*m__1*((5*x[18]*sin(x[8])*(x[12]^2 - x[13]^2/5)*cos(x[7]))/8 + sin(x[7])*((x[17] - g__1)*x[12]^2 + ((17*x[17])/8 - (17*g__1)/8)*x[13]^2 - x[19]/4)*x[12])*L*sin(x[6]))*cos(x[6]) + 24*m__1*(-(5*x[12]*x[13]*x[18]*cos(x[8])^2*sin(x[4]))/4 + ((x[17] - g__1)*x[12]^2 + ((17*x[17])/8 - (17*g__1)/8)*x[13]^2 - x[19]/4)*x[12]*sin(x[8])*cos(x[8]) + x[14]*sin(x[8])*(x[12] - x[13])*(x[12] + x[13])*(x[17] - g__1)*cos(x[7])/24 + (3*x[12]*x[13]*x[18]*sin(x[4]))/4 + x[15]*sin(x[6])*sin(x[7])*sin(x[8])*(x[12] - x[13])*(x[12] + x[13])*(x[17] - g__1)/24)*L)*cos(x[4])^2 - 18*x[12]*((8*(x[17] - g__1)*(cos(x[7])^2 - 2)*(-sin(x[4])*(m__1 + m__p)*(x[17] - g__1)*sin(x[5])/12 + L*x[12]*x[13]*m__1)*cos(x[8])^2/3 + m__1*sin(x[4])*((2*x[12]*x[18]*cos(x[7])^2*sin(x[8]))/3 + (x[17] - g__1)*(x[12]^2 - x[13]^2/3)*sin(x[7])*cos(x[7]) - (4*x[12]*x[18]*sin(x[8]))/3)*L*cos(x[8]) - (x[17] - g__1)*((-sin(x[4])*(m__1 + m__p)*(x[17] - g__1)*sin(x[5])/9 + L*x[12]*x[13]*m__1)*cos(x[7])^2 + (2*sin(x[4])*(m__1 + m__p)*(x[17] - g__1)*sin(x[5]))/9 - (8*L*x[12]*x[13]*m__1)/3))*cos(x[6])^2 + ((4*L*x[12]*x[18]*cos(x[8])^2*sin(x[4])*sin(x[7])*m__1)/3 - 16*(x[17] - g__1)*sin(x[7])*sin(x[8])*(-sin(x[4])*(m__1 + m__p)*(x[17] - g__1)*sin(x[5])/12 + L*x[12]*x[13]*m__1)*cos(x[8])/3 + m__1*sin(x[4])*((x[17] - g__1)*(x[12]^2 - x[13]^2/3)*sin(x[8])*cos(x[7]) - (2*x[12]*x[18]*sin(x[7]))/3)*L)*sin(x[6])*cos(x[6]) + 8*(x[17] - g__1)*(-sin(x[4])*(m__1 + m__p)*(x[17] - g__1)*sin(x[5])/12 + L*x[12]*x[13]*m__1)*cos(x[8])^2/3 + (2*L*x[12]*x[18]*cos(x[8])*sin(x[4])*sin(x[8])*m__1)/3 - 5*(x[17] - g__1)*(-sin(x[4])*(m__1 + m__p)*(x[17] - g__1)*sin(x[5])/15 + L*x[12]*x[13]*m__1)/3)*cos(x[4]) - 12*(x[17] - g__1)*m__1*(cos(x[8])*sin(x[8])*(cos(x[7])^2 - 2)*cos(x[6])^2 + 2*(cos(x[8])^2 - 1/2)*sin(x[7])*sin(x[6])*cos(x[6]) + cos(x[8])*sin(x[8]))*x[12]^3*L)*cos(x[5])^4 + (-36*m__1*sin(x[5])*x[13]^2*((-x[18]*(cos(x[7])^2 - 2)*cos(x[8])^2/6 + x[18]*cos(x[7])^2/3 + x[12]*sin(x[7])*sin(x[8])*(x[17] - g__1)*cos(x[7]) - x[18]/6)*cos(x[6])^2 - cos(x[8])*(x[12]*cos(x[7])*(x[17] - g__1) - x[18]*sin(x[7])*sin(x[8])/3)*sin(x[6])*cos(x[6]) - x[18]*cos(x[8])^2/6)*L*cos(x[4])^5 + 36*((((x[17] - g__1)*(cos(x[7])^2 - 2)*(L*x[12]*x[13]*sin(x[4])*sin(x[5])*m__1 - ((m__1 + m__p)*(x[17] - g__1))/9)*cos(x[8])^2)/2 + m__1*sin(x[5])*((2*x[12]*x[18]*cos(x[7])^2*sin(x[8]))/3 + sin(x[7])*((x[17] - g__1)*x[12]^2 + (x[17]/2 - g__1/2)*x[13]^2 - x[19]/6)*cos(x[7]) - (4*x[12]*x[18]*sin(x[8]))/3)*L*cos(x[8]) - (x[17] - g__1)*(L*x[12]*x[13]*sin(x[4])*sin(x[5])*m__1 - ((m__1 + m__p)*(x[17] - g__1))/9)*cos(x[7])^2 + L*x[13]*x[18]*cos(x[7])*sin(x[4])*sin(x[5])*sin(x[7])*sin(x[8])*m__1/3 + ((x[17] - g__1)*(L*x[12]*x[13]*sin(x[4])*sin(x[5])*m__1 - ((m__1 + m__p)*(x[17] - g__1))/9))/2)*cos(x[6])^2 + ((4*L*x[12]*x[18]*cos(x[8])^2*sin(x[5])*sin(x[7])*m__1)/3 + (-L*x[13]*x[18]*cos(x[7])*sin(x[4])*sin(x[5])*m__1/3 - (x[17] - g__1)*sin(x[7])*(L*x[12]*x[13]*sin(x[4])*sin(x[5])*m__1 - ((m__1 + m__p)*(x[17] - g__1))/9)*sin(x[8]))*cos(x[8]) + m__1*sin(x[5])*L*(sin(x[8])*((x[17] - g__1)*x[12]^2 + (x[17]/2 - g__1/2)*x[13]^2 - x[19]/6)*cos(x[7]) - (2*x[12]*x[18]*sin(x[7]))/3))*sin(x[6])*cos(x[6]) + cos(x[8])*((x[17] - g__1)*(L*x[12]*x[13]*sin(x[4])*sin(x[5])*m__1 - ((m__1 + m__p)*(x[17] - g__1))/9)*cos(x[8]) + (4*L*x[12]*x[18]*sin(x[5])*sin(x[8])*m__1)/3)/2)*x[13]*cos(x[4])^4 + ((-6*m__1*sin(x[5])*(cos(x[7])^2 - 2)*L*(x[12]^2*x[18] + 5/2*x[13]^2*x[18] - 1/6*x[20])*cos(x[8])^2 + (36*(m__1*sin(x[4])*x[13]*L*((x[17] - g__1)*x[12]^2 + (x[17]/2 - g__1/2)*x[13]^2 - x[19]/6)*sin(x[5]) - (11*x[12]*(x[17] - g__1)^2*(m__1 + m__p))/36)*sin(x[8])*cos(x[7])^2 - 24*sin(x[7])*x[18]*(L*x[12]*x[13]*sin(x[4])*sin(x[5])*m__1 - ((m__1 + m__p)*(x[17] - g__1))/6)*cos(x[7]) - 72*(m__1*sin(x[4])*x[13]*L*((x[17] - g__1)*x[12]^2 + (x[17]/2 - g__1/2)*x[13]^2 - x[19]/6)*sin(x[5]) - (11*x[12]*(x[17] - g__1)^2*(m__1 + m__p))/36)*sin(x[8]))*cos(x[8]) + 3*L*x[13]^2*x[18]*cos(x[7])^2*sin(x[5])*m__1 + 27*(x[17] - g__1)*x[13]*sin(x[7])*(L*x[12]*x[13]*sin(x[5])*m__1 - (4*sin(x[4])*(m__1 + m__p)*(x[17] - g__1))/27)*sin(x[8])*cos(x[7]) - 6*m__1*sin(x[5])*L*(x[12]^2*x[18] + 5/2*x[13]^2*x[18] - 1/6*x[20]))*cos(x[6])^2 + (72*(m__1*sin(x[4])*x[13]*L*((x[17] - g__1)*x[12]^2 + (x[17]/2 - g__1/2)*x[13]^2 - x[19]/6)*sin(x[5]) - (11*x[12]*(x[17] - g__1)^2*(m__1 + m__p))/36)*sin(x[7])*sin(x[6])*cos(x[8])^2 + (-27*(x[17] - g__1)*x[13]*(L*x[12]*x[13]*sin(x[5])*m__1 - (4*sin(x[4])*(m__1 + m__p)*(x[17] - g__1))/27)*sin(x[6])*cos(x[7]) + 12*m__1*sin(x[7])*L*(x[12]^2*x[18] + 5/2*x[13]^2*x[18] - 1/6*x[20])*sin(x[8])*sin(x[6])*sin(x[5]) + x[15]*(x[17] - g__1)^2*(m__1 + m__p))*cos(x[8]) - 36*((2*x[18]*sin(x[8])*(L*x[12]*x[13]*sin(x[4])*sin(x[5])*m__1 - ((m__1 + m__p)*(x[17] - g__1))/6)*cos(x[7]))/3 + (m__1*sin(x[4])*x[13]*L*((x[17] - g__1)*x[12]^2 + (x[17]/2 - g__1/2)*x[13]^2 - x[19]/6)*sin(x[5]) - (11*x[12]*(x[17] - g__1)^2*(m__1 + m__p))/36)*sin(x[7]))*sin(x[6]))*cos(x[6]) - 6*m__1*sin(x[5])*L*(x[12]^2*x[18] + 5/2*x[13]^2*x[18] - 1/6*x[20])*cos(x[8])^2 + 36*(m__1*sin(x[4])*x[13]*L*((x[17] - g__1)*x[12]^2 + (x[17]/2 - g__1/2)*x[13]^2 - x[19]/6)*sin(x[5]) - (11*x[12]*(x[17] - g__1)^2*(m__1 + m__p))/36)*sin(x[8])*cos(x[8]) + x[14]*sin(x[8])*(x[17] - g__1)^2*(m__1 + m__p)*cos(x[7]) + 6*m__1*(x[12]^2*x[18] + 3/2*x[13]^2*x[18] - 1/6*x[20])*L*sin(x[5]) + x[15]*sin(x[6])*sin(x[7])*sin(x[8])*(x[17] - g__1)^2*(m__1 + m__p))*cos(x[4])^3 + ((-6*(cos(x[7])^2 - 2)*(m__1*sin(x[4])*x[12]*((x[17] - g__1)*x[12]^2 + ((15*x[17])/2 - (15*g__1)/2)*x[13]^2 - x[19]/2)*L*sin(x[5]) - (11*x[13]*(x[17] - g__1)^2*(m__1 + m__p))/6)*cos(x[8])^2 + (-24*(L*x[12]*x[13]*sin(x[5])*m__1 - sin(x[4])*(m__1 + m__p)*(x[17] - g__1)/6)*x[18]*sin(x[8])*cos(x[7])^2 - 42*(x[17] - g__1)*(L*x[13]*m__1*(x[12]^2 - x[13]^2/7)*sin(x[5]) - (11*x[12]*sin(x[4])*(m__1 + m__p)*(x[17] - g__1))/42)*sin(x[7])*cos(x[7]) + 48*(L*x[12]*x[13]*sin(x[5])*m__1 - sin(x[4])*(m__1 + m__p)*(x[17] - g__1)/6)*x[18]*sin(x[8]))*cos(x[8]) + 9*(x[17] - g__1)*(L*x[12]*x[13]*sin(x[4])*sin(x[5])*m__1 + 5*(m__1 + m__p)*(x[17] - g__1)/9)*x[13]*cos(x[7])^2 + 3*L*x[13]^2*x[18]*cos(x[7])*sin(x[4])*sin(x[5])*sin(x[7])*sin(x[8])*m__1 - 6*m__1*sin(x[4])*x[12]*((x[17] - g__1)*x[12]^2 + ((15*x[17])/2 - (15*g__1)/2)*x[13]^2 - x[19]/2)*L*sin(x[5]) + 11*x[13]*(x[17] - g__1)^2*(m__1 + m__p))*cos(x[6])^2 + (-48*(L*x[12]*x[13]*sin(x[5])*m__1 - sin(x[4])*(m__1 + m__p)*(x[17] - g__1)/6)*sin(x[7])*x[18]*sin(x[6])*cos(x[8])^2 + 12*(-L*x[13]^2*x[18]*cos(x[7])*sin(x[4])*sin(x[5])*m__1/4 + sin(x[7])*(m__1*sin(x[4])*x[12]*((x[17] - g__1)*x[12]^2 + ((15*x[17])/2 - (15*g__1)/2)*x[13]^2 - x[19]/2)*L*sin(x[5]) - (11*x[13]*(x[17] - g__1)^2*(m__1 + m__p))/6)*sin(x[8]))*sin(x[6])*cos(x[8]) - 42*(x[17] - g__1)*(L*x[13]*m__1*(x[12]^2 - x[13]^2/7)*sin(x[5]) - (11*x[12]*sin(x[4])*(m__1 + m__p)*(x[17] - g__1))/42)*sin(x[8])*sin(x[6])*cos(x[7]) + 24*m__1*(-x[13]*x[15]*sin(x[8])*(x[17] - g__1)*sin(x[4])/24 + x[12]*x[18]*sin(x[6])*sin(x[7]))*x[13]*L*sin(x[5]) - 4*x[18]*sin(x[4])*sin(x[6])*sin(x[7])*(m__1 + m__p)*(x[17] - g__1))*cos(x[6]) + (-6*m__1*sin(x[4])*x[12]*((x[17] - g__1)*x[12]^2 + ((15*x[17])/2 - (15*g__1)/2)*x[13]^2 - x[19]/2)*L*sin(x[5]) + 11*x[13]*(x[17] - g__1)^2*(m__1 + m__p))*cos(x[8])^2 + (L*x[13]^2*x[14]*sin(x[4])*sin(x[5])*m__1*(x[17] - g__1)*cos(x[7]) - 24*m__1*(-x[13]*x[15]*sin(x[6])*sin(x[7])*(x[17] - g__1)*sin(x[4])/24 + x[12]*x[18]*sin(x[8]))*x[13]*L*sin(x[5]) + 4*x[18]*sin(x[4])*sin(x[8])*(m__1 + m__p)*(x[17] - g__1))*cos(x[8]) + x[15]*sin(x[4])*sin(x[6])*(x[17] - g__1)^2*(m__1 + m__p)*cos(x[7]) + 6*m__1*sin(x[4])*((x[17] - g__1)*x[12]^2 + ((9*x[17])/2 - (9*g__1)/2)*x[13]^2 - x[19]/2)*x[12]*L*sin(x[5]) - 9*(x[17] - g__1)^2*(x[14]*sin(x[4])*sin(x[7])/9 + x[13])*(m__1 + m__p))*cos(x[4])^2 + ((6*L*sin(x[5])*m__1*x[12]^2*x[18]*(cos(x[7])^2 - 2)*cos(x[8])^2 - 36*(x[17] - g__1)*((L*x[12]*x[13]*sin(x[4])*sin(x[5])*m__1 - 2*(m__1 + m__p)*(x[17] - g__1)/9)*x[12]*sin(x[8])*cos(x[7])^2 - x[18]*sin(x[7])*(m__1 + m__p)*cos(x[7])/9 - 2*(L*x[12]*x[13]*sin(x[4])*sin(x[5])*m__1 - 2*(m__1 + m__p)*(x[17] - g__1)/9)*x[12]*sin(x[8]))*cos(x[8]) + 6*(x[17] - g__1)*x[13]*sin(x[7])*(L*x[12]*x[13]*sin(x[5])*m__1 - (7*sin(x[4])*(m__1 + m__p)*(x[17] - g__1))/6)*sin(x[8])*cos(x[7]) + 6*L*sin(x[5])*m__1*x[12]^2*x[18])*cos(x[6])^2 + (-72*(x[17] - g__1)*(L*x[12]*x[13]*sin(x[4])*sin(x[5])*m__1 - 2*(m__1 + m__p)*(x[17] - g__1)/9)*sin(x[7])*x[12]*sin(x[6])*cos(x[8])^2 + (-6*(x[17] - g__1)*x[13]*(L*x[12]*x[13]*sin(x[5])*m__1 - (7*sin(x[4])*(m__1 + m__p)*(x[17] - g__1))/6)*sin(x[6])*cos(x[7]) - 12*L*x[12]^2*x[18]*sin(x[5])*sin(x[6])*sin(x[7])*sin(x[8])*m__1 - x[15]*(x[17] - g__1)^2*(m__1 + m__p))*cos(x[8]) + 36*(x[17] - g__1)*(x[18]*sin(x[8])*(m__1 + m__p)*cos(x[7])/9 + (L*x[12]*x[13]*sin(x[4])*sin(x[5])*m__1 - 2*(m__1 + m__p)*(x[17] - g__1)/9)*sin(x[7])*x[12])*sin(x[6]))*cos(x[6]) + 6*L*sin(x[5])*cos(x[8])^2*m__1*x[12]^2*x[18] - 36*(x[17] - g__1)*(L*x[12]*x[13]*sin(x[4])*sin(x[5])*m__1 - 2*(m__1 + m__p)*(x[17] - g__1)/9)*x[12]*sin(x[8])*cos(x[8]) - x[14]*sin(x[8])*(x[17] - g__1)^2*(m__1 + m__p)*cos(x[7]) - 6*L*sin(x[5])*m__1*x[12]^2*x[18] - x[15]*sin(x[6])*sin(x[7])*sin(x[8])*(x[17] - g__1)^2*(m__1 + m__p))*cos(x[4]) + 6*(x[17] - g__1)*sin(x[4])*x[12]*((L*x[12]^2*sin(x[5])*m__1*(cos(x[7])^2 - 2)*cos(x[8])^2 + (7*cos(x[7])*sin(x[7])*(m__1 + m__p)*(x[17] - g__1)*cos(x[8]))/6 + L*x[12]^2*sin(x[5])*m__1)*cos(x[6])^2 - 2*(L*x[12]^2*cos(x[8])*sin(x[5])*sin(x[7])*m__1 - (7*cos(x[7])*(m__1 + m__p)*(x[17] - g__1))/12)*sin(x[8])*sin(x[6])*cos(x[6]) + L*x[12]^2*sin(x[5])*m__1*(cos(x[8]) - 1)*(cos(x[8]) + 1)))*cos(x[5])^3 + (6*(((cos(x[7])^2 - 2)*cos(x[8])^2 - 2*cos(x[7])^2 + 1)*cos(x[6])^2 - 2*cos(x[6])*sin(x[7])*cos(x[8])*sin(x[6])*sin(x[8]) + cos(x[8])^2)*(x[17] - g__1)*m__1*x[13]^3*L*cos(x[4])^5 + 36*m__1*x[13]^2*(((x[12]*sin(x[8])*(x[17] - g__1)*cos(x[7])^2 - x[18]*cos(x[7])*sin(x[7])/3 - 2*x[12]*sin(x[8])*(x[17] - g__1))*cos(x[8]) + x[13]*cos(x[7])*sin(x[4])*sin(x[7])*sin(x[8])*(x[17] - g__1)/3)*cos(x[6])^2 + 2*(x[12]*sin(x[7])*(x[17] - g__1)*cos(x[8])^2 - x[13]*sin(x[4])*cos(x[7])*(x[17] - g__1)*cos(x[8])/6 - x[18]*cos(x[7])*sin(x[8])/6 - x[12]*sin(x[7])*(x[17] - g__1)/2)*sin(x[6])*cos(x[6]) + x[12]*cos(x[8])*sin(x[8])*(x[17] - g__1))*L*cos(x[4])^4 - 18*x[13]*((m__1*((x[17] - g__1)*x[12]^2 + ((7*x[17])/6 - (7*g__1)/6)*x[13]^2 - x[19]/6)*(cos(x[7])^2 - 2)*L*cos(x[8])^2 + ((2*L*x[13]*x[18]*cos(x[7])^2*sin(x[4])*sin(x[8])*m__1)/3 + 2*(x[17] - g__1)*(-11*(m__1 + m__p)*(x[17] - g__1)*sin(x[5])/36 + L*x[12]*x[13]*sin(x[4])*m__1)*sin(x[7])*cos(x[7]) - (4*L*x[13]*x[18]*sin(x[4])*sin(x[8])*m__1)/3)*cos(x[8]) + m__1*(((x[17] - g__1)*x[13]^2*cos(x[7])^2)/6 + (x[17] - g__1)*x[12]^2 + ((7*x[17])/6 - (7*g__1)/6)*x[13]^2 - x[19]/6)*L)*cos(x[6])^2 - 2*(-(2*L*x[13]*x[18]*cos(x[8])^2*sin(x[4])*sin(x[7])*m__1)/3 + m__1*((x[17] - g__1)*x[12]^2 + ((7*x[17])/6 - (7*g__1)/6)*x[13]^2 - x[19]/6)*sin(x[7])*L*sin(x[8])*cos(x[8]) - (x[17] - g__1)*(-11*(m__1 + m__p)*(x[17] - g__1)*sin(x[5])/36 + L*x[12]*x[13]*sin(x[4])*m__1)*sin(x[8])*cos(x[7]) + L*x[13]*x[18]*sin(x[4])*sin(x[7])*m__1/3)*sin(x[6])*cos(x[6]) + (((x[17] - g__1)*x[12]^2 + ((7*x[17])/6 - (7*g__1)/6)*x[13]^2 - x[19]/6)*cos(x[8])^2 + (2*x[13]*x[18]*cos(x[8])*sin(x[4])*sin(x[8]))/3 + (g__1 - x[17])*x[12]^2 + (-(5*x[17])/6 + (5*g__1)/6)*x[13]^2 + x[19]/6)*m__1*L)*cos(x[4])^3 + ((12*(-((m__1 + m__p)*(x[17] - g__1)*sin(x[5]))/3 + L*x[12]*x[13]*sin(x[4])*m__1)*(cos(x[7])^2 - 2)*x[18]*cos(x[8])^2 - 36*x[13]*((x[17] - g__1)*sin(x[8])*(-(11*sin(x[4])*(m__1 + m__p)*(x[17] - g__1)*sin(x[5]))/36 + L*x[12]*x[13]*m__1)*cos(x[7])^2 + L*x[13]*x[18]*cos(x[7])*sin(x[7])*m__1/12 - 2*(x[17] - g__1)*sin(x[8])*(-(11*sin(x[4])*(m__1 + m__p)*(x[17] - g__1)*sin(x[5]))/36 + L*x[12]*x[13]*m__1))*cos(x[8]) - 4*x[18]*sin(x[5])*(m__1 + m__p)*(x[17] - g__1)*cos(x[7])^2 + 9*(x[17] - g__1)*(-(10*x[12]*(m__1 + m__p)*(x[17] - g__1)*sin(x[5]))/9 + L*x[13]^3*sin(x[4])*m__1)*sin(x[7])*sin(x[8])*cos(x[7]) + 12*(-((m__1 + m__p)*(x[17] - g__1)*sin(x[5]))/3 + L*x[12]*x[13]*sin(x[4])*m__1)*x[18])*cos(x[6])^2 + (-72*(x[17] - g__1)*x[13]*sin(x[7])*(-(11*sin(x[4])*(m__1 + m__p)*(x[17] - g__1)*sin(x[5]))/36 + L*x[12]*x[13]*m__1)*sin(x[6])*cos(x[8])^2 + (-9*(x[17] - g__1)*(-(10*x[12]*(m__1 + m__p)*(x[17] - g__1)*sin(x[5]))/9 + L*x[13]^3*sin(x[4])*m__1)*sin(x[6])*cos(x[7]) + 8*x[18]*sin(x[6])*sin(x[7])*sin(x[8])*(m__1 + m__p)*(x[17] - g__1)*sin(x[5]) - 24*m__1*(x[12]*x[18]*sin(x[4])*sin(x[6])*sin(x[7])*sin(x[8]) - x[13]*x[15]*(x[17] - g__1)/24)*x[13]*L)*cos(x[8]) + 36*(-L*x[13]*x[18]*cos(x[7])*sin(x[8])*m__1/12 + (x[17] - g__1)*sin(x[7])*(-(11*sin(x[4])*(m__1 + m__p)*(x[17] - g__1)*sin(x[5]))/36 + L*x[12]*x[13]*m__1))*x[13]*sin(x[6]))*cos(x[6]) + 12*(-((m__1 + m__p)*(x[17] - g__1)*sin(x[5]))/3 + L*x[12]*x[13]*sin(x[4])*m__1)*x[18]*cos(x[8])^2 - 36*(x[17] - g__1)*x[13]*sin(x[8])*(-(11*sin(x[4])*(m__1 + m__p)*(x[17] - g__1)*sin(x[5]))/36 + L*x[12]*x[13]*m__1)*cos(x[8]) + L*x[13]^2*x[14]*sin(x[8])*m__1*(x[17] - g__1)*cos(x[7]) + 4*x[18]*sin(x[5])*(m__1 + m__p)*(x[17] - g__1) - 12*m__1*x[13]*L*(x[12]*x[18]*sin(x[4]) - x[13]*x[15]*sin(x[6])*sin(x[7])*sin(x[8])*(x[17] - g__1)/12))*cos(x[4])^2 + 18*(x[17] - g__1)*(((cos(x[7])^2 - 2)*x[12]*(-sin(x[4])*(m__1 + m__p)*(x[17] - g__1)*sin(x[5])/3 + L*x[12]*x[13]*m__1)*cos(x[8])^2 - ((-7*(m__1 + m__p)*(x[17] - g__1)*sin(x[5])/6 + L*x[12]*x[13]*sin(x[4])*m__1)*x[13]*sin(x[7])*cos(x[7])*cos(x[8]))/3 - (5*x[12]*sin(x[4])*sin(x[5])*(m__1 + m__p)*(x[17] - g__1)*cos(x[7])^2)/9 + (2*x[18]*sin(x[4])*sin(x[5])*sin(x[7])*sin(x[8])*(m__1 + m__p)*cos(x[7]))/9 + x[12]*(-sin(x[4])*(m__1 + m__p)*(x[17] - g__1)*sin(x[5])/3 + L*x[12]*x[13]*m__1))*cos(x[6])^2 + (-2*(x[18]*sin(x[4])*sin(x[5])*(m__1 + m__p)*cos(x[7])/9 + sin(x[7])*x[12]*(-sin(x[4])*(m__1 + m__p)*(x[17] - g__1)*sin(x[5])/3 + L*x[12]*x[13]*m__1)*sin(x[8]))*sin(x[6])*cos(x[8]) - (((-7*(m__1 + m__p)*(x[17] - g__1)*sin(x[5])/6 + L*x[12]*x[13]*sin(x[4])*m__1)*x[13]*sin(x[6])*cos(x[7]) + x[15]*sin(x[4])*sin(x[5])*(m__1 + m__p)*(x[17] - g__1)/6)*sin(x[8]))/3)*cos(x[6]) + x[12]*(-sin(x[4])*(m__1 + m__p)*(x[17] - g__1)*sin(x[5])/3 + L*x[12]*x[13]*m__1)*cos(x[8])^2 + sin(x[4])*sin(x[5])*(m__1 + m__p)*(x[17] - g__1)*(sin(x[6])*sin(x[7])*x[15] + x[14]*cos(x[7]))*cos(x[8])/18 - x[12]*(-sin(x[4])*(m__1 + m__p)*(x[17] - g__1)*sin(x[5])/3 + L*x[12]*x[13]*m__1))*cos(x[4]) + 7*x[12]*cos(x[6])*cos(x[7])*sin(x[5])*(x[17] - g__1)^2*(m__1 + m__p)*(cos(x[6])*sin(x[7])*sin(x[8]) - sin(x[6])*cos(x[8])))*cos(x[5])^2 + (-12*L*x[13]^3*cos(x[6])*cos(x[7])*sin(x[5])*m__1*(cos(x[6])*sin(x[7])*cos(x[8]) + sin(x[6])*sin(x[8]))*(x[17] - g__1)*cos(x[4])^4 - 12*m__1*sin(x[5])*x[13]^2*(((-x[18]*cos(x[7])^2/2 + x[18])*cos(x[8])^2 + x[13]*sin(x[4])*sin(x[8])*(cos(x[7])^2 - 2)*(x[17] - g__1)*cos(x[8]) - x[18]/2)*cos(x[6])^2 + 2*(x[13]*sin(x[4])*(x[17] - g__1)*cos(x[8])^2 + x[18]*cos(x[8])*sin(x[8])/2 - x[13]*sin(x[4])*(x[17] - g__1)/2)*sin(x[7])*sin(x[6])*cos(x[6]) - x[18]*cos(x[8])^2/2 + x[13]*sin(x[4])*sin(x[8])*(x[17] - g__1)*cos(x[8]) + x[18]/2)*L*cos(x[4])^3 + 18*(x[17] - g__1)*x[13]*(((L*x[12]*x[13]*sin(x[4])*sin(x[5])*m__1 - ((m__1 + m__p)*(x[17] - g__1))/3)*(cos(x[7])^2 - 2)*cos(x[8])^2 - L*x[13]^2*cos(x[7])*cos(x[8])*sin(x[5])*sin(x[7])*m__1/3 - 5*(m__1 + m__p)*(x[17] - g__1)*cos(x[7])^2/9 + L*x[12]*x[13]*sin(x[4])*sin(x[5])*m__1 - ((m__1 + m__p)*(x[17] - g__1))/3)*cos(x[6])^2 - 2*sin(x[8])*((L*x[12]*x[13]*sin(x[4])*sin(x[5])*m__1 - ((m__1 + m__p)*(x[17] - g__1))/3)*sin(x[7])*cos(x[8]) + L*x[13]^2*cos(x[7])*sin(x[5])*m__1/6)*sin(x[6])*cos(x[6]) + (cos(x[8]) + 1)*(L*x[12]*x[13]*sin(x[4])*sin(x[5])*m__1 - ((m__1 + m__p)*(x[17] - g__1))/3)*(cos(x[8]) - 1))*cos(x[4])^2 + 10*(x[17] - g__1)*(sin(x[7])*cos(x[7])*(-(2*x[18]*cos(x[8]))/5 + x[13]*sin(x[4])*sin(x[8])*(x[17] - g__1))*cos(x[6])^2 + (-(x[17] - g__1)*(x[13]*cos(x[7])*sin(x[4])*sin(x[6]) - x[15]/10)*cos(x[8]) - (2*x[18]*cos(x[7])*sin(x[6])*sin(x[8]))/5)*cos(x[6]) + sin(x[8])*(x[17] - g__1)*(sin(x[6])*sin(x[7])*x[15] + x[14]*cos(x[7]))/10)*(m__1 + m__p)*cos(x[4]) - 7*x[12]*cos(x[6])*cos(x[7])*sin(x[4])*(x[17] - g__1)^2*(m__1 + m__p)*(cos(x[6])*sin(x[7])*cos(x[8]) + sin(x[6])*sin(x[8])))*cos(x[5]) + 6*(x[17] - g__1)*cos(x[4])*x[13]*(((1 + (cos(x[7])^2 - 2)*cos(x[8])^2)*cos(x[6])^2 - 2*cos(x[6])*sin(x[7])*cos(x[8])*sin(x[6])*sin(x[8]) + cos(x[8])^2 - 1)*m__1*x[13]^2*L*cos(x[4])^2 - (7*cos(x[6])*cos(x[7])*sin(x[5])*(m__1 + m__p)*(cos(x[6])*sin(x[7])*cos(x[8]) + sin(x[6])*sin(x[8]))*(x[17] - g__1))/6))/(cos(x[4])^4*m__1*(((((cos(x[7])^2 - 2)*cos(x[8])^2 - 2*cos(x[7])^2 + 1)*cos(x[6])^2 - 2*cos(x[6])*sin(x[7])*cos(x[8])*sin(x[6])*sin(x[8]) + cos(x[8])^2)*cos(x[4])^2 + 2*cos(x[6])*cos(x[7])*sin(x[4])*(cos(x[6])*sin(x[7])*sin(x[8]) - sin(x[6])*cos(x[8]))*cos(x[4]) - 2*(cos(x[8])^2 - 1/2)*(cos(x[7])^2 - 2)*cos(x[6])^2 + 4*cos(x[6])*sin(x[7])*cos(x[8])*sin(x[6])*sin(x[8]) - 2*cos(x[8])^2 + 1)*cos(x[5])^2 - 2*sin(x[5])*(cos(x[6])*cos(x[7])*(cos(x[6])*sin(x[7])*cos(x[8]) + sin(x[6])*sin(x[8]))*cos(x[4]) + sin(x[4])*(cos(x[8])*sin(x[8])*(cos(x[7])^2 - 2)*cos(x[6])^2 + 2*(cos(x[8])^2 - 1/2)*sin(x[7])*sin(x[6])*cos(x[6]) + cos(x[8])*sin(x[8])))*cos(x[5]) + (1 + (cos(x[7])^2 - 2)*cos(x[8])^2)*cos(x[6])^2 - 2*cos(x[6])*sin(x[7])*cos(x[8])*sin(x[6])*sin(x[8]) + cos(x[8])^2 - 1)*L*cos(x[5])^4)