dx1/dt = x2;
dx2/dt = x3;
dx3/dt = x4;
dx4/dt = x5;
dx5/dt = (-siga1*x1^2*beta,SI*phi*M*Mar*Ks^2*siga2^2 + siga1*x1^2*beta,SI*M*Mar*Ks^2*siga2^2 - siga1*x1^2*phi*M^2*Mar*Ks^2*siga2*beta,SA + siga1*x1^2*phi*M^2*Ks^2*siga2^2*beta,SA + siga1*x1^2*M^2*Mar*Ks^2*siga2*beta,SA - siga1*x1*x2*beta,SI*phi*M*Mar*Ks^2*siga2 - 2*siga1*x1*x2*beta,SI*phi*M*Mar*Ks*siga2^2 - siga1*x1*x2*beta,SI*phi*M*Ks^2*siga2^2 - siga1*x1*x2*beta,SI*phi*Mar*Ks^2*siga2^2 + siga1*x1*x2*beta,SI*M*Mar*Ks^2*siga2 + 2*siga1*x1*x2*beta,SI*M*Mar*Ks*siga2^2 + siga1*x1*x2*beta,SI*M*Ks^2*siga2^2 + siga1*x1*x2*beta,SI*Mar*Ks^2*siga2^2 - siga1*x1*x2*phi*M^2*Mar*Ks^2*beta,SA - 2*siga1*x1*x2*phi*M^2*Mar*Ks*siga2*beta,SA + siga1*x1*x2*phi*M^2*Ks^2*siga2*beta,SA + 2*siga1*x1*x2*phi*M^2*Ks*siga2^2*beta,SA - 2*siga1*x1*x2*phi*M*Mar*Ks^2*siga2*beta,SA + 2*siga1*x1*x2*phi*M*Ks^2*siga2^2*beta,SA + siga1*x1*x2*M^2*Mar*Ks^2*beta,SA + 2*siga1*x1*x2*M^2*Mar*Ks*siga2*beta,SA + siga1*x1*x2*M^2*Ks^2*siga2*beta,SA + 2*siga1*x1*x2*M*Mar*Ks^2*siga2*beta,SA - siga1*x1*beta,SI*x4*phi*Mar*Ks*siga2 + siga1*x1*beta,SI*x4*Mar*Ks*siga2 - siga1*x1*beta,SI*phi*M*x3*Mar*Ks*siga2 - siga1*x1*beta,SI*phi*M*x3*Ks*siga2^2 - siga1*x1*beta,SI*phi*x3*Mar*Ks^2*siga2 - siga1*x1*beta,SI*phi*x3*Mar*Ks*siga2^2 + siga1*x1*beta,SI*M*x3*Mar*Ks*siga2 + 
siga1*x1*beta,SI*M*x3*Ks*siga2^2 + siga1*x1*beta,SI*x3*Mar*Ks^2*siga2 + siga1*x1*beta,SI*x3*Mar*Ks*siga2^2 - siga1*x1*x4*phi*M*Mar*Ks*beta,SA + siga1*x1*x4*phi*M*Ks*siga2*beta,SA + siga1*x1*x4*M*Mar*Ks*beta,SA + siga1*x1*x4*M*Ks*siga2*beta,SA - siga1*x1*phi*M^2*x3*Mar*Ks*beta,SA + siga1*x1*phi*M^2*x3*Ks*siga2*beta,SA - siga1*x1*phi*M^2*Mar*x6*Ks*siga2^2 - siga1*x1*phi*M^2*Ks*beta*siga2^2 - siga1*x1*phi*M*x3*Mar*Ks^2*beta,SA - 2*siga1*x1*phi*M*x3*Mar*Ks*siga2*beta,SA + siga1*x1*phi*M*x3*Ks^2*siga2*beta,SA + 2*siga1*x1*phi*M*x3*Ks*siga2^2*beta,SA + siga1*x1*M^2*x3*Mar*Ks*beta,SA + siga1*x1*M^2*x3*Ks*siga2*beta,SA + siga1*x1*M^2*Mar*x6*Ks*siga2^2 + siga1*x1*M^2*Ks*beta*siga2^2 + siga1*x1*M*x3*Mar*Ks^2*beta,SA + 2*siga1*x1*M*x3*Mar*Ks*siga2*beta,SA + siga1*x1*M*x3*Ks^2*siga2*beta,SA - siga1*x2^2*beta,SI*phi*M*Mar*Ks*siga2 - siga1*x2^2*beta,SI*phi*M*Mar*siga2^2 - siga1*x2^2*beta,SI*phi*M*Ks^2*siga2 - siga1*x2^2*beta,SI*phi*M*Ks*siga2^2 - siga1*x2^2*beta,SI*phi*Mar*Ks*siga2^2 - siga1*x2^2*beta,SI*phi*Ks^2*siga2^2 + siga1*x2^2*beta,SI*M*Mar*Ks*siga2 + siga1*x2^2*beta,SI*M*Mar*siga2^2 + siga1*x2^2*beta,SI*M*Ks^2*siga2 + siga1*x2^2*beta,SI*M*Ks*siga2^2 + siga1*x2^2*beta,SI*Mar*Ks*siga2^2 + siga1*x2^2*beta,SI*Ks^2*siga2^2 
- siga1*x2^2*phi*M^2*Mar*Ks*beta,SA - siga1*x2^2*phi*M^2*Mar*siga2*beta,SA + 
siga1*x2^2*phi*M^2*Ks*siga2*beta,SA + siga1*x2^2*phi*M^2*siga2^2*beta,SA - siga1*x2^2*phi*M*Mar*Ks^2*beta,SA - 2*siga1*x2^2*phi*M*Mar*Ks*siga2*beta,SA + siga1*x2^2*phi*M*Ks^2*siga2*beta,SA + 2*siga1*x2^2*phi*M*Ks*siga2^2*beta,SA - 
siga1*x2^2*phi*Mar*Ks^2*siga2*beta,SA + siga1*x2^2*phi*Ks^2*siga2^2*beta,SA + siga1*x2^2*M^2*Mar*Ks*beta,SA + siga1*x2^2*M^2*Mar*siga2*beta,SA + siga1*x2^2*M^2*Ks^2*beta,SA + siga1*x2^2*M^2*Ks*siga2*beta,SA + siga1*x2^2*M*Mar*Ks^2*beta,SA + 2*siga1*x2^2*M*Mar*Ks*siga2*beta,SA + siga1*x2^2*M*Ks^2*siga2*beta,SA + siga1*x2^2*Mar*Ks^2*siga2*beta,SA - siga1*x2*beta,SI*x4*phi*Mar*siga2 - 
siga1*x2*beta,SI*x4*phi*Ks*siga2 + siga1*x2*beta,SI*x4*Mar*siga2 + siga1*x2*beta,SI*x4*Ks*siga2 - siga1*x2*beta,SI*phi*M*x3*Mar*siga2 - 2*siga1*x2*beta,SI*phi*M*x3*Ks*siga2 - siga1*x2*beta,SI*phi*M*x3*siga2^2 - siga1*x2*beta,SI*phi*x3*Mar*Ks*siga2 - siga1*x2*beta,SI*phi*x3*Mar*siga2^2 - siga1*x2*beta,SI*phi*x3*Ks^2*siga2 - 2*siga1*x2*beta,SI*phi*x3*Ks*siga2^2 + siga1*x2*beta,SI*M*x3*Mar*siga2 + 2*siga1*x2*beta,SI*M*x3*Ks*siga2 + siga1*x2*beta,SI*M*x3*siga2^2 + siga1*x2*beta,SI*x3*Mar*Ks*siga2 + siga1*x2*beta,SI*x3*Mar*siga2^2 + siga1*x2*beta,SI*x3*Ks^2*siga2 + 2*siga1*x2*beta,SI*x3*Ks*siga2^2 - siga1*x2*x4*phi*M*Mar*beta,SA + siga1*x2*x4*phi*M*siga2*beta,SA - siga1*x2*x4*phi*Mar*Ks*beta,SA + siga1*x2*x4*phi*Ks*siga2*beta,SA + siga1*x2*x4*M*Mar*beta,SA + 2*siga1*x2*x4*M*Ks*beta,SA + siga1*x2*x4*M*siga2*beta,SA + siga1*x2*x4*Mar*Ks*beta,SA + siga1*x2*x4*Ks*siga2*beta,SA - siga1*x2*phi*M^2*x3*Mar*beta,SA + siga1*x2*phi*M^2*x3*siga2*beta,SA - siga1*x2*phi*M^2*Mar*x6*Ks*siga2 - siga1*x2*phi*M^2*Mar*x6*siga2^2 - siga1*x2*phi*M^2*Ks*beta*siga2 - siga1*x2*phi*M^2*Ks*siga2^2 - siga1*x2*phi*M^2*beta*siga2^2 - 3*siga1*x2*phi*M*x3*Mar*Ks*beta,SA - 2*siga1*x2*phi*M*x3*Mar*siga2*beta,SA + 3*siga1*x2*phi*M*x3*Ks*siga2*beta,SA + 2*siga1*x2*phi*M*x3*siga2^2*beta,SA - siga1*x2*phi*M*Mar*x6*Ks*siga2^2 - siga1*x2*phi*M*Ks*beta*siga2^2 - siga1*x2*phi*x3*Mar*Ks^2*beta,SA - 2*siga1*x2*phi*x3*Mar*Ks*siga2*beta,SA + siga1*x2*phi*x3*Ks^2*siga2*beta,SA + 2*siga1*x2*phi*x3*Ks*siga2^2*beta,SA + siga1*x2*M^2*x3*Mar*beta,SA + 2*siga1*x2*M^2*x3*Ks*beta,SA + siga1*x2*M^2*x3*siga2*beta,SA + siga1*x2*M^2*Mar*x6*Ks*siga2 + siga1*x2*M^2*Mar*x6*siga2^2 + siga1*x2*M^2*Ks*beta*siga2 + siga1*x2*M^2*Ks*siga2^2 + siga1*x2*M^2*beta*siga2^2 + 3*siga1*x2*M*x3*Mar*Ks*beta,SA + 2*siga1*x2*M*x3*Mar*siga2*beta,SA + 2*siga1*x2*M*x3*Ks^2*beta,SA + 3*siga1*x2*M*x3*Ks*siga2*beta,SA + siga1*x2*M*Mar*x6*Ks*siga2^2 + siga1*x2*M*Ks*beta*siga2^2 + siga1*x2*x3*Mar*Ks^2*beta,SA + 2*siga1*x2*x3*Mar*Ks*siga2*beta,SA + siga1*x2*x3*Ks^2*siga2*beta,SA - siga1*beta,SI*x4*phi*x3*siga2 + siga1*beta,SI*x4*x3*siga2 - siga1*beta,SI*phi*M*x3^2*siga2 - siga1*beta,SI*phi*x3^2*Ks*siga2 - siga1*beta,SI*phi*x3^2*siga2^2 + siga1*beta,SI*M*x3^2*siga2 + siga1*beta,SI*x3^2*Ks*siga2 + siga1*beta,SI*x3^2*siga2^2 + siga1*x4^2*beta,SA - siga1*x4*phi*M^2*siga2 - siga1*x4*phi*M*Mar*x6*siga2 - siga1*x4*phi*M*Ks*siga2 - siga1*x4*phi*M*beta*siga2 - siga1*x4*phi*M*siga2^2 - siga1*x4*phi*x3*Mar*beta,SA + siga1*x4*phi*x3*siga2*beta,SA + siga1*x4*M^2*siga2 + 2*siga1*x4*M*x3*beta,SA + siga1*x4*M*Mar*x6*siga2 + siga1*x4*M*Ks*siga2 + siga1*x4*M*beta*siga2 + siga1*x4*M*siga2^2 + siga1*x4*x3*Mar*beta,SA + 2*siga1*x4*x3*Ks*beta,SA + siga1*x4*x3*siga2*beta,SA - siga1*phi*M^2*x3*Mar*x6*siga2 - siga1*phi*M^2*x3*Ks*siga2 - siga1*phi*M^2*x3*beta*siga2 - siga1*phi*M^2*x3*siga2^2 - siga1*phi*M*x5*siga2 - siga1*phi*M*x3^2*Mar*beta,SA + siga1*phi*M*x3^2*siga2*beta,SA - siga1*phi*M*x3*Mar*x6*Ks*siga2 - siga1*phi*M*x3*Mar*x6*siga2^2 - siga1*phi*M*x3*Ks*beta*siga2 - siga1*phi*M*x3*Ks*siga2^2 - siga1*phi*M*x3*beta*siga2^2 - siga1*phi*x3^2*Mar*Ks*beta,SA - siga1*phi*x3^2*Mar*siga2*beta,SA + siga1*phi*x3^2*Ks*siga2*beta,SA + siga1*phi*x3^2*siga2^2*beta,SA + siga1*M^2*x3^2*beta,SA + siga1*M^2*x3*Mar*x6*siga2 + siga1*M^2*x3*Ks*siga2 + siga1*M^2*x3*beta*siga2 + siga1*M^2*x3*siga2^2 + siga1*M*x5*siga2 + siga1*M*x3^2*Mar*beta,SA + 2*siga1*M*x3^2*Ks*beta,SA + siga1*M*x3^2*siga2*beta,SA + siga1*M*x3*Mar*x6*Ks*siga2 + siga1*M*x3*Mar*x6*siga2^2 + siga1*M*x3*Ks*beta*siga2 + siga1*M*x3*Ks*siga2^2 + siga1*M*x3*beta*siga2^2 + siga1*x3^2*Mar*Ks*beta,SA + siga1*x3^2*Mar*siga2*beta,SA + siga1*x3^2*Ks^2*beta,SA + siga1*x3^2*Ks*siga2*beta,SA - x1*x2*beta,SI*phi*M*Mar*Ks^2*siga2^2 + x1*x2*beta,SI*M*Mar*Ks^2*siga2^2 - x1*x2*phi*M^2*Mar*Ks^2*siga2*beta,SA + x1*x2*phi*M^2*Ks^2*siga2^2*beta,SA + x1*x2*M^2*Mar*Ks^2*siga2*beta,SA - x1*beta,SI*x4*phi*M*Mar*Ks*siga2 - x1*beta,SI*x4*phi*Mar*Ks^2*siga2 - x1*beta,SI*x4*phi*Mar*Ks*siga2^2 + x1*beta,SI*x4*M*Mar*Ks*siga2 + x1*beta,SI*x4*Mar*Ks^2*siga2 + x1*beta,SI*x4*Mar*Ks*siga2^2 - x1*beta,SI*phi*alpa*Mar*Ks*siga2 - x1*beta,SI*phi*M*x3*Mar*Ks^2*siga2 - x1*beta,SI*phi*M*x3*Mar*Ks*siga2^2 - x1*beta,SI*phi*x5*Mar*Ks*siga2 - x1*beta,SI*phi*x3*Mar*Ks^2*siga2^2 + x1*beta,SI*alpa*Mar*Ks*siga2 + x1*beta,SI*M*x3*Mar*Ks^2*siga2 + x1*beta,SI*M*x3*Mar*Ks*siga2^2 + x1*beta,SI*x5*Mar*Ks*siga2 + x1*beta,SI*x3*Mar*Ks^2*siga2^2 - x1*x4*phi*M^2*Mar*Ks*beta,SA + x1*x4*phi*M^2*Ks*siga2*beta,SA - x1*x4*phi*M*Mar*Ks^2*beta,SA - x1*x4*phi*M*Mar*Ks*siga2*beta,SA + x1*x4*phi*M*Ks^2*siga2*beta,SA + x1*x4*phi*M*Ks*siga2^2*beta,SA + x1*x4*M^2*Mar*Ks*beta,SA + x1*x4*M*Mar*Ks^2*beta,SA + x1*x4*M*Mar*Ks*siga2*beta,SA - x1*phi*alpa*M*Mar*Ks*beta,SA + x1*phi*alpa*M*Ks*siga2*beta,SA - x1*phi*M^2*x3*Mar*Ks^2*beta,SA - x1*phi*M^2*x3*Mar*Ks*siga2*beta,SA + x1*phi*M^2*x3*Ks^2*siga2*beta,SA + 
x1*phi*M^2*x3*Ks*siga2^2*beta,SA - x1*phi*M*x5*Mar*Ks*beta,SA + x1*phi*M*x5*Ks*siga2*beta,SA - x1*phi*M*x3*Mar*Ks^2*siga2*beta,SA + x1*phi*M*x3*Ks^2*siga2^2*beta,SA + x1*alpa*M*Mar*Ks*beta,SA + x1*M^2*x3*Mar*Ks^2*beta,SA + x1*M^2*x3*Mar*Ks*siga2*beta,SA + x1*M*x5*Mar*Ks*beta,SA + x1*M*x3*Mar*Ks^2*siga2*beta,SA - x2^2*beta,SI*phi*M*Mar*Ks*siga2^2 - x2^2*beta,SI*phi*M*Ks^2*siga2^2 + x2^2*beta,SI*M*Mar*Ks*siga2^2 + x2^2*beta,SI*M*Ks^2*siga2^2 - x2^2*phi*M^2*Mar*Ks*siga2*beta,SA + x2^2*phi*M^2*Ks*siga2^2*beta,SA - x2^2*phi*M*Mar*Ks^2*siga2*beta,SA + x2^2*phi*M*Ks^2*siga2^2*beta,SA + x2^2*M^2*Mar*Ks*siga2*beta,SA 
+ x2^2*M^2*Ks^2*siga2*beta,SA + x2^2*M*Mar*Ks^2*siga2*beta,SA - x2*beta,SI*x4*phi*M*Mar*siga2 - x2*beta,SI*x4*phi*M*Ks*siga2 - x2*beta,SI*x4*phi*Mar*Ks*siga2 - x2*beta,SI*x4*phi*Mar*siga2^2 - x2*beta,SI*x4*phi*Ks^2*siga2 - x2*beta,SI*x4*phi*Ks*siga2^2 + x2*beta,SI*x4*M*Mar*siga2 + x2*beta,SI*x4*M*Ks*siga2 + x2*beta,SI*x4*Mar*Ks*siga2 + x2*beta,SI*x4*Mar*siga2^2 + x2*beta,SI*x4*Ks^2*siga2 + x2*beta,SI*x4*Ks*siga2^2 - x2*beta,SI*phi*alpa*Mar*siga2 - x2*beta,SI*phi*alpa*Ks*siga2 - x2*beta,SI*phi*M*x3*Mar*Ks*siga2 - x2*beta,SI*phi*M*x3*Mar*siga2^2 - x2*beta,SI*phi*M*x3*Ks^2*siga2 - 2*x2*beta,SI*phi*M*x3*Ks*siga2^2 - x2*beta,SI*phi*x5*Mar*siga2 - x2*beta,SI*phi*x5*Ks*siga2 - x2*beta,SI*phi*x3*Mar*Ks*siga2^2 - x2*beta,SI*phi*x3*Ks^2*siga2^2 + x2*beta,SI*alpa*Mar*siga2 + x2*beta,SI*alpa*Ks*siga2 + x2*beta,SI*M*x3*Mar*Ks*siga2 + x2*beta,SI*M*x3*Mar*siga2^2 + x2*beta,SI*M*x3*Ks^2*siga2 + 2*x2*beta,SI*M*x3*Ks*siga2^2 + x2*beta,SI*x5*Mar*siga2 + x2*beta,SI*x5*Ks*siga2 + x2*beta,SI*x3*Mar*Ks*siga2^2 + x2*beta,SI*x3*Ks^2*siga2^2 - x2*x4*phi*M^2*Mar*beta,SA + x2*x4*phi*M^2*siga2*beta,SA - 2*x2*x4*phi*M*Mar*Ks*beta,SA - x2*x4*phi*M*Mar*siga2*beta,SA + 
2*x2*x4*phi*M*Ks*siga2*beta,SA + x2*x4*phi*M*siga2^2*beta,SA - x2*x4*phi*Mar*Ks^2*beta,SA - x2*x4*phi*Mar*Ks*siga2*beta,SA + x2*x4*phi*Ks^2*siga2*beta,SA 
+ x2*x4*phi*Ks*siga2^2*beta,SA + x2*x4*M^2*Mar*beta,SA + x2*x4*M^2*Ks*beta,SA + 2*x2*x4*M*Mar*Ks*beta,SA + x2*x4*M*Mar*siga2*beta,SA + x2*x4*M*Ks^2*beta,SA + 2*x2*x4*M*Ks*siga2*beta,SA + x2*x4*Mar*Ks^2*beta,SA + x2*x4*Mar*Ks*siga2*beta,SA - x2*phi*alpa*M*Mar*beta,SA + x2*phi*alpa*M*siga2*beta,SA - x2*phi*alpa*Mar*Ks*beta,SA + x2*phi*alpa*Ks*siga2*beta,SA - x2*phi*M^2*x3*Mar*Ks*beta,SA - x2*phi*M^2*x3*Mar*siga2*beta,SA + x2*phi*M^2*x3*Ks*siga2*beta,SA + x2*phi*M^2*x3*siga2^2*beta,SA - x2*phi*M^2*Mar*x6*Ks*siga2^2 - x2*phi*M^2*Ks*beta*siga2^2 - x2*phi*M*x5*Mar*beta,SA + x2*phi*M*x5*siga2*beta,SA - x2*phi*M*x3*Mar*Ks^2*beta,SA - 3*x2*phi*M*x3*Mar*Ks*siga2*beta,SA + x2*phi*M*x3*Ks^2*siga2*beta,SA + 3*x2*phi*M*x3*Ks*siga2^2*beta,SA - x2*phi*x5*Mar*Ks*beta,SA + x2*phi*x5*Ks*siga2*beta,SA - x2*phi*x3*Mar*Ks^2*siga2*beta,SA + x2*phi*x3*Ks^2*siga2^2*beta,SA + x2*alpa*M*Mar*beta,SA + x2*alpa*M*Ks*beta,SA + x2*alpa*Mar*Ks*beta,SA + x2*M^2*x3*Mar*Ks*beta,SA + x2*M^2*x3*Mar*siga2*beta,SA + x2*M^2*x3*Ks^2*beta,SA + 2*x2*M^2*x3*Ks*siga2*beta,SA + x2*M^2*Mar*x6*Ks*siga2^2 + x2*M^2*Ks*beta*siga2^2 + x2*M*x5*Mar*beta,SA + x2*M*x5*Ks*beta,SA + x2*M*x3*Mar*Ks^2*beta,SA + 3*x2*M*x3*Mar*Ks*siga2*beta,SA + 2*x2*M*x3*Ks^2*siga2*beta,SA 
+ x2*x5*Mar*Ks*beta,SA + x2*x3*Mar*Ks^2*siga2*beta,SA - beta,SI*x4*phi*M*x3*siga2 - beta,SI*x4*phi*x3*Ks*siga2 - beta,SI*x4*phi*x3*siga2^2 + beta,SI*x4*M*x3*siga2 + beta,SI*x4*x3*Ks*siga2 + beta,SI*x4*x3*siga2^2 - beta,SI*phi*alpa*x3*siga2 - beta,SI*phi*M*x3^2*Ks*siga2 - beta,SI*phi*M*x3^2*siga2^2 - beta,SI*phi*x5*x3*siga2 - beta,SI*phi*x3^2*Ks*siga2^2 + beta,SI*alpa*x3*siga2 + beta,SI*M*x3^2*Ks*siga2 + beta,SI*M*x3^2*siga2^2 + beta,SI*x5*x3*siga2 + beta,SI*x3^2*Ks*siga2^2 + x4^2*M*beta,SA + x4^2*Ks*beta,SA + x4^2*siga2*beta,SA - x4*phi*M^2*Mar*x6*siga2 - x4*phi*M^2*Ks*siga2 - x4*phi*M^2*beta*siga2 - x4*phi*M^2*siga2^2 - x4*phi*M*x3*Mar*beta,SA + x4*phi*M*x3*siga2*beta,SA - x4*phi*M*Mar*x6*Ks*siga2 - x4*phi*M*Mar*x6*siga2^2 - x4*phi*M*Ks*beta*siga2 - x4*phi*M*Ks*siga2^2 - x4*phi*M*beta*siga2^2 - x4*phi*x3*Mar*Ks*beta,SA - x4*phi*x3*Mar*siga2*beta,SA + x4*phi*x3*Ks*siga2*beta,SA + x4*phi*x3*siga2^2*beta,SA + x4*alpa*beta,SA + x4*M^2*x3*beta,SA + x4*M^2*Mar*x6*siga2 + x4*M^2*Ks*siga2 + x4*M^2*beta*siga2 + x4*M^2*siga2^2 + x4*M*x3*Mar*beta,SA + 3*x4*M*x3*Ks*beta,SA + 2*x4*M*x3*siga2*beta,SA + x4*M*Mar*x6*Ks*siga2 + x4*M*Mar*x6*siga2^2 + x4*M*Ks*beta*siga2 + x4*M*Ks*siga2^2 + x4*M*beta*siga2^2 + x4*x5*beta,SA + x4*x3*Mar*Ks*beta,SA + x4*x3*Mar*siga2*beta,SA + x4*x3*Ks^2*beta,SA + 2*x4*x3*Ks*siga2*beta,SA - phi*alpa*M*Mar*x6*siga2 - phi*alpa*M*beta*siga2 - phi*alpa*x3*Mar*beta,SA + phi*alpa*x3*siga2*beta,SA - phi*M^2*x5*siga2 - phi*M^2*x3*Mar*x6*Ks*siga2 - phi*M^2*x3*Mar*x6*siga2^2 - phi*M^2*x3*Ks*beta*siga2 - phi*M^2*x3*Ks*siga2^2 - phi*M^2*x3*beta*siga2^2 - phi*M*x5*Mar*x6*siga2 - phi*M*x5*Ks*siga2 - phi*M*x5*beta*siga2 - phi*M*x5*siga2^2 - phi*M*x3^2*Mar*Ks*beta,SA - 
phi*M*x3^2*Mar*siga2*beta,SA + phi*M*x3^2*Ks*siga2*beta,SA + phi*M*x3^2*siga2^2*beta,SA - phi*M*x3*Mar*x6*Ks*siga2^2 - phi*M*x3*Ks*beta*siga2^2 - phi*x5*x3*Mar*beta,SA + phi*x5*x3*siga2*beta,SA - phi*x3^2*Mar*Ks*siga2*beta,SA + phi*x3^2*Ks*siga2^2*beta,SA + alpa*M*x3*beta,SA + alpa*M*Mar*x6*siga2 + alpa*M*beta*siga2 + alpa*x3*Mar*beta,SA + alpa*x3*Ks*beta,SA + M^2*x5*siga2 + M^2*x3^2*Ks*beta,SA + M^2*x3^2*siga2*beta,SA + M^2*x3*Mar*x6*Ks*siga2 + M^2*x3*Mar*x6*siga2^2 + M^2*x3*Ks*beta*siga2 + M^2*x3*Ks*siga2^2 + M^2*x3*beta*siga2^2 + 
M*x5*x3*beta,SA + M*x5*Mar*x6*siga2 + M*x5*Ks*siga2 + M*x5*beta*siga2 + M*x5*siga2^2 + M*x3^2*Mar*Ks*beta,SA + M*x3^2*Mar*siga2*beta,SA + M*x3^2*Ks^2*beta,SA + 2*M*x3^2*Ks*siga2*beta,SA + M*x3*Mar*x6*Ks*siga2^2 + M*x3*Ks*beta*siga2^2 + x5*x3*Mar*beta,SA + x5*x3*Ks*beta,SA + x3^2*Mar*Ks*siga2*beta,SA + x3^2*Ks^2*siga2*beta,SA) / (phi*M*siga2 - M*siga2);
dx6/dt = (-x1*beta,SI*phi*Mar*Ks*siga2 + x1*beta,SI*Mar*Ks*siga2 - x1*phi*M*Mar*Ks*beta,SA + x1*phi*M*Ks*siga2*beta,SA + x1*M*Mar*Ks*beta,SA - x2*beta,SI*phi*Mar*siga2 - x2*beta,SI*phi*Ks*siga2 + x2*beta,SI*Mar*siga2 + x2*beta,SI*Ks*siga2 - x2*phi*M*Mar*beta,SA + x2*phi*M*siga2*beta,SA - x2*phi*Mar*Ks*beta,SA + x2*phi*Ks*siga2*beta,SA + x2*M*Mar*beta,SA + x2*M*Ks*beta,SA + x2*Mar*Ks*beta,SA - beta,SI*phi*x3*siga2 + beta,SI*x3*siga2 + x4*beta,SA - phi*M*Mar*x6*siga2 - phi*M*beta*siga2 - phi*x3*Mar*beta,SA + phi*x3*siga2*beta,SA + M*x3*beta,SA + M*Mar*x6*siga2 + M*beta*siga2 + x3*Mar*beta,SA + x3*Ks*beta,SA) / 
(phi*M*siga2 - M*siga2);
y1 = x1
