MATLAB Answers

How can I symbolically solve this system of equations?

2 views (last 30 days)
onamaewa
onamaewa on 24 May 2019
Answered: Star Strider on 24 May 2019
syms B1 B2 g m1 m2 K1 K2 P v1 v2 x1 x2
x1(dot) = v1
v1(dot) = (1/m1) * (-K1 * x1 - B1 * v1 + K1 * x2 + P)
x2(dot) = v2
v2(dot) = (1/m2) * (K1 * x1 - (K1 + K2) * x2 - B2 * v2 - m2 * g)

  0 Comments

Sign in to comment.

Answers (1)

Star Strider
Star Strider on 24 May 2019
Try this:
syms B1 B2 g m1 m2 K1 K2 P v1(t) v2(t) x1(t) x2(t) x10 x20 v10 v20
Eq1 = diff(x1) == v1;
Eq2 = diff(v1) == (1/m1) * (-K1 * x1 - B1 * v1 + K1 * x2 + P);
Eq3 = diff(x2) == v2;
Eq4 = diff(v2) == (1/m2) * (K1 * x1 - (K1 + K2) * x2 - B2 * v2 - m2 * g);
S = dsolve(Eq1,Eq2,Eq3,Eq4, x1(0)==x10, x2(0)==x20, v1(0)==v10, v2(0)==v20);
it then produces:
V1 =
exp(t*root(#X^4*m1*m2 + #X^3*B2*m1 + #X^3*B1*m2 + #X^2*K2*m1 + #X^2*K1*m2 + #X^2*K1*m1 + #X^2*B1*B2 + #X*B2*K1 + #X*B1*K2 + #X*B1*K1 + K1*K2, #X, 4))*((exp(-t*root(#X^4*m1*m2 + #X^3*B2*m1 + #X^3*B1*m2 + #X^2*K2*m1 + #X^2*K1*m2 + #X^2*K1*m1 + #X^2*B1*B2 + #X*B2*K1 + #X*B1*K2 + #X*B1*K1 + K1*K2, #X, 4))*(B1*K1*P*m2^2 - B2^3*g*m1^2 + B2*K1*P*m1*m2 + K1*P*m1*m2^2*root(#X^4*m1*m2 + #X^3*B2*m1 + #X^3*B1*m2 + #X^2*K2*m1 + #X^2*K1*m2 + #X^2*K1*m1 + #X^2*B1*B2 + #X*B2*K1 + #X*B1*K2 + #X*B1*K1 + K1*K2, #X, 1) + K1*P*m1*m2^2*root(#X^4*m1*m2 + #X^3*B2*m1 + #X^3*B1*m2 + #X^2*K2*m1 + #X^2*K1*m2 + #X^2*K1*m1 + #X^2*B1*B2 + #X*B2*K1 + #X*B1*K2 + #X*B1*K1 + K1*K2, #X, 2) + K1*P*m1*m2^2*root(#X^4*m1*m2 + #X^3*B2*m1 + #X^3*B1*m2 + #X^2*K2*m1 + #X^2*K1*m2 + #X^2*K1*m1 + #X^2*B1*B2 + #X*B2*K1 + #X*B1*K2 + #X*B1*K1 + K1*K2, #X, 3) + 2*B2*K1*g*m1^2*m2 + 2*B2*K2*g*m1^2*m2 - g*m2*B2^2*m1^2*root(#X^4*m1*m2 + #X^3*B2*m1 + #X^3*B1*m2 + #X^2*K2*m1 + #X^2*K1*m2 + #X^2*K1*m1 + #X^2*B1*B2 + #X*B2*K1 + #X*B1*K2 + #X*B1*K1 + K1*K2, #X, 1) - g*m2*B2^2*m1^2*root(#X^4*m1*m2 + #X^3*B2*m1 + #X^3*B1*m2 + #X^2*K2*m1 + #X^2*K1*m2 + #X^2*K1*m1 + #X^2*B1*B2 + #X*B2*K1 + #X*B1*K2 + #X*B1*K1 + K1*K2, #X, 2) - g*m2*B2^2*m1^2*root(#X^4*m1*m2 + #X^3*B2*m1 + #X^3*B1*m2 + #X^2*K2*m1 + #X^2*K1*m2 + #X^2*K1*m1 + #X^2*B1*B2 + #X*B2*K1 + #X*B1*K2 + #X*B1*K1 + K1*K2, #X, 3) + K1*g*m1^2*m2^2*root(#X^4*m1*m2 + #X^3*B2*m1 + #X^3*B1*m2 + #X^2*K2*m1 + #X^2*K1*m2 + #X^2*K1*m1 + #X^2*B1*B2 + #X*B2*K1 + #X*B1*K2 + #X*B1*K1 + K1*K2, #X, 1) + K1*g*m1^2*m2^2*root(#X^4*m1*m2 + #X^3*B2*m1 + #X^3*B1*m2 + #X^2*K2*m1 + #X^2*K1*m2 + #X^2*K1*m1 + #X^2*B1*B2 + #X*B2*K1 + #X*B1*K2 + #X*B1*K1 + K1*K2, #X, 2) + K2*g*m1^2*m2^2*root(#X^4*m1*m2 + #X^3*B2*m1 + #X^3*B1*m2 + #X^2*K2*m1 + #X^2*K1*m2 + #X^2*K1*m1 + #X^2*B1*B2 + #X*B2*K1 + #X*B1*K2 + #X*B1*K1 + K1*K2, #X, 1) + K1*g*m1^2*m2^2*root(#X^4*m1*m2 + #X^3*B2*m1 + #X^3*B1*m2 + #X^2*K2*m1 + #X^2*K1*m2 + #X^2*K1*m1 + #X^2*B1*B2 + #X*B2*K1 + #X*B1*K2 + #X*B1*K1 + K1*K2, #X, 3) + K2*g*m1^2*m2^2*root(#X^4*m1*m2 + #X^3*B2*m1 + #X^3*B1*m2 + #X^2*K2*m1 + #X^2*K1*m2 + #X^2*K1*m1 + #X^2*B1*B2 + #X*B2*K1 + #X*B1*K2 + #X*B1*K1 + K1*K2, #X, 2) + K2*g*m1^2*m2^2*root(#X^4*m1*m2 + #X^3*B2*m1 + #X^3*B1*m2 + #X^2*K2*m1 + #X^2*K1*m2 + #X^2*K1*m1 + #X^2*B1*B2 + #X*B2*K1 + #X*B1*K2 + #X*B1*K1 + K1*K2, #X, 3) - g*m1^2*m2^3*root(#X^4*m1*m2 + #X^3*B2*m1 + #X^3*B1*m2 + #X^2*K2*m1 + #X^2*K1*m2 + #X^2*K1*m1 + #X^2*B1*B2 + #X*B2*K1 + #X*B1*K2 + #X*B1*K1 + K1*K2, #X, 1)*root(#X^4*m1*m2 + #X^3*B2*m1 + #X^3*B1*m2 + #X^2*K2*m1 + #X^2*K1*m2 + #X^2*K1*m1 + #X^2*B1*B2 + #X*B2*K1 + #X*B1*K2 + #X*B1*K1 + K1*K2, #X, 2)*root(#X^4*m1*m2 + #X^3*B2*m1 + #X^3*B1*m2 + #X^2*K2*m1 + #X^2*K1*m2 + #X^2*K1*m1 + #X^2*B1*B2 + #X*B2*K1 + #X*B1*K2 + #X*B1*K1 + K1*K2, #X, 3) - B2*g*m1^2*m2^2*root(#X^4*m1*m2 + #X^3*B2*m1 + #X^3*B1*m2 + #X^2*K2*m1 + #X^2*K1*m2 + #X^2*K1*m1 + #X^2*B1*B2 + #X*B2*K1 + #X*B1*K2 + #X*B1*K1 + K1*K2, #X, 1)*root(#X^4*m1*m2 + #X^3*B2*m1 + #X^3*B1*m2 + #X^2*K2*m1 + #X^2*K1*m2 + #X^2*K1*m1 + #X^2*B1*B2 + #X*B2*K1 + #X*B1*K2 + #X*B1*K1 + K1*K2, #X, 2) - B2*g*m1^2*m2^2*root(#X^4*m1*m2 + #X^3*B2*m1 + #X^3*B1*m2 + #X^2*K2*m1 + #X^2*K1*m2 + #X^2*K1*m1 + #X^2*B1*B2 + #X*B2*K1 + #X*B1*K2 + #X*B1*K1 + K1*K2, #X, 1)*root(#X^4*m1*m2 + #X^3*B2*m1 + #X^3*B1*m2 + #X^2*K2*m1 + #X^2*K1*m2 + #X^2*K1*m1 + #X^2*B1*B2 + #X*B2*K1 + #X*B1*K2 + #X*B1*K1 + K1*K2, #X, 3) - B2*g*m1^2*m2^2*root(#X^4*m1*m2 + #X^3*B2*m1 + #X^3*B1*m2 + #X^2*K2*m1 + #X^2*K1*m2 + #X^2*K1*m1 + #X^2*B1*B2 + #X*B2*K1 + #X*B1*K2 + #X*B1*K1 + K1*K2, #X, 2)*root(#X^4*m1*m2 + #X^3*B2*m1 + #X^3*B1*m2 + #X^2*K2*m1 + #X^2*K1*m2 + #X^2*K1*m1 + #X^2*B1*B2 + #X*B2*K1 + #X*B1*K2 + #X*B1*K1 + K1*K2, #X, 3)))/(root(#X^4*m1*m2 + #X^3*B2*m1 + #X^3*B1*m2 + #X^2*K2*m1 + #X^2*K1*m2 + #X^2*K1*m1 + #X^2*B1*B2 + #X*B2*K1 + #X*B1*K2 + #X*B1*K1 + K1*K2, #X, 4)*(m1^2*m2^3*root(#X^4*m1*m2 + #X^3*B2*m1 + #X^3*B1*m2 + #X^2*K2*m1 + #X^2*K1*m2 + #X^2*K1*m1 + #X^2*B1*B2 + #X*B2*K1 + #X*B1*K2 + #X*B1*K1 + K1*K2, #X, 4)^3 - root(#X^4*m1*m2 + #X^3*B2*m1 + #X^3*B1*m2 + #X^2*K2*m1 + #X^2*K1*m2 + #X^2*K1*m1 + #X^2*B1*B2 + #X*B2*K1 + #X*B1*K2 + #X*B1*K1 + K1*K2, #X, 4)^2*m1^2*m2^3*root(#X^4*m1*m2 + #X^3*B2*m1 + #X^3*B1*m2 + #X^2*K2*m1 + #X^2*K1*m2 + #X^2*K1*m1 + #X^2*B1*B2 + #X*B2*K1 + #X*B1*K2 + #X*B1*K1 + K1*K2, #X, 1) - root(#X^4*m1*m2 + #X^3*B2*m1 + #X^3*B1*m2 + #X^2*K2*m1 + #X^2*K1*m2 + #X^2*K1*m1 + #X^2*B1*B2 + #X*B2*K1 + #X*B1*K2 + #X*B1*K1 + K1*K2, #X, 4)^2*m1^2*m2^3*root(#X^4*m1*m2 + #X^3*B2*m1 + #X^3*B1*m2 + #X^2*K2*m1 + #X^2*K1*m2 + #X^2*K1*m1 + #X^2*B1*B2 + #X*B2*K1 + #X*B1*K2 + #X*B1*K1 + K1*K2, #X, 2) - root(#X^4*m1*m2 + #X^3*B2*m1 + #X^3*B1*m2 + #X^2*K2*m1 + #X^2*K1*m2 + #X^2*K1*m1 + #X^2*B1*B2 + #X*B2*K1 + #X*B1*K2 + #X*B1*K1 + K1*K2, #X, 4)^2*m1^2*m2^3*root(#X^4*m1*m2 + #X^3*B2*m1 + #X^3*B1*m2 + #X^2*K2*m1 + #X^2*K1*m2 + #X^2*K1*m1 + #X^2*B1*B2 + #X*B2*K1 + #X*B1*K2 + #X*B1*K1 + K1*K2, #X, 3) - m1^2*m2^3*root(#X^4*m1*m2 + #X^3*B2*m1 + #X^3*B1*m2 + #X^2*K2*m1 + #X^2*K1*m2 + #X^2*K1*m1 + #X^2*B1*B2 + #X*B2*K1 + #X*B1*K2 + #X*B1*K1 + K1*K2, #X, 1)*root(#X^4*m1*m2 + #X^3*B2*m1 + #X^3*B1*m2 + #X^2*K2*m1 + #X^2*K1*m2 + #X^2*K1*m1 + #X^2*B1*B2 + #X*B2*K1 + #X*B1*K2 + #X*B1*K1 + K1*K2, #X, 2)*root(#X^4*m1*m2 + #X^3*B2*m1 + #X^3*B1*m2 + #X^2*K2*m1 + #X^2*K1*m2 + #X^2*K1*m1 + #X^2*B1*B2 + #X*B2*K1 + #X*B1*K2 + #X*B1*K1 + K1*K2, #X, 3) + m1^2*m2^3*root(#X^4*m1*m2 + #X^3*B2*m1 + #X^3*B1*m2 + #X^2*K2*m1 + #X^2*K1*m2 + #X^2*K1*m1 + #X^2*B1*B2 + #X*B2*K1 + #X*B1*K2 + #X*B1*K1 + K1*K2, #X, 1)*root(#X^4*m1*m2 + #X^3*B2*m1 + #X^3*B1*m2 + #X^2*K2*m1 + #X^2*K1*m2 + #X^2*K1*m1 + #X^2*B1*B2 + #X*B2*K1 + #X*B1*K2 + #X*B1*K1 + K1*K2, #X, 2)*root(#X^4*m1*m2 + #X^3*B2*m1 + #X^3*B1*m2 + #X^2*K2*m1 + #X^2*K1*m2 + #X^2*K1*m1 + #X^2*B1*B2 + #X*B2*K1 + #X*B1*K2 + #X*B1*K1 + K1*K2, #X, 4) + m1^2*m2^3*root(#X^4*m1*m2 + #X^3*B2*m1 + #X^3*B1*m2 + #X^2*K2*m1 + #X^2*K1*m2 + #X^2*K1*m1 + #X^2*B1*B2 + #X*B2*K1 + #X*B1*K2 + #X*B1*K1 + K1*K2, #X, 1)*root(#X^4*m1*m2 + #X^3*B2*m1 + #X^3*B1*m2 + #X^2*K2*m1 + #X^2*K1*m2 + #X^2*K1*m1 + #X^2*B1*B2 + #X*B2*K1 + #X*B1*K2 + #X*B1*K1 + K1*K2, #X, 3)*root(#X^4*m1*m2 + #X^3*B2*m1 + #X^3*B1*m2 + #X^2*K2*m1 + #X^2*K1*m2 + #X^2*K1*m1 + #X^2*B1*B2 + #X*B2*K1 + #X*B1*K2 + #X*B1*K1 + K1*K2, #X, 4) + m1^2*m2^3*root(#X^4*m1*m2 + #X^3*B2*m1 + #X^3*B1*m2 + #X^2*K2*m1 + #X^2*K1*m2 + #X^2*K1*m1 + #X^2*B1*B2 + #X*B2*K1 + #X*B1*K2 + #X*B1*K1 + K1*K2, #X, 2)*root(#X^4*m1*m2 + #X^3*B2*m1 + #X^3*B1*m2 + #X^2*K2*m1 + #X^2*K1*m2 + #X^2*K1*m1 + #X^2*B1*B2 + #X*B2*K1 + #X*B1*K2 + #X*B1*K1 + K1*K2, #X, 3)*root(#X^4*m1*m2 + #X^3*B2*m1 + #X^3*B1*m2 + #X^2*K2*m1 + #X^2*K1*m2 + #X^2*K1*m1 + #X^2*B1*B2 + #X*B2*K1 + #X*B1*K2 + #X*B1*K1 + K1*K2, #X, 4))) - (B1*K1*P*m2^2 - B2^3*g*m1^2 + v20*B2^3*m1^2*root(#X^4*m1*m2 + #X^3*B2*m1 + #X^3*B1*m2 + #X^2*K2*m1 + #X^2*K1*m2 + #X^2*K1*m1 + #X^2*B1*B2 + #X*B2*K1 + #X*B1*K2 + #X*B1*K1 + K1*K2, #X, 4) + B2*K1*P*m1*m2 + K1*P*m1*m2^2*root(#X^4*m1*m2 + #X^3*B2*m1 + #X^3*B1*m2 + #X^2*K2*m1 + #X^2*K1*m2 + #X^2*K1*m1 + #X^2*B1*B2 + #X*B2*K1 + #X*B1*K2 + #X*B1*K1 + K1*K2, #X, 1) + K1*P*m1*m2^2*root(#X^4*m1*m2 + #X^3*B2*m1 + #X^3*B1*m2 + #X^2*K2*m1 + #X^2*K1*m2 + #X^2*K1*m1 + #X^2*B1*B2 + #X*B2*K1 + #X*B1*K2 + #X*B1*K1 + K1*K2, #X, 2) + K1*P*m1*m2^2*root(#X^4*m1*m2 + #X^3*B2*m1 + #X^3*B1*m2 + #X^2*K2*m1 + #X^2*K1*m2 + #X^2*K1*m1 + #X^2*B1*B2 + #X*B2*K1 + #X*B1*K2 + #X*B1*K1 + K1*K2, #X, 3) + 2*B2*K1*g*m1^2*m2 + 2*B2*K2*g*m1^2*m2 - K1*x10*B2^2*m1^2*root(#X^4*m1*m2 + #X^3*B2*m1 + #X^3*B1*m2 + #X^2*K2*m1 + #X^2*K1*m2 + #X^2*K1*m1 + #X^2*B1*B2 + #X*B2*K1 + #X*B1*K2 + #X*B1*K1 + K1*K2, #X, 4) + K1*x20*B2^2*m1^2*root(#X^4*m1*m2 + #X^3*B2*m1 + #X^3*B1*m2 + #X^2*K2*m1 + #X^2*K1*m2 + #X^2*K1*m1 + #X^2*B1*B2 + #X*B2*K1 + #X*B1*K2 + #X*B1*K1 + K1*K2, #X, 4) + K2*x20*B2^2*m1^2*root(#X^4*m1*m2 + #X^3*B2*m1 + #X^3*B1*m2 + #X^2*K2*m1 + #X^2*K1*m2 + #X^2*K1*m1 + #X^2*B1*B2 + #X*B2*K1 + #X*B1*K2 + #X*B1*K1 + K1*K2, #X, 4) - g*m2*B2^2*m1^2*root(#X^4*m1*m2 + #X^3*B2*m1 + #X^3*B1*m2 + #X^2*K2*m1 + #X^2*K1*m2 + #X^2*K1*m1 + #X^2*B1*B2 + #X*B2*K1 + #X*B1*K2 + #X*B1*K1 + K1*K2, #X, 1) - g*m2*B2^2*m1^2*root(#X^4*m1*m2 + #X^3*B2*m1 + #X^3*B1*m2 + #X^2*K2*m1 + #X^2*K1*m2 + #X^2*K1*m1 + #X^2*B1*B2 + #X*B2*K1 + #X*B1*K2 + #X*B1*K1 + K1*K2, #X, 2) - g*m2*B2^2*m1^2*root(#X^4*m1*m2 + #X^3*B2*m1 + #X^3*B1*m2 + #X^2*K2*m1 + #X^2*K1*m2 + #X^2*K1*m1 + #X^2*B1*B2 + #X*B2*K1 + #X*B1*K2 + #X*B1*K1 + K1*K2, #X, 3) + K1*g*m1^2*m2^2*root(#X^4*m1*m2 + #X^3*B2*m1 + #X^3*B1*m2 + #X^2*K2*m1 + #X^2*K1*m2 + #X^2*K1*m1 + #X^2*B1*B2 + #X*B2*K1 + #X*B1*K2 + #X*B1*K1 + K1*K2, #X, 1) + K1*g*m1^2*m2^2*root(#X^4*m1*m2 + #X^3*B2*m1 + #X^3*B1*m2 + #X^2*K2*m1 + #X^2*K1*m2 + #X^2*K1*m1 + #X^2*B1*B2 + #X*B2*K1 + #X*B1*K2 + #X*B1*K1 + K1*K2, #X, 2) + K2*g*m1^2*m2^2*root(#X^4*m1*m2 + #X^3*B2*m1 + #X^3*B1*m2 + #X^2*K2*m1 + #X^2*K1*m2 + #X^2*K1*m1 + #X^2*B1*B2 + #X*B2*K1 + #X*B1*K2 + #X*B1*K1 + K1*K2, #X, 1) + K1*g*m1^2*m2^2*root(#X^4*m1*m2 + #X^3*B2*m1 + #X^3*B1*m2 + #X^2*K2*m1 + #X^2*K1*m2 + #X^2*K1*m1 + #X^2*B1*B2 + #X*B2*K1 + #X*B1*K2 + #X*B1*K1 + K1*K2, #X, 3) + K2*g*m1^2*m2^2*root(#X^4*m1*m2 + #X^3*B2*m1 + #X^3*B1*m2 + #X^2*K2*m1 + #X^2*K1*m2 + #X^2*K1*m1 + #X^2*B1*B2 + #X*B2*K1 + #X*B1*K2 + #X*B1*K1 + K1*K2, #X, 2) + K2*g*m1^2*m2^2*root(#X^4*m1*m2 + #X^3*B2*m1 + #X^3*B1*m2 + #X^2*K2*m1 + #X^2*K1*m2 + #X^2*K1*m1 + #X^2*B1*B2 + #X*B2*K1 + #X*B1*K2 + #X*B1*K1 + K1*K2, #X, 3) + m1*x10*K1^2*m2^2*root(#X^4*m1*m2 + #X^3*B2*m1 + #X^3*B1*m2 + #X^2*K2*m1 + #X^2*K1*m2 + #X^2*K1*m1 + #X^2*B1*B2 + #X*B2*K1 + #X*B1*K2 + #X*B1*K1 + K1*K2, #X, 4) + m2*x10*K1^2*m1^2*root(#X^4*m1*m2 + #X^3*B2*m1 + #X^3*B1*m2 + #X^2*K2*m1 + #X^2*K1*m2 + #X^2*K1*m1 + #X^2*B1*B2 + #X*B2*K1 + #X*B1*K2 + #X*B1*K1 + K1*K2, #X, 4) - m1*x20*K1^2*m2^2*root(#X^4*m1*m2 + #X^3*B2*m1 + #X^3*B1*m2 + #X^2*K2*m1 + #X^2*K1*m2 + #X^2*K1*m1 + #X^2*B1*B2 + #X*B2*K1 + #X*B1*K2 + #X*B1*K1 + K1*K2, #X, 4) - m2*x20*K1^2*m1^2*root(#X^4*m1*m2 + #X^3*B2*m1 + #X^3*B1*m2 + #X^2*K2*m1 + #X^2*K1*m2 + #X^2*K1*m1 + #X^2*B1*B2 + #X*B2*K1 + #X*B1*K2 + #X*B1*K1 + K1*K2, #X, 4) - m2*x20*K2^2*m1^2*root(#X^4*m1*m2 + #X^3*B2*m1 + #X^3*B1*m2 + #X^2*K2*m1 + #X^2*K1*m2 + #X^2*K1*m1 + #X^2*B1*B2 + #X*B2*K1 + #X*B1*K2 + #X*B1*K1 + K1*K2, #X, 4) - g*m1^2*m2^3*root(#X^4*m1*m2 + #X^3*B2*m1 + #X^3*B1*m2 + #X^2*K2*m1 + #X^2*K1*m2 + #X^2*K1*m1 + #X^2*B1*B2 + #X*B2*K1 + #X*B1*K2 + #X*B1*K1 + K1*K2, #X, 1)*root(#X^4*m1*m2 + #X^3*B2*m1 + #X^3*B1*m2 + #X^2*K2*m1 + #X^2*K1*m2 + #X^2*K1*m1 + #X^2*B1*B2 + #X*B2*K1 + #X*B1*K2 + #X*B1*K1 + K1*K2, #X, 2)*root(#X^4*m1*m2 + #X^3*B2*m1 + #X^3*B1*m2 + #X^2*K2*m1 + #X^2*K1*m2 + #X^2*K1*m1 + #X^2*B1*B2 + #X*B2*K1 + #X*B1*K2 + #X*B1*K1 + K1*K2, #X, 3) + B1*K1*m1*v10*m2^2*root(#X^4*m1*m2 + #X^3*B2*m1 + #X^3*B1*m2 + #X^2*K2*m1 + #X^2*K1*m2 + #X^2*K1*m1 + #X^2*B1*B2 + #X*B2*K1 + #X*B1*K2 + #X*B1*K1 + K1*K2, #X, 4) + B2*K1*m2*v10*m1^2*root(#X^4*m1*m2 + #X^3*B2*m1 + #X^3*B1*m2 + #X^2*K2*m1 + #X^2*K1*m2 + #X^2*K1*m1 + #X^2*B1*B2 + #X*B2*K1 + #X*B1*K2 + #X*B1*K1 + K1*K2, #X, 4) - 2*B2*K1*m2*v20*m1^2*root(#X^4*m1*m2 + #X^3*B2*m1 + #X^3*B1*m2 + #X^2*K2*m1 + #X^2*K1*m2 + #X^2*K1*m1 + #X^2*B1*B2 + #X*B2*K1 + #X*B1*K2 + #X*B1*K1 + K1*K2, #X, 4) - 2*B2*K2*m2*v20*m1^2*root(#X^4*m1*m2 + #X^3*B2*m1 + #X^3*B1*m2 + #X^2*K2*m1 + #X^2*K1*m2 + #X^2*K1*m1 + #X^2*B1*B2 + #X*B2*K1 + #X*B1*K2 + #X*B1*K1 + K1*K2, #X, 4) + K1*K2*m2*x10*m1^2*root(#X^4*m1*m2 + #X^3*B2*m1 + #X^3*B1*m2 + #X^2*K2*m1 + #X^2*K1*m2 + #X^2*K1*m1 + #X^2*B1*B2 + #X*B2*K1 + #X*B1*K2 + #X*B1*K1 + K1*K2, #X, 4) - 2*K1*K2*m2*x20*m1^2*root(#X^4*m1*m2 + #X^3*B2*m1 + #X^3*B1*m2 + #X^2*K2*m1 + #X^2*K1*m2 + #X^2*K1*m1 + #X^2*B1*B2 + #X*B2*K1 + #X*B1*K2 + #X*B1*K1 + K1*K2, #X, 4) - B2*g*m1^2*m2^2*root(#X^4*m1*m2 + #X^3*B2*m1 + #X^3*B1*m2 + #X^2*K2*m1 + #X^2*K1*m2 + #X^2*K1*m1 + #X^2*B1*B2 + #X*B2*K1 + #X*B1*K2 + #X*B1*K1 + K1*K2, #X, 1)*root(#X^4*m1*m2 + #X^3*B2*m1 + #X^3*B1*m2 + #X^2*K2*m1 + #X^2*K1*m2 + #X^2*K1*m1 + #X^2*B1*B2 + #X*B2*K1 + #X*B1*K2 + #X*B1*K1 + K1*K2, #X, 2) - B2*g*m1^2*m2^2*root(#X^4*m1*m2 + #X^3*B2*m1 + #X^3*B1*m2 + #X^2*K2*m1 + #X^2*K1*m2 + #X^2*K1*m1 + #X^2*B1*B2 + #X*B2*K1 + #X*B1*K2 + #X*B1*K1 + K1*K2, #X, 1)*root(#X^4*m1*m2 + #X^3*B2*m1 + #X^3*B1*m2 + #X^2*K2*m1 + #X^2*K1*m2 + #X^2*K1*m1 + #X^2*B1*B2 + #X*B2*K1 + #X*B1*K2 + #X*B1*K1 + K1*K2, #X, 3) - B2*g*m1^2*m2^2*root(#X^4*m1*m2 + #X^3*B2*m1 + #X^3*B1*m2 + #X^2*K2*m1 + #X^2*K1*m2 + #X^2*K1*m1 + #X^2*B1*B2 + #X*B2*K1 + #X*B1*K2 + #X*B1*K1 + K1*K2, #X, 2)*root(#X^4*m1*m2 + #X^3*B2*m1 + #X^3*B1*m2 + #X^2*K2*m1 + #X^2*K1*m2 + #X^2*K1*m1 + #X^2*B1*B2 + #X*B2*K1 + #X*B1*K2 + #X*B1*K1 + K1*K2, #X, 3) + m2*v20*B2^2*m1^2*root(#X^4*m1*m2 + #X^3*B2*m1 + #X^3*B1*m2 + #X^2*K2*m1 + #X^2*K1*m2 + #X^2*K1*m1 + #X^2*B1*B2 + #X*B2*K1 + #X*B1*K2 + #X*B1*K1 + K1*K2, #X, 1)*root(#X^4*m1*m2 + #X^3*B2*m1 + #X^3*B1*m2 + #X^2*K2*m1 + #X^2*K1*m2 + #X^2*K1*m1 + #X^2*B1*B2 + #X*B2*K1 + #X*B1*K2 + #X*B1*K1 + K1*K2, #X, 4) + m2*v20*B2^2*m1^2*root(#X^4*m1*m2 + #X^3*B2*m1 + #X^3*B1*m2 + #X^2*K2*m1 + #X^2*K1*m2 + #X^2*K1*m1 + #X^2*B1*B2 + #X*B2*K1 + #X*B1*K2 + #X*B1*K1 + K1*K2, #X, 2)*root(#X^4*m1*m2 + #X^3*B2*m1 + #X^3*B1*m2 + #X^2*K2*m1 + #X^2*K1*m2 + #X^2*K1*m1 + #X^2*B1*B2 + #X*B2*K1 + #X*B1*K2 + #X*B1*K1 + K1*K2, #X, 4) + m2*v20*B2^2*m1^2*root(#X^4*m1*m2 + #X^3*B2*m1 + #X^3*B1*m2 + #X^2*K2*m1 + #X^2*K1*m2 + #X^2*K1*m1 + #X^2*B1*B2 + #X*B2*K1 + #X*B1*K2 + #X*B1*K1 + K1*K2, #X, 3)*root(#X^4*m1*m2 + #X^3*B2*m1 + #X^3*B1*m2 + #X^2*K2*m1 + #X^2*K1*m2 + #X^2*K1*m1 + #X^2*B1*B2 + #X*B2*K1 + #X*B1*K2 + #X*B1*K1 + K1*K2, #X, 4) + K1*v10*m1^2*m2^2*root(#X^4*m1*m2 + #X^3*B2*m1 + #X^3*B1*m2 + #X^2*K2*m1 + #X^2*K1*m2 + #X^2*K1*m1 + #X^2*B1*B2 + #X*B2*K1 + #X*B1*K2 + #X*B1*K1 + K1*K2, #X, 1)*root(#X^4*m1*m2 + #X^3*B2*m1 + #X^3*B1*m2 + #X^2*K2*m1 + #X^2*K1*m2 + #X^2*K1*m1 + #X^2*B1*B2 + #X*B2*K1 + #X*B1*K2 + #X*B1*K1 + K1*K2, #X, 4) + K1*v10*m1^2*m2^2*root(#X^4*m1*m2 + #X^3*B2*m1 + #X^3*B1*m2 + #X^2*K2*m1 + #X^2*K1*m2 + #X^2*K1*m1 + #X^2*B1*B2 + #X*B2*K1 + #X*B1*K2 + #X*B1*K1 + K1*K2, #X, 2)*root(#X^4*m1*m2 + #X^3*B2*m1 + #X^3*B1*m2 + #X^2*K2*m1 + #X^2*K1*m2 + #X^2*K1*m1 + #X^2*B1*B2 + #X*B2*K1 + #X*B1*K2 + #X*B1*K1 + K1*K2, #X, 4) + K1*v10*m1^2*m2^2*root(#X^4*m1*m2 + #X^3*B2*m1 + #X^3*B1*m2 + #X^2*K2*m1 + #X^2*K1*m2 + #X^2*K1*m1 + #X^2*B1*B2 + #X*B2*K1 + #X*B1*K2 + #X*B1*K1 + K1*K2, #X, 3)*root(#X^4*m1*m2 + #X^3*B2*m1 + #X^3*B1*m2 + #X^2*K2*m1 + #X^2*K1*m2 + #X^2*K1*m1 + #X^2*B1*B2 + #X*B2*K1 + #X*B1*K2 + #X*B1*K1 + K1*K2, #X, 4) - K1*v20*m1^2*m2^2*root(#X^4*m1*m2 + #X^3*B2*m1 + #X^3*B1*m2 + #X^2*K2*m1 + #X^2*K1*m2 + #X^2*K1*m1 + #X^2*B1*B2 + #X*B2*K1 + #X*B1*K2 + #X*B1*K1 + K1*K2, #X, 1)*root(#X^4*m1*m2 + #X^3*B2*m1 + #X^3*B1*m2 + #X^2*K2*m1 + #X^2*K1*m2 + #X^2*K1*m1 + #X^2*B1*B2 + #X*B2*K1 + #X*B1*K2 + #X*B1*K1 + K1*K2, #X, 4) - K1*v20*m1^2*m2^2*root(#X^4*m1*m2 + #X^3*B2*m1 + #X^3*B1*m2 + #X^2*K2*m1 + #X^2*K1*m2 + #X^2*K1*m1 + #X^2*B1*B2 + #X*B2*K1 + #X*B1*K2 + #X*B1*K1 + K1*K2, #X, 2)*root(#X^4*m1*m2 + #X^3*B2*m1 + #X^3*B1*m2 + #X^2*K2*m1 + #X^2*K1*m2 + #X^2*K1*m1 + #X^2*B1*B2 + #X*B2*K1 + #X*B1*K2 + #X*B1*K1 + K1*K2, #X, 4) - K2*v20*m1^2*m2^2*root(#X^4*m1*m2 + #X^3*B2*m1 + #X^3*B1*m2 + #X^2*K2*m1 + #X^2*K1*m2 + #X^2*K1*m1 + #X^2*B1*B2 + #X*B2*K1 + #X*B1*K2 + #X*B1*K1 + K1*K2, #X, 1)*root(#X^4*m1*m2 + #X^3*B2*m1 + #X^3*B1*m2 + #X^2*K2*m1 + #X^2*K1*m2 + #X^2*K1*m1 + #X^2*B1*B2 + #X*B2*K1 + #X*B1*K2 + #X*B1*K1 + K1*K2, #X, 4) - K1*v20*m1^2*m2^2*root(#X^4*m1*m2 + #X^3*B2*m1 + #X^3*B1*m2 + #X^2*K2*m1 + #X^2*K1*m2 + #X^2*K1*m1 + #X^2*B1*B2 + #X*B2*K1 + #X*B1*K2 + #X*B1*K1 + K1*K2, #X, 3)*root(#X^4*m1*m2 + #X^3*B2*m1 + #X^3*B1*m2 + #X^2*K2*m1 + #X^2*K1*m2 + #X^2*K1*m1 + #X^2*B1*B2 + #X*B2*K1 + #X*B1*K2 + #X*B1*K1 + K1*K2, #X, 4) - K2*v20*m1^2*m2^2*root(#X^4*m1*m2 + #X^3*B2*m1 + #X^3*B1*m2 + #X^2*K2*m1 + #X^2*K1*m2 + #X^2*K1*m1 + #X^2*B1*B2 + #X*B2*K1 + #X*B1*K2 + #X*B1*K1 + K1*K2, #X, 2)*root(#X^4*m1*m2 + #X^3*B2*m1 + #X^3*B1*m2 + #X^2*K2*m1 + #X^2*K1*m2 + #X^2*K1*m1 + #X^2*B1*B2 + #X*B2*K1 + #X*B1*K2 + #X*B1*K1 + K1*K2, #X, 4) - K2*v20*m1^2*m2^2*root(#X^4*m1*m2 + #X^3*B2*m1 + #X^3*B1*m2 + #X^2*K2*m1 + #X^2*K1*m2 + #X^2*K1*m1 + #X^2*B1*B2 + #X*B2*K1 + #X*B1*K2 + #X*B1*K1 + K1*K2, #X, 3)*root(#X^4*m1*m2 + #X^3*B2*m1 + #X^3*B1*m2 + #X^2*K2*m1 + #X^2*K1*m2 + #X^2*K1*m1 + #X^2*B1*B2 + #X*B2*K1 + #X*B1*K2 + #X*B1*K1 + K1*K2, #X, 4) + v20*m1^2*m2^3*root(#X^4*m1*m2 + #X^3*B2*m1 + #X^3*B1*m2 + #X^2*K2*m1 + #X^2*K1*m2 + #X^2*K1*m1 + #X^2*B1*B2 + #X*B2*K1 + #X*B1*K2 + #X*B1*K1 + K1*K2, #X, 1)*root(#X^4*m1*m2 + #X^3*B2*m1 + #X^3*B1*m2 + #X^2*K2*m1 + #X^2*K1*m2 + #X^2*K1*m1 + #X^2*B1*B2 + #X*B2*K1 + #X*B1*K2 + #X*B1*K1 + K1*K2, #X, 2)*root(#X^4*m1*m2 + #X^3*B2*m1 + #X^3*B1*m2 + #X^2*K2*m1 + #X^2*K1*m2 + #X^2*K1*m1 + #X^2*B1*B2 + #X*B2*K1 + #X*B1*K2 + #X*B1*K1 + K1*K2, #X, 3)*root(#X^4*m1*m2 + #X^3*B2*m1 + #X^3*B1*m2 + #X^2*K2*m1 + #X^2*K1*m2 + #X^2*K1*m1 + #X^2*B1*B2 + #X*B2*K1 + #X*B1*K2 + #X*B1*K1 + K1*K2, #X, 4) - B2*K1*m2*x10*m1^2*root(#X^4*m1*m2 + #X^3*B2*m1 + #X^3*B1*m2 + #X^2*K2*m1 + #X^2*K1*m2 + #X^2*K1*m1 + #X^2*B1*B2 + #X*B2*K1 + #X*B1*K2 + #X*B1*K1 + K1*K2, #X, 1)*root(#X^4*m1*m2 + #X^3*B2*m1 + #X^3*B1*m2 + #X^2*K2*m1 + #X^2*K1*m2 + #X^2*K1*m1 + #X^2*B1*B2 + #X*B2*K1 + #X*B1*K2 + #X*B1*K1 + K1*K2, #X, 4) - B2*K1*m2*x10*m1^2*root(#X^4*m1*m2 + #X^3*B2*m1 + #X^3*B1*m2 + #X^2*K2*m1 + #X^2*K1*m2 + #X^2*K1*m1 + #X^2*B1*B2 + #X*B2*K1 + #X*B1*K2 + #X*B1*K1 + K1*K2, #X, 2)*root(#X^4*m1*m2 + #X^3*B2*m1 + #X^3*B1*m2 + #X^2*K2*m1 + #X^2*K1*m2 + #X^2*K1*m1 + #X^2*B1*B2 + #X*B2*K1 + #X*B1*K2 + #X*B1*K1 + K1*K2, #X, 4) - B2*K1*m2*x10*m1^2*root(#X^4*m1*m2 + #X^3*B2*m1 + #X^3*B1*m2 + #X^2*K2*m1 + #X^2*K1*m2 + #X^2*K1*m1 + #X^2*B1*B2 + #X*B2*K1 + #X*B1*K2 + #X*B1*K1 + K1*K2, #X, 3)*root(#X^4*m1*m2 + #X^3*B2*m1 + #X^3*B1*m2 + #X^2*K2*m1 + #X^2*K1*m2 + #X^2*K1*m1 + #X^2*B1*B2 + #X*B2*K1 + #X*B1*K2 + #X*B1*K1 + K1*K2, #X, 4) + B2*K1*m2*x20*m1^2*root(#X^4*m1*m2 + #X^3*B2*m1 + #X^3*B1*m2 + #X^2*K2*m1 + #X^2*K1*m2 + #X^2*K1*m1 + #X^2*B1*B2 + #X*B2*K1 + #X*B1*K2 + #X*B1*K1 + K1*K2, #X, 1)*root(#X^4*m1*m2 + #X^3*B2*m1 + #X^3*B1*m2 + #X^2*K2*m1 + #X^2*K1*m2 + #X^2*K1*m1 + #X^2*B1*B2 + #X*B2*K1 + #X*B1*K2 + #X*B1*K1 + K1*K2, #X, 4) + B2*K1*m2*x20*m1^2*root(#X^4*m1*m2 + #X^3*B2*m1 + #X^3*B1*m2 + #X^2*K2*m1 + #X^2*K1*m2 + #X^2*K1*m1 + #X^2*B1*B2 + #X*B2*K1 + #X*B1*K2 + #X*B1*K1 + K1*K2, #X, 2)*root(#X^4*m1*m2 + #X^3*B2*m1 + #X^3*B1*m2 + #X^2*K2*m1 + #X^2*K1*m2 + #X^2*K1*m1 + #X^2*B1*B2 + #X*B2*K1 + #X*B1*K2 + #X*B1*K1 + K1*K2, #X, 4) + B2*K2*m2*x20*m1^2*root(#X^4*m1*m2 + #X^3*B2*m1 + #X^3*B1*m2 + #X^2*K2*m1 + #X^2*K1*m2 + #X^2*K1*m1 + #X^2*B1*B2 + #X*B2*K1 + #X*B1*K2 + #X*B1*K1 + K1*K2, #X, 1)*root(#X^4*m1*m2 + #X^3*B2*m1 + #X^3*B1*m2 + #X^2*K2*m1 + #X^2*K1*m2 + #X^2*K1*m1 + #X^2*B1*B2 + #X*B2*K1 + #X*B1*K2 + #X*B1*K1 + K1*K2, #X, 4) + B2*K1*m2*x20*m1^2*root(#X^4*m1*m2 + #X^3*B2*m1 + #X^3*B1*m2 + #X^2*K2*m1 + #X^2*K1*m2 + #X^2*K1*m1 + #X^2*B1*B2 + #X*B2*K1 + #X*B1*K2 + #X*B1*K1 + K1*K2, #X, 3)*root(#X^4*m1*m2 + #X^3*B2*m1 + #X^3*B1*m2 + #X^2*K2*m1 + #X^2*K1*m2 + #X^2*K1*m1 + #X^2*B1*B2 + #X*B2*K1 + #X*B1*K2 + #X*B1*K1 + K1*K2, #X, 4) + B2*K2*m2*x20*m1^2*root(#X^4*m1*m2 + #X^3*B2*m1 + #X^3*B1*m2 + #X^2*K2*m1 + #X^2*K1*m2 + #X^2*K1*m1 + #X^2*B1*B2 + #X*B2*K1 + #X*B1*K2 + #X*B1*K1 + K1*K2, #X, 2)*root(#X^4*m1*m2 + #X^3*B2*m1 + #X^3*B1*m2 + #X^2*K2*m1 + #X^2*K1*m2 + #X^2*K1*m1 + #X^2*B1*B2 + #X*B2*K1 + #X*B1*K2 + #X*B1*K1 + K1*K2, #X, 4) + B2*K2*m2*x20*m1^2*root(#X^4*m1*m2 + #X^3*B2*m1 + #X^3*B1*m2 + #X^2*K2*m1 + #X^2*K1*m2 + #X^2*K1*m1 + #X^2*B1*B2 + #X*B2*K1 + #X*B1*K2 + #X*B1*K1 + K1*K2, #X, 3)*root(#X^4*m1*m2 + #X^3*B2*m1 + #X^3*B1*m2 + #X^2*K2*m1 + #X^2*K1*m2 + #X^2*K1*m1 + #X^2*B1*B2 + #X*B2*K1 + #X*B1*K2 + #X*B1*K1 + K1*K2, #X, 4) + B2*v20*m1^2*m2^2*root(#X^4*m1*m2 + #X^3*B2*m1 + #X^3*B1*m2 + #X^2*K2*m1 + #X^2*K1*m2 + #X^2*K1*m1 + #X^2*B1*B2 + #X*B2*K1 + #X*B1*K2 + #X*B1*K1 + K1*K2, #X, 1)*root(#X^4*m1*m2 + #X^3*B2*m1 + #X^3*B1*m2 + #X^2*K2*m1 + #X^2*K1*m2 + #X^2*K1*m1 + #X^2*B1*B2 + #X*B2*K1 + #X*B1*K2 + #X*B1*K1 + K1*K2, #X, 2)*root(#X^4*m1*m2 + #X^3*B2*m1 + #X^3*B1*m2 + #X^2*K2*m1 + #X^2*K1*m2 + #X^2*K1*m1 + #X^2*B1*B2 + #X*B2*K1 + #X*B1*K2 + #X*B1*K1 + K1*K2, #X, 4) + B2*v20*m1^2*m2^2*root(#X^4*m1*m2 + #X^3*B2*m1 + #X^3*B1*m2 + #X^2*K2*m1 + #X^2*K1*m2 + #X^2*K1*m1 + #X^2*B1*B2 + #X*B2*K1 + #X*B1*K2 + #X*B1*K1 + K1*K2, #X, 1)*root(#X^4*m1*m2 + #X^3*B2*m1 + #X^3*B1*m2 + #X^2*K2*m1 + #X^2*K1*m2 + #X^2*K1*m1 + #X^2*B1*B2 + #X*B2*K1 + #X*B1*K2 + #X*B1*K1 + K1*K2, #X, 3)*root(#X^4*m1*m2 + #X^3*B2*m1 + #X^3*B1*m2 + #X^2*K2*m1 + #X^2*K1*m2 + #X^2*K1*m1 + #X^2*B1*B2 + #X*B2*K1 + #X*B1*K2 + #X*B1*K1 + K1*K2, #X, 4) + B2*v20*m1^2*m2^2*root(#X^4*m1*m2 + #X^3*B2*m1 + #X^3*B1*m2 + #X^2*K2*m1 + #X^2*K1*m2 + #X^2*K1*m1 + #X^2*B1*B2 + #X*B2*K1 + #X*B1*K2 + #X*B1*K1 + K1*K2, #X, 2)*root(#X^4*m1*m2 + #X^3*B2*m1 + #X^3*B1*m2 + #X^2*K2*m1 + #X^2*K1*m2 + #X^2*K1*m1 + #X^2*B1*B2 + #X*B2*K1 + #X*B1*K2 + #X*B1*K1 + K1*K2, #X, 3)*root(#X^4*m1*m2 + #X^3*B2*m1 + #X^3*B1*m2 + #X^2*K2*m1 + #X^2*K1*m2 + #X^2*K1*m1 + #X^2*B1*B2 + #X*B2*K1 + #X*B1*K2 + #X*B1*K1 + K1*K2, #X, 4) - K1*x10*m1^2*m2^2*root(#X^4*m1*m2 + #X^3*B2*m1 + #X^3*B1*m2 + #X^2*K2*m1 + #X^2*K1*m2 + #X^2*K1*m1 + #X^2*B1*B2 + #X*B2*K1 + #X*B1*K2 + #X*B1*K1 + K1*K2, #X, 1)*root(#X^4*m1*m2 + #X^3*B2*m1 + #X^3*B1*m2 + #X^2*K2*m1 + #X^2*K1*m2 + #X^2*K1*m1 + #X^2*B1*B2 + #X*B2*K1 + #X*B1*K2 + #X*B1*K1 + K1*K2, #X, 2)*root(#X^4*m1*m2 + #X^3*B2*m1 + #X^3*B1*m2 + #X^2*K2*m1 + #X^2*K1*m2 + #X^2*K1*m1 + #X^2*B1*B2 + #X*B2*K1 + #X*B1*K2 + #X*B1*K1 + K1*K2, #X, 4) - K1*x10*m1^2*m2^2*root(#X^4*m1*m2 + #X^3*B2*m1 + #X^3*B1*m2 + #X^2*K2*m1 + #X^2*K1*m2 + #X^2*K1*m1 + #X^2*B1*B2 + #X*B2*K1 + #X*B1*K2 + #X*B1*K1 + K1*K2, #X, 1)*root(#X^4*m1*m2 + #X^3*B2*m1 + #X^3*B1*m2 + #X^2*K2*m1 + #X^2*K1*m2 + #X^2*K1*m1 + #X^2*B1*B2 + #X*B2*K1 + #X*B1*K2 + #X*B1*K1 + K1*K2, #X, 3)*root(#X^4*m1*m2 + #X^3*B2*m1 + #X^3*B1*m2 + #X^2*K2*m1 + #X^2*K1*m2 + #X^2*K1*m1 + #X^2*B1*B2 + #X*B2*K1 + #X*B1*K2 + #X*B1*K1 + K1*K2, #X, 4) - K1*x10*m1^2*m2^2*root(#X^4*m1*m2 + #X^3*B2*m1 + #X^3*B1*m2 + #X^2*K2*m1 + #X^2*K1*m2 + #X^2*K1*m1 + #X^2*B1*B2 + #X*B2*K1 + #X*B1*K2 + #X*B1*K1 + K1*K2, #X, 2)*root(#X^4*m1*m2 + #X^3*B2*m1 + #X^3*B1*m2 + #X^2*K2*m1 + #X^2*K1*m2 + #X^2*K1*m1 + #X^2*B1*B2 + #X*B2*K1 + #X*B1*K2 + #X*B1*K1 + K1*K2, #X, 3)*root(#X^4*m1*m2 + #X^3*B2*m1 + #X^3*B1*m2 + #X^2*K2*m1 + #X^2*K1*m2 + #X^2*K1*m1 + #X^2*B1*B2 + #X*B2*K1 + #X*B1*K2 + #X*B1*K1 + K1*K2, #X, 4) + K1*x20*m1^2*m2^2*root(#X^4*m1*m2 + #X^3*B2*m1 + #X^3*B1*m2 + #X^2*K2*m1 + #X^2*K1*m2 + #X^2*K1*m1 + #X^2*B1*B2 + #X*B2*K1 + #X*B1*K2 + #X*B1*K1 + K1*K2, #X, 1)*root(#X^4*m1*m2 + #X^3*B2*m1 + #X^3*B1*m2 + #X^2*K2*m1 + #X^2*K1*m2 + #X^2*K1*m1 + #X^2*B1*B2 + #X*B2*K1 + #X*B1*K2 + #X*B1*K1 + K1*K2, #X, 2)*root(#X^4*m1*m2 + #X^3*B2*m1 + #X^3*B1*m2 + #X^2*K2*m1 + #X^2*K1*m2 + #X^2*K1*m1 + #X^2*B1*B2 + #X*B2*K1 + #X*B1*K2 + #X*B1*K1 + K1*K2, #X, 4) + K1*x20*m1^2*m2^2*root(#X^4*m1*m2 + #X^3*B2*m1 + #X^3*B1*m2 + #X^2*K2*m1 + #X^2*K1*m2 + #X^2*K1*m1 + #X^2*B1*B2 + #X*B2*K1 + #X*B1*K2 + #X*B1*K1 + K1*K2, #X, 1)*root(#X^4*m1*m2 + #X^3*B2*m1 + #X^3*B1*m2 + #X^2*K2*m1 + #X^2*K1*m2 + #X^2*K1*m1 + #X^2*B1*B2 + #X*B2*K1 + #X*B1*K2 + #X*B1*K1 + K1*K2, #X, 3)*root(#X^4*m1*m2 + #X^3*B2*m1 + #X^3*B1*m2 + #X^2*K2*m1 + #X^2*K1*m2 + #X^2*K1*m1 + #X^2*B1*B2 + #X*B2*K1 + #X*B1*K2 + #X*B1*K1 + K1*K2, #X, 4) + K2*x20*m1^2*m2^2*root(#X^4*m1*m2 + #X^3*B2*m1 + #X^3*B1*m2 + #X^2*K2*m1 + #X^2*K1*m2 + #X^2*K1*m1 + #X^2*B1*B2 + #X*B2*K1 + #X*B1*K2 + #X*B1*K1 + K1*K2, #X, 1)*root(#X^4*m1*m2 + #X^3*B2*m1 + #X^3*B1*m2 + #X^2*K2*m1 + #X^2*K1*m2 + #X^2*K1*m1 + #X^2*B1*B2 + #X*B2*K1 + #X*B1*K2 + #X*B1*K1 + K1*K2, #X, 2)*root(#X^4*m1*m2 + #X^3*B2*m1 + #X^3*B1*m2 + #X^2*K2*m1 + #X^2*K1*m2 + #X^2*K1*m1 + #X^2*B1*B2 + #X*B2*K1 + #X*B1*K2 + #X*B1*K1 + K1*K2, #X, 4) + K1*x20*m1^2*m2^2*root(#X^4*m1*m2 + #X^3*B2*m1 + #X^3*B1*m2 + #X^2*K2*m1 + #X^2*K1*m2 + #X^2*K1*m1 + #X^2*B1*B2 + #X*B2*K1 + #X*B1*K2 + #X*B1*K1 + K1*K2, #X, 2)*root(#X^4*m1*m2 + #X^3*B2*m1 + #X^3*B1*m2 + #X^2*K2*m1 + #X^2*K1*m2 + #X^2*K1*m1 + #X^2*B1*B2 + #X*B2*K1 + #X*B1*K2 + #X*B1*K1 + K1*K2, #X, 3)*root(#X^4*m1*m2 + #X^3*B2*m1 + #X^3*B1*m2 + #X^2*K2*m1 + #X^2*K1*m2 + #X^2*K1*m1 + #X^2*B1*B2 + #X*B2*K1 + #X*B1*K2 + #X*B1*K1 + K1*K2, #X, 4) + K2*x20*m1^2*m2^2*root(#X^4*m1*m2 + #X^3*B2*m1 + #X^3*B1*m2 + #X^2*K2*m1 + #X^2*K1*m2 + #X^2*K1*m1 + #X^2*B1*B2 + #X*B2*K1 + #X*B1*K2 + #X*B1*K1 + K1*K2, #X, 1)*root(#X^4*m1*m2 + #X^3*B2*m1 + #X^3*B1*m2 + #X^2*K2*m1 + #X^2*K1*m2 + #X^2*K1*m1 + #X^2*B1*B2 + #X*B2*K1 + #X*B1*K2 + #X*B1*K1 + K1*K2, #X, 3)*root(#X^4*m1*m2 + #X^3*B2*m1 + #X^3*B1*m2 + #X^2*K2*m1 + #X^2*K1*m2 + #X^2*K1*m1 + #X^2*B1*B2 + #X*B2*K1 + #X*B1*K2 + #X*B1*K1 + K1*K2, #X, 4) + K2*x20*m1^2*m2^2*root(#X^4*m1*m2 + #X^3*B2*m1 + #X^3*B1*m2 + #X^2*K2*m1 + #X^2*K1*m2 + #X^2*K1*m1 + #X^2*B1*B2 + #X*B2*K1 + #X*B1*K2 + #X*B1*K1 + K1*K2, #X, 2)*root(#X^4*m1*m2 + #X^3*B2*m1 + #X^3*B1*m2 + #X^2*K2*m1 + #X^2*K1*m2 + #X^2*K1*m1 + #X^2*B1*B2 + #X*B2*K1 + #X*B1*K2 + #X*B1*K1 + K1*K2, #X, 3)*root(#X^4*m1*m2 + #X^3*B2*m1 + #X^3*B1*m2 + #X^2*K2*m1 + #X^2*K1*m2 + #X^2*K1*m1 + #X^2*B1*B2 + #X*B2*K1 + #X*B1*K2 + #X*B1*K1 + K1*K2, #X, 4))/(m1*root(#X^4*m1*m2 + #X^3*B2*m1 + #X^3*B1*m2 + #X^2*K2*m1 + #X^2*K1*m2 + #X^2*K1*m1 + #X^2*B1*B2 + #X*B2*K1 + #X*B1*K2 + #X*B1*K1 + K1*K2, #X, 4)*(m1*m2^3*root(#X^4*m1*m2 + #X^3*B2*m1 + #X^3*B1*m2 + #X^2*K2*m1 + #X^2*K1*m2 + #X^2*K1*m1 + #X^2*B1*B2 + #X*B2*K1 + #X*B1*K2 + #X*B1*K1 + K1*K2, #X, 4)^3 - m1*root(#X^4*m1*m2 + #X^3*B2*m1 + #X^3*B1*m2 + #X^2*K2*m1 + #X^2*K1*m2 + #X^2*K1*m1 + #X^2*B1*B2 + #X*B2*K1 + #X*B1*K2 + #X*B1*K1 + K1*K2, #X, 4)^2*m2^3*root(#X^4*m1*m2 + #X^3*B2*m1 + #X^3*B1*m2 + #X^2*K2*m1 + #X^2*K1*m2 + #X^2*K1*m1 + #X^2*B1*B2 + #X*B2*K1 + #X*B1*K2 + #X*B1*K1 + K1*K2, #X, 1) - m1*root(#X^4*m1*m2 + #X^3*B2*m1 + #X^3*B1*m2 + #X^2*K2*m1 + #X^2*K1*m2 + #X^2*K1*m1 + #X^2*B1*B2 + #X*B2*K1 + #X*B1*K2 + #X*B1*K1 + K1*K2, #X, 4)^2*m2^3*root(#X^4*m1*m2 + #X^3*B2*m1 + #X^3*B1*m2 + #X^2*K2*m1 + #X^2*K1*m2 + #X^2*K1*m1 + #X^2*B1*B2 + #X*B2*K1 + #X*B1*K2 + #X*B1*K1 + K1*K2, #X, 2) - m1*root(#X^4*m1*m2 + #X^3*B2*m1 + #X^3*B1*m2 + #X^2*K2*m1 + #X^2*K1*m2 + #X^2*K1*m1 + #X^2*B1*B2 + #X*B2*K1 + #X*B1*K2 + #X*B1*K1 + K1*K2, #X, 4)^2*m2^3*root(#X^4*m1*m2 + #X^3*B2*m1 + #X^3*B1*m2 + #X^2*K2*m1 + #X^2*K1*m2 + #X^2*K1*m1 + #X^2*B1*B2 + #X*B2*K1 + #X*B1*K2 + #X*B1*K1 + K1*K2, #X, 3) - m1*m2^3*root(#X^4*m1*m2 + #X^3*B2*m1 + #X^3*B1*m2 + #X^2*K2*m1 + #X^2*K1*m2 + #X^2*K1*m1 + #X^2*B1*B2 + #X*B2*K1 + #X*B1*K2 + #X*B1*K1 + K1*K2, #X, 1)*root(#X^4*m1*m2 + #X^3*B2*m1 + #X^3*B1*m2 + #X^2*K2*m1 + #X^2*K1*m2 + #X^2*K1*m1 + #X^2*B1*B2 + #X*B2*K1 + #X*B1*K2 + #X*B1*K1 + K1*K2, #X, 2)*root(#X^4*m1*m2 + #X^3*B2*m1 + #X^3*B1*m2 + #X^2*K2*m1 + #X^2*K1*m2 + #X^2*K1*m1 + #X^2*B1*B2 + #X*B2*K1 + #X*B1*K2 + #X*B1*K1 + K1*K2, #X, 3) + m1*m2^3*root(#X^4*m1*m2 + #X^3*B2*m1 + #X^3*B1*m2 + #X^2*K2*m1 + #X^2*K1*m2 + #X^2*K1*m1 + #X^2*B1*B2 + #X*B2*K1 + #X*B1*K2 + #X*B1*K1 + K1*K2, #X, 1)*root(#X^4*m1*m2 + #X^3*B2*m1 + #X^3*B1*m2 + #X^2*K2*m1 + #X^2*K1*m2 + #X^2*K1*m1 + #X^2*B1*B2 + #X*B2*K1 + #X*B1*K2 + #X*B1*K1 + K1*K2, #X, 2)*root(#X^4*m1*m2 + #X^3*B2*m1 + #X^3*B1*m2 + #X^2*K2*m1 + #X^2*K1*m2 + #X^2*K1*m1 + #X^2*B1*B2 + #X*B2*K1 + #X*B1*K2 + #X*B1*K1 + K1*K2, #X, 4) + m1*m2^3*root(#X^4*m1*m2 + #X^3*B2*m1 + #X^3*B1*m2 + #X^2*K2*m1 + #X^2*K1*m2 + #X^2*K1*m1 + #X^2*B1*B2 + #X*B2*K1 + #X*B1*K2 + #X*B1*K1 + K1*K2, #X, 1)*root(#X^4*m1*m2 + #X^3*B2*m1 + #X^3*B1*m2 + #X^2*K2*m1 + #X^2*K1*m2 + #X^2*K1*m1 + #X^2*B1*B2 + #X*B2*K1 + #X*B1*K2 + #X*B1*K1 + K1*K2, #X, 3)*root(#X^4*m1*m2 + #X^3*B2*m1 + #X^3*B1*m2 + #X^2*K2*m1 + #X^2*K1*m2 + #X^2*K1*m1 + #X^2*B1*B2 + #X*B2*K1 + #X*B1*K2 + #X*B1*K1 + K1*K2, #X, 4) + m1*m2^3*root(#X^4*m1*m2 + #X^3*B2*m1 + #X^3*B1*m2 + #X^2*K2*m1 + #X^2*K1*m2 + #X^2*K1*m1 + #X^2*B1*B2 + #X*B2*K1 + #X*B1*K2 + #X*B1*K1 + K1*K2, #X, 2)*root(#X^4*m1*m2 + #X^3*B2*m1 + #X^3*B1*m2 + #X^2*K2*m1 + #X^2*K1*m2 + #X^2*K1*m1 + #X^2*B1*B2 + #X*B2*K1 + #X*B1*K2 + #X*B1*K1 + K1*K2, #X, 3)*root(#X^4*m1*m2 + #X^3*B2*m1 + #X^3*B1*m2 + #X^2*K2*m1 + #X^2*K1*m2 + #X^2*K1*m1 + #X^2*B1*B2 + #X*B2*K1 + #X*B1*K2 + #X*B1*K1 + K1*K2, #X, 4))))*((K1 + K2)/K1 + (B2*root(#X^4*m1*m2 + #X^3*B2*m1 + #X^3*B1*m2 + #X^2*K2*m1 + #X^2*K1*m2 + #X^2*K1*m1 + #X^2*B1*B2 + #X*B2*K1 + #X*B1*K2 + #X*B1*K1 + K1*K2, #X, 4))/K1 + (m2*root(#X^4*m1*m2 + #X^3*B2*m1 + #X^3*B1*m2 + #X^2*K2*m1 + #X^2*K1*m2 + #X^2*K1*m1 + #X^2*B1*B2 + #X*B2*K1 + #X*B1*K2 + #X*B1*K1 + K1*K2, #X, 4)^2)/K1) + exp(t*root(#X^4*m1*m2 + #X^3*B2*m1 + #X^3*B1*m2 + #X^2*K2*m1 + #X^2*K1*m2 + #X^2*K1*m1 + #X^2*B1*B2 + #X*B2*K1 + #X*B1*K2 + #X*B1*K1 + K1*K2, #X, 1))*((exp(-t*root(#X^4*m1*m2 + #X^3*B2*m1 + #X^3*B1*m2 + #X^2*K2*m1 + #X^2*K1*m2 + #X^2*K1*m1 + #X^2*B1*B2 + #X*B2*K1 + #X*B1*K2 + #X*B1*K1 + K1*K2, #X, 1))*(B1*K1*P*m2^2 - B2^3*g*m1^2 + B2*K1*P*m1*m2 + K1*P*m1*m2^2*root(#X^4*m1*m2 + #X^3*B2*m1 + #X^3*B1*m2 + #X^2*K2*m1 + #X^2*K1*m2 + #X^2*K1*m1 + #X^2*B1*B2 + #X*B2*K1 + #X*B1*K2 + #X*B1*K1 + K1*K2, #X, 2) + K1*P*m1*m2^2*root(#X^4*m1*m2 + #X^3*B2*m1 + #X^3*B1*m2 + #X^2*K2*m1 + #X^2*K1*m2 + #X^2*K1*m1 + #X^2*B1*B2 + #X*B2*K1 + #X*B1*K2 + #X*B1*K1 + K1*K2, #X, 3) + K1*P*m1*m2^2*root(#X^4*m1*m2 + #X^3*B2*m1 + #X^3*B1*m2 + #X^2*K2*m1 + #X^2*K1*m2 + #X^2*K1*m1 + #X^2*B1*B2 + #X*B2*K1 + #X*B1*K2 + #X*B1*K1 + K1*K2, #X, 4) + 2*B2*K1*g*m1^2*m2 + 2*B2*K2*g*m1^2*m2 - g*m2*B2^2*m1^2*root(#X^4*m1*m2 + #X^3*B2*m1 + #X^3*B1*m2 + #X^2*K2*m1 + #X^2*K1*m2 + #X^2*K1*m1 + #X^2*B1*B2 + #X*B2*K1 + #X*B1*K2 + #X*B1*K1 + K1*K2, #X, 2) - g*m2*B2^2*m1^2*root(#X^4*m1*m2 + #X^3*B2*m1 + #X^3*B1*m2 + #X^2*K2*m1 + #X^2*K1*m2 + #X^2*K1*m1 + #X^2*B1*B2 + #X*B2*K1 + #X*B1*K2 + #X*B1*K1 + K1*K2, #X, 3) - g*m2*B2^2*m1^2*root(#X^4*m1*m2 + #X^3*B2*m1 + #X^3*B1*m2 + #X^2*K2*m1 + #X^2*K1*m2 + #X^2*K1*m1 + #X^2*B1*B2 + #X*B2*K1 + #X*B1*K2 + #X*B1*K1 + K1*K2, #X, 4) + K1*g*m1^2*m2^2*root(#X^4*m1*m2 + #X^3*B2*m1 + #X^3*B1*m2 + #X^2*K2*m1 + #X^2*K1*m2 + #X^2*K1*m1 + #X^2*B1*B2 + #X*B2*K1 + #X*B1*K2 + #X*B1*K1 + K1*K2, #X, 2) + K1*g*m1^2*m2^2*root(#X^4*m1*m2 + #X^3*B2*m1 + #X^3*B1*m2 + #X^2*K2*m1 + #X^2*K1*m2 + #X^2*K1*m1 + #X^2*B1*B2 + #X*B2*K1 + #X*B1*K2 + #X*B1*K1 + K1*K2, #X, 3) + K2*g*m1^2*m2^2*root(#X^4*m1*m2 + #X^3*B2*m1 + #X^3*B1*m2 + #X^2*K2*m1 + #X^2*K1*m2 + #X^2*K1*m1 + #X^2*B1*B2 + #X*B2*K1 + #X*B1*K2 + #X*B1*K1 + K1*K2, #X, 2) + K1*g*m1^2*m2^2*root(#X^4*m1*m2 + #X^3*B2*m1 + #X^3*B1*m2 + #X^2*K2*m1 + #X^2*K1*m2 + #X^2*K1*m1 + #X^2*B1*B2 + #X*B2*K1 + #X*B1*K2 + #X*B1*K1 + K1*K2, #X, 4) + K2*g*m1^2*m2^2*root(#X^4*m1*m2 + #X^3*B2*m1 + #X^3*B1*m2 + #X^2*K2*m1 + #X^2*K1*m2 + #X^2*K1*m1 + #X^2*B1*B2 + #X*B2*K1 + #X*B1*K2 + #X*B1*K1 + K1*K2, #X, 3) + K2*g*m1^2*m2^2*root(#X^4*m1*m2 + #X^3*B2*m1 + #X^3*B1*m2 + #X^2*K2*m1 + #X^2*K1*m2 + #X^2*K1*m1 + #X^2*B1*B2 + #X*B2*K1 + #X*B1*K2 + #X*B1*K1 + K1*K2, #X, 4) - g*m1^2*m2^3*root(#X^4*m1*m2 + #X^3*B2*m1 + #X^3*B1*m2 + #X^2*K2*m1 + #X^2*K1*m2 + #X^2*K1*m1 + #X^2*B1*B2 + #X*B2*K1 + #X*B1*K2 + #X*B1*K1 + K1*K2, #X, 2)*root(#X^4*m1*m2 + #X^3*B2*m1 + #X^3*B1*m2 + #X^2*K2*m1 + #X^2*K1*m2 + #X^2*K1*m1 + #X^2*B1*B2 + #X*B2*K1 + #X*B1*K2 + #X*B1*K1 + K1*K2, #X, 3)*root(#X^4*m1*m2 + #X^3*B2*m1 + #X^3*B1*m2 + #X^2*K2*m1 + #X^2*K1*m2 + #X^2*K1*m1 + #X^2*B1*B2 + #X*B2*K1 + #X*B1*K2 + #X*B1*K1 + K1*K2, #X, 4) - B2*g*m1^2*m2^2*root(#X^4*m1*m2 + #X^3*B2*m1 + #X^3*B1*m2 + #X^2*K2*m1 + #X^2*K1*m2 + #X^2*K1*m1 + #X^2*B1*B2 + #X*B2*K1 + #X*B1*K2 + #X*B1*K1 + K1*K2, #X, 2)*root(#X^4*m1*m2 + #X^3*B2*m1 + #X^3*B1*m2 + #X^2*K2*m1 + #X^2*K1*m2 + #X^2*K1*m1 + #X^2*B1*B2 + #X*B2*K1 + #X*B1*K2 + #X*B1*K1 + K1*K2, #X, 3) - B2*g*m1^2*m2^2*root(#X^4*m1*m2 + #X^3*B2*m1 + #X^3*B1*m2 + #X^2*K2*m1 + #X^2*K1*m2 + #X^2*K1*m1 + #X^2*B1*B2 + #X*B2*K1 + #X*B1*K2 + #X*B1*K1 + K1*K2, #X, 2)*root(#X^4*m1*m2 + #X^3*B2*m1 + #X^3*B1*m2 + #X^2*K2*m1 + #X^2*K1*m2 + #X^2*K1*m1 + #X^2*B1*B2 + #X*B2*K1 + #X*B1*K2 + #X*B1*K1 + K1*K2, #X, 4) - B2*g*m1^2*m2^2*root(#X^4*m1*m2 + #X^3*B2*m1 + #X^3*B1*m2 + #X^2*K2*m1 + #X^2*K1*m2 + #X^2*K1*m1 + #X^2*B1*B2 + #X*B2*K1 + #X*B1*K2 + #X*B1*K1 + K1*K2, #X, 3)*root(#X^4*m1*m2 + #X^3*B2*m1 + #X^3*B1*m2 + #X^2*K2*m1 + #X^2*K1*m2 + #X^2*K1*m1 + #X^2*B1*B2 + #X*B2*K1 + #X*B1*K2 + #X*B1*K1 + K1*K2, #X, 4)))/(root(#X^4*m1*m2 + #X^3*B2*m1 + #X^3*B1*m2 + #X^2*K2*m1 + #X^2*K1*m2 + #X^2*K1*m1 + #X^2*B1*B2 + #X*B2*K1 + #X*B1*K2 + #X*B1*K1 + K1*K2, #X, 1)*(m1^2*m2^3*root(#X^4*m1*m2 + #X^3*B2*m1 + #X^3*B1*m2 + #X^2*K2*m1 + #X^2*K1*m2 + #X^2*K1*m1 + #X^2*B1*B2 + #X*B2*K1 + #X*B1*K2 + #X*B1*K1 + K1*K2, #X, 1)^3 - root(#X^4*m1*m2 + #X^3*B2*m1 + #X^3*B1*m2 + #X^2*K2*m1 + #X^2*K1*m2 + #X^2*K1*m1 + #X^2*B1*B2 + #X*B2*K1 + #X*B1*K2 + #X*B1*K1 + K1*K2, #X, 1)^2*m1^2*m2^3*root(#X^4*m1*m2 + #X^3*B2*m1 + #X^3*B1*m2 + #X^2*K2*m1 + #X^2*K1*m2 + #X^2*K1*m1 + #X^2*B1*B2 + #X*B2*K1 + #X*B1*K2 + #X*B1*K1 + K1*K2, #X, 2) - root(#X^4*m1*m2 + #X^3*B2*m1 + #X^3*B1*m2 + #X^2*K2*m1 + #X^2*K1*m2 + #X^2*K1*m1 + #X^2*B1*B2 + #X*B2*K1 + #X*B1*K2 + #X*B1*K1 + K1*K2, #X, 1)^2*m1^2*m2^3*root(#X^4*m1*m2 + #X^3*B2*m1 + #X^3*B1*m2 + #X^2*K2*m1 + #X^2*K1*m2 + #X^2*K1*m1 + #X^2*B1*B2 + #X*B2*K1 + #X*B1*K2 + #X*B1*K1 + K1*K2, #X, 3) - root(#X^4*m1*m2 + #X^3*B2*m1 + #X^3*B1*m2 + #X^2*K2*m1 + #X^2*K1*m2 + #X^2*K1*m1 + #X^2*B1*B2 + #X*B2*K1 + #X*B1*K2 + #X*B1*K1 + K1*K2, #X, 1... Output truncated. Text exceeds maximum line length for Command Window display.
... and similarly for the others.
I gave up on trying to simplify it after 30 minutes with no output:
V1 = simplify(S.v1, 'Steps', 500)
V2 = simplify(S.v2, 'Steps', 500)
X1 = simplify(S.x1, 'Steps', 500)
X2 = simplify(S.x2, 'Steps', 500)
so I leave it to you to determine what the ‘X’ parameters are.
It’s likely better to use the odeToVectorField function and then the matlabFunction function on the first output of odeToVectorField than trying to do this symbolically.

  0 Comments

Sign in to comment.

Community Treasure Hunt

Find the treasures in MATLAB Central and discover how the community can help you!

Start Hunting!