diff(P0(t), t) = P1(t),
diff(P1(t), t) = P2(t),
diff(P2(t), t) = P3(t),
diff(P3(t), t) = P4(t),
diff(P5(t), t) = (-P0(t)*beta_SI*phi*Mar*Ks*siga2 + P0(t)*beta_SI*Mar*Ks*siga2 - P0(t)*phi*M*Mar*Ks*beta_SA + P0(t)*phi*M*Ks*siga2*beta_SA + P0(t)*M*Mar*Ks*beta_SA - P1(t)*beta_SI*phi*Mar*siga2 - P1(t)*beta_SI*phi*Ks*siga2 + P1(t)*beta_SI*Mar*siga2 + P1(t)*beta_SI*Ks*siga2 - P1(t)*phi*M*Mar*beta_SA + P1(t)*phi*M*siga2*beta_SA - P1(t)*phi*Mar*Ks*beta_SA + P1(t)*phi*Ks*siga2*beta_SA + P1(t)*M*Mar*beta_SA + P1(t)*M*Ks*beta_SA + P1(t)*Mar*Ks*beta_SA - beta_SI*phi*P2(t)*siga2 + beta_SI*P2(t)*siga2 + P3(t)*beta_SA - phi*M*Mar*P5(t)*siga2 - phi*M*beta*siga2 - phi*P2(t)*Mar*beta_SA + phi*P2(t)*siga2*beta_SA + M*P2(t)*beta_SA + M*Mar*P5(t)*siga2 + M*beta*siga2 + P2(t)*Mar*beta_SA + P2(t)*Ks*beta_SA) / (phi*M*siga2 - M*siga2),
diff(P4(t), t) = (-siga1*P0(t)^2*beta_SI*phi*M*Mar*Ks^2*siga2^2 + siga1*P0(t)^2*beta_SI*M*Mar*Ks^2*siga2^2 - siga1*P0(t)^2*phi*M^2*Mar*Ks^2*siga2*beta_SA + siga1*P0(t)^2*phi*M^2*Ks^2*siga2^2*beta_SA + siga1*P0(t)^2*M^2*Mar*Ks^2*siga2*beta_SA - siga1*P0(t)*P1(t)*beta_SI*phi*M*Mar*Ks^2*siga2 - 2*siga1*P0(t)*P1(t)*beta_SI*phi*M*Mar*Ks*siga2^2 - siga1*P0(t)*P1(t)*beta_SI*phi*M*Ks^2*siga2^2 - siga1*P0(t)*P1(t)*beta_SI*phi*Mar*Ks^2*siga2^2 + siga1*P0(t)*P1(t)*beta_SI*M*Mar*Ks^2*siga2 + 2*siga1*P0(t)*P1(t)*beta_SI*M*Mar*Ks*siga2^2 + siga1*P0(t)*P1(t)*beta_SI*M*Ks^2*siga2^2 + siga1*P0(t)*P1(t)*beta_SI*Mar*Ks^2*siga2^2 - siga1*P0(t)*P1(t)*phi*M^2*Mar*Ks^2*beta_SA - 2*siga1*P0(t)*P1(t)*phi*M^2*Mar*Ks*siga2*beta_SA + siga1*P0(t)*P1(t)*phi*M^2*Ks^2*siga2*beta_SA + 2*siga1*P0(t)*P1(t)*phi*M^2*Ks*siga2^2*beta_SA - 2*siga1*P0(t)*P1(t)*phi*M*Mar*Ks^2*siga2*beta_SA + 2*siga1*P0(t)*P1(t)*phi*M*Ks^2*siga2^2*beta_SA + siga1*P0(t)*P1(t)*M^2*Mar*Ks^2*beta_SA + 2*siga1*P0(t)*P1(t)*M^2*Mar*Ks*siga2*beta_SA + siga1*P0(t)*P1(t)*M^2*Ks^2*siga2*beta_SA + 2*siga1*P0(t)*P1(t)*M*Mar*Ks^2*siga2*beta_SA - siga1*P0(t)*beta_SI*P3(t)*phi*Mar*Ks*siga2 + siga1*P0(t)*beta_SI*P3(t)*Mar*Ks*siga2 - siga1*P0(t)*beta_SI*phi*M*P2(t)*Mar*Ks*siga2 - siga1*P0(t)*beta_SI*phi*M*P2(t)*Ks*siga2^2 - siga1*P0(t)*beta_SI*phi*P2(t)*Mar*Ks^2*siga2 - siga1*P0(t)*beta_SI*phi*P2(t)*Mar*Ks*siga2^2 + siga1*P0(t)*beta_SI*M*P2(t)*Mar*Ks*siga2 + siga1*P0(t)*beta_SI*M*P2(t)*Ks*siga2^2 + siga1*P0(t)*beta_SI*P2(t)*Mar*Ks^2*siga2 + siga1*P0(t)*beta_SI*P2(t)*Mar*Ks*siga2^2 - siga1*P0(t)*P3(t)*phi*M*Mar*Ks*beta_SA + siga1*P0(t)*P3(t)*phi*M*Ks*siga2*beta_SA + siga1*P0(t)*P3(t)*M*Mar*Ks*beta_SA + siga1*P0(t)*P3(t)*M*Ks*siga2*beta_SA - siga1*P0(t)*phi*M^2*P2(t)*Mar*Ks*beta_SA + siga1*P0(t)*phi*M^2*P2(t)*Ks*siga2*beta_SA - siga1*P0(t)*phi*M^2*Mar*P5(t)*Ks*siga2^2 - siga1*P0(t)*phi*M^2*Ks*beta*siga2^2 - siga1*P0(t)*phi*M*P2(t)*Mar*Ks^2*beta_SA - 2*siga1*P0(t)*phi*M*P2(t)*Mar*Ks*siga2*beta_SA + siga1*P0(t)*phi*M*P2(t)*Ks^2*siga2*beta_SA + 2*siga1*P0(t)*phi*M*P2(t)*Ks*siga2^2*beta_SA + siga1*P0(t)*M^2*P2(t)*Mar*Ks*beta_SA + siga1*P0(t)*M^2*P2(t)*Ks*siga2*beta_SA + siga1*P0(t)*M^2*Mar*P5(t)*Ks*siga2^2 + siga1*P0(t)*M^2*Ks*beta*siga2^2 + siga1*P0(t)*M*P2(t)*Mar*Ks^2*beta_SA + 2*siga1*P0(t)*M*P2(t)*Mar*Ks*siga2*beta_SA + siga1*P0(t)*M*P2(t)*Ks^2*siga2*beta_SA - siga1*P1(t)^2*beta_SI*phi*M*Mar*Ks*siga2 - siga1*P1(t)^2*beta_SI*phi*M*Mar*siga2^2 - siga1*P1(t)^2*beta_SI*phi*M*Ks^2*siga2 - siga1*P1(t)^2*beta_SI*phi*M*Ks*siga2^2 - siga1*P1(t)^2*beta_SI*phi*Mar*Ks*siga2^2 - siga1*P1(t)^2*beta_SI*phi*Ks^2*siga2^2 + siga1*P1(t)^2*beta_SI*M*Mar*Ks*siga2 + siga1*P1(t)^2*beta_SI*M*Mar*siga2^2 + siga1*P1(t)^2*beta_SI*M*Ks^2*siga2 + siga1*P1(t)^2*beta_SI*M*Ks*siga2^2 + siga1*P1(t)^2*beta_SI*Mar*Ks*siga2^2 + siga1*P1(t)^2*beta_SI*Ks^2*siga2^2 - siga1*P1(t)^2*phi*M^2*Mar*Ks*beta_SA - siga1*P1(t)^2*phi*M^2*Mar*siga2*beta_SA + siga1*P1(t)^2*phi*M^2*Ks*siga2*beta_SA + siga1*P1(t)^2*phi*M^2*siga2^2*beta_SA - siga1*P1(t)^2*phi*M*Mar*Ks^2*beta_SA - 2*siga1*P1(t)^2*phi*M*Mar*Ks*siga2*beta_SA + siga1*P1(t)^2*phi*M*Ks^2*siga2*beta_SA + 2*siga1*P1(t)^2*phi*M*Ks*siga2^2*beta_SA - siga1*P1(t)^2*phi*Mar*Ks^2*siga2*beta_SA + siga1*P1(t)^2*phi*Ks^2*siga2^2*beta_SA + siga1*P1(t)^2*M^2*Mar*Ks*beta_SA + siga1*P1(t)^2*M^2*Mar*siga2*beta_SA + siga1*P1(t)^2*M^2*Ks^2*beta_SA + siga1*P1(t)^2*M^2*Ks*siga2*beta_SA + siga1*P1(t)^2*M*Mar*Ks^2*beta_SA + 2*siga1*P1(t)^2*M*Mar*Ks*siga2*beta_SA + siga1*P1(t)^2*M*Ks^2*siga2*beta_SA + siga1*P1(t)^2*Mar*Ks^2*siga2*beta_SA - siga1*P1(t)*beta_SI*P3(t)*phi*Mar*siga2 - siga1*P1(t)*beta_SI*P3(t)*phi*Ks*siga2 + siga1*P1(t)*beta_SI*P3(t)*Mar*siga2 + siga1*P1(t)*beta_SI*P3(t)*Ks*siga2 - siga1*P1(t)*beta_SI*phi*M*P2(t)*Mar*siga2 - 2*siga1*P1(t)*beta_SI*phi*M*P2(t)*Ks*siga2 - siga1*P1(t)*beta_SI*phi*M*P2(t)*siga2^2 - siga1*P1(t)*beta_SI*phi*P2(t)*Mar*Ks*siga2 - siga1*P1(t)*beta_SI*phi*P2(t)*Mar*siga2^2 - siga1*P1(t)*beta_SI*phi*P2(t)*Ks^2*siga2 - 2*siga1*P1(t)*beta_SI*phi*P2(t)*Ks*siga2^2 + siga1*P1(t)*beta_SI*M*P2(t)*Mar*siga2 + 2*siga1*P1(t)*beta_SI*M*P2(t)*Ks*siga2 + siga1*P1(t)*beta_SI*M*P2(t)*siga2^2 + siga1*P1(t)*beta_SI*P2(t)*Mar*Ks*siga2 + siga1*P1(t)*beta_SI*P2(t)*Mar*siga2^2 + siga1*P1(t)*beta_SI*P2(t)*Ks^2*siga2 + 2*siga1*P1(t)*beta_SI*P2(t)*Ks*siga2^2 - siga1*P1(t)*P3(t)*phi*M*Mar*beta_SA + siga1*P1(t)*P3(t)*phi*M*siga2*beta_SA - siga1*P1(t)*P3(t)*phi*Mar*Ks*beta_SA + siga1*P1(t)*P3(t)*phi*Ks*siga2*beta_SA + siga1*P1(t)*P3(t)*M*Mar*beta_SA + 2*siga1*P1(t)*P3(t)*M*Ks*beta_SA + siga1*P1(t)*P3(t)*M*siga2*beta_SA + siga1*P1(t)*P3(t)*Mar*Ks*beta_SA + siga1*P1(t)*P3(t)*Ks*siga2*beta_SA - siga1*P1(t)*phi*M^2*P2(t)*Mar*beta_SA + siga1*P1(t)*phi*M^2*P2(t)*siga2*beta_SA - siga1*P1(t)*phi*M^2*Mar*P5(t)*Ks*siga2 - siga1*P1(t)*phi*M^2*Mar*P5(t)*siga2^2 - siga1*P1(t)*phi*M^2*Ks*beta*siga2 - siga1*P1(t)*phi*M^2*Ks*siga2^2 - siga1*P1(t)*phi*M^2*beta*siga2^2 - 3*siga1*P1(t)*phi*M*P2(t)*Mar*Ks*beta_SA - 2*siga1*P1(t)*phi*M*P2(t)*Mar*siga2*beta_SA + 3*siga1*P1(t)*phi*M*P2(t)*Ks*siga2*beta_SA + 2*siga1*P1(t)*phi*M*P2(t)*siga2^2*beta_SA - siga1*P1(t)*phi*M*Mar*P5(t)*Ks*siga2^2 - siga1*P1(t)*phi*M*Ks*beta*siga2^2 - siga1*P1(t)*phi*P2(t)*Mar*Ks^2*beta_SA - 2*siga1*P1(t)*phi*P2(t)*Mar*Ks*siga2*beta_SA + siga1*P1(t)*phi*P2(t)*Ks^2*siga2*beta_SA + 2*siga1*P1(t)*phi*P2(t)*Ks*siga2^2*beta_SA + siga1*P1(t)*M^2*P2(t)*Mar*beta_SA + 2*siga1*P1(t)*M^2*P2(t)*Ks*beta_SA + siga1*P1(t)*M^2*P2(t)*siga2*beta_SA + siga1*P1(t)*M^2*Mar*P5(t)*Ks*siga2 + siga1*P1(t)*M^2*Mar*P5(t)*siga2^2 + siga1*P1(t)*M^2*Ks*beta*siga2 + siga1*P1(t)*M^2*Ks*siga2^2 + siga1*P1(t)*M^2*beta*siga2^2 + 3*siga1*P1(t)*M*P2(t)*Mar*Ks*beta_SA + 2*siga1*P1(t)*M*P2(t)*Mar*siga2*beta_SA + 2*siga1*P1(t)*M*P2(t)*Ks^2*beta_SA + 3*siga1*P1(t)*M*P2(t)*Ks*siga2*beta_SA + siga1*P1(t)*M*Mar*P5(t)*Ks*siga2^2 + siga1*P1(t)*M*Ks*beta*siga2^2 + siga1*P1(t)*P2(t)*Mar*Ks^2*beta_SA + 2*siga1*P1(t)*P2(t)*Mar*Ks*siga2*beta_SA + siga1*P1(t)*P2(t)*Ks^2*siga2*beta_SA - siga1*beta_SI*P3(t)*phi*P2(t)*siga2 + siga1*beta_SI*P3(t)*P2(t)*siga2 - siga1*beta_SI*phi*M*P2(t)^2*siga2 - siga1*beta_SI*phi*P2(t)^2*Ks*siga2 - siga1*beta_SI*phi*P2(t)^2*siga2^2 + siga1*beta_SI*M*P2(t)^2*siga2 + siga1*beta_SI*P2(t)^2*Ks*siga2 + siga1*beta_SI*P2(t)^2*siga2^2 + siga1*P3(t)^2*beta_SA - siga1*P3(t)*phi*M^2*siga2 - siga1*P3(t)*phi*M*Mar*P5(t)*siga2 - siga1*P3(t)*phi*M*Ks*siga2 - siga1*P3(t)*phi*M*beta*siga2 - siga1*P3(t)*phi*M*siga2^2 - siga1*P3(t)*phi*P2(t)*Mar*beta_SA + siga1*P3(t)*phi*P2(t)*siga2*beta_SA + siga1*P3(t)*M^2*siga2 + 2*siga1*P3(t)*M*P2(t)*beta_SA + siga1*P3(t)*M*Mar*P5(t)*siga2 + siga1*P3(t)*M*Ks*siga2 + siga1*P3(t)*M*beta*siga2 + siga1*P3(t)*M*siga2^2 + siga1*P3(t)*P2(t)*Mar*beta_SA + 2*siga1*P3(t)*P2(t)*Ks*beta_SA + siga1*P3(t)*P2(t)*siga2*beta_SA - siga1*phi*M^2*P2(t)*Mar*P5(t)*siga2 - siga1*phi*M^2*P2(t)*Ks*siga2 - siga1*phi*M^2*P2(t)*beta*siga2 - siga1*phi*M^2*P2(t)*siga2^2 - siga1*phi*M*P4(t)*siga2 - siga1*phi*M*P2(t)^2*Mar*beta_SA + siga1*phi*M*P2(t)^2*siga2*beta_SA - siga1*phi*M*P2(t)*Mar*P5(t)*Ks*siga2 - siga1*phi*M*P2(t)*Mar*P5(t)*siga2^2 - siga1*phi*M*P2(t)*Ks*beta*siga2 - siga1*phi*M*P2(t)*Ks*siga2^2 - siga1*phi*M*P2(t)*beta*siga2^2 - siga1*phi*P2(t)^2*Mar*Ks*beta_SA - siga1*phi*P2(t)^2*Mar*siga2*beta_SA + siga1*phi*P2(t)^2*Ks*siga2*beta_SA + siga1*phi*P2(t)^2*siga2^2*beta_SA + siga1*M^2*P2(t)^2*beta_SA + siga1*M^2*P2(t)*Mar*P5(t)*siga2 + siga1*M^2*P2(t)*Ks*siga2 + siga1*M^2*P2(t)*beta*siga2 + siga1*M^2*P2(t)*siga2^2 + siga1*M*P4(t)*siga2 + siga1*M*P2(t)^2*Mar*beta_SA + 2*siga1*M*P2(t)^2*Ks*beta_SA + siga1*M*P2(t)^2*siga2*beta_SA + siga1*M*P2(t)*Mar*P5(t)*Ks*siga2 + siga1*M*P2(t)*Mar*P5(t)*siga2^2 + siga1*M*P2(t)*Ks*beta*siga2 + siga1*M*P2(t)*Ks*siga2^2 + siga1*M*P2(t)*beta*siga2^2 + siga1*P2(t)^2*Mar*Ks*beta_SA + siga1*P2(t)^2*Mar*siga2*beta_SA + siga1*P2(t)^2*Ks^2*beta_SA + siga1*P2(t)^2*Ks*siga2*beta_SA - P0(t)*P1(t)*beta_SI*phi*M*Mar*Ks^2*siga2^2 + P0(t)*P1(t)*beta_SI*M*Mar*Ks^2*siga2^2 - P0(t)*P1(t)*phi*M^2*Mar*Ks^2*siga2*beta_SA + P0(t)*P1(t)*phi*M^2*Ks^2*siga2^2*beta_SA + P0(t)*P1(t)*M^2*Mar*Ks^2*siga2*beta_SA - P0(t)*beta_SI*P3(t)*phi*M*Mar*Ks*siga2 - P0(t)*beta_SI*P3(t)*phi*Mar*Ks^2*siga2 - P0(t)*beta_SI*P3(t)*phi*Mar*Ks*siga2^2 + P0(t)*beta_SI*P3(t)*M*Mar*Ks*siga2 + P0(t)*beta_SI*P3(t)*Mar*Ks^2*siga2 + P0(t)*beta_SI*P3(t)*Mar*Ks*siga2^2 - P0(t)*beta_SI*phi*alpa*Mar*Ks*siga2 - P0(t)*beta_SI*phi*M*P2(t)*Mar*Ks^2*siga2 - P0(t)*beta_SI*phi*M*P2(t)*Mar*Ks*siga2^2 - P0(t)*beta_SI*phi*P4(t)*Mar*Ks*siga2 - P0(t)*beta_SI*phi*P2(t)*Mar*Ks^2*siga2^2 + P0(t)*beta_SI*alpa*Mar*Ks*siga2 + P0(t)*beta_SI*M*P2(t)*Mar*Ks^2*siga2 + P0(t)*beta_SI*M*P2(t)*Mar*Ks*siga2^2 + P0(t)*beta_SI*P4(t)*Mar*Ks*siga2 + P0(t)*beta_SI*P2(t)*Mar*Ks^2*siga2^2 - P0(t)*P3(t)*phi*M^2*Mar*Ks*beta_SA + P0(t)*P3(t)*phi*M^2*Ks*siga2*beta_SA - P0(t)*P3(t)*phi*M*Mar*Ks^2*beta_SA - P0(t)*P3(t)*phi*M*Mar*Ks*siga2*beta_SA + P0(t)*P3(t)*phi*M*Ks^2*siga2*beta_SA + P0(t)*P3(t)*phi*M*Ks*siga2^2*beta_SA + P0(t)*P3(t)*M^2*Mar*Ks*beta_SA + P0(t)*P3(t)*M*Mar*Ks^2*beta_SA + P0(t)*P3(t)*M*Mar*Ks*siga2*beta_SA - P0(t)*phi*alpa*M*Mar*Ks*beta_SA + P0(t)*phi*alpa*M*Ks*siga2*beta_SA - P0(t)*phi*M^2*P2(t)*Mar*Ks^2*beta_SA - P0(t)*phi*M^2*P2(t)*Mar*Ks*siga2*beta_SA + P0(t)*phi*M^2*P2(t)*Ks^2*siga2*beta_SA + P0(t)*phi*M^2*P2(t)*Ks*siga2^2*beta_SA - P0(t)*phi*M*P4(t)*Mar*Ks*beta_SA + P0(t)*phi*M*P4(t)*Ks*siga2*beta_SA - P0(t)*phi*M*P2(t)*Mar*Ks^2*siga2*beta_SA + P0(t)*phi*M*P2(t)*Ks^2*siga2^2*beta_SA + P0(t)*alpa*M*Mar*Ks*beta_SA + P0(t)*M^2*P2(t)*Mar*Ks^2*beta_SA + P0(t)*M^2*P2(t)*Mar*Ks*siga2*beta_SA + P0(t)*M*P4(t)*Mar*Ks*beta_SA + P0(t)*M*P2(t)*Mar*Ks^2*siga2*beta_SA - P1(t)^2*beta_SI*phi*M*Mar*Ks*siga2^2 - P1(t)^2*beta_SI*phi*M*Ks^2*siga2^2 + P1(t)^2*beta_SI*M*Mar*Ks*siga2^2 + P1(t)^2*beta_SI*M*Ks^2*siga2^2 - P1(t)^2*phi*M^2*Mar*Ks*siga2*beta_SA + P1(t)^2*phi*M^2*Ks*siga2^2*beta_SA - P1(t)^2*phi*M*Mar*Ks^2*siga2*beta_SA + P1(t)^2*phi*M*Ks^2*siga2^2*beta_SA + P1(t)^2*M^2*Mar*Ks*siga2*beta_SA + P1(t)^2*M^2*Ks^2*siga2*beta_SA + P1(t)^2*M*Mar*Ks^2*siga2*beta_SA - P1(t)*beta_SI*P3(t)*phi*M*Mar*siga2 - P1(t)*beta_SI*P3(t)*phi*M*Ks*siga2 - P1(t)*beta_SI*P3(t)*phi*Mar*Ks*siga2 - P1(t)*beta_SI*P3(t)*phi*Mar*siga2^2 - P1(t)*beta_SI*P3(t)*phi*Ks^2*siga2 - P1(t)*beta_SI*P3(t)*phi*Ks*siga2^2 + P1(t)*beta_SI*P3(t)*M*Mar*siga2 + P1(t)*beta_SI*P3(t)*M*Ks*siga2 + P1(t)*beta_SI*P3(t)*Mar*Ks*siga2 + P1(t)*beta_SI*P3(t)*Mar*siga2^2 + P1(t)*beta_SI*P3(t)*Ks^2*siga2 + P1(t)*beta_SI*P3(t)*Ks*siga2^2 - P1(t)*beta_SI*phi*alpa*Mar*siga2 - P1(t)*beta_SI*phi*alpa*Ks*siga2 - P1(t)*beta_SI*phi*M*P2(t)*Mar*Ks*siga2 - P1(t)*beta_SI*phi*M*P2(t)*Mar*siga2^2 - P1(t)*beta_SI*phi*M*P2(t)*Ks^2*siga2 - 2*P1(t)*beta_SI*phi*M*P2(t)*Ks*siga2^2 - P1(t)*beta_SI*phi*P4(t)*Mar*siga2 - P1(t)*beta_SI*phi*P4(t)*Ks*siga2 - P1(t)*beta_SI*phi*P2(t)*Mar*Ks*siga2^2 - P1(t)*beta_SI*phi*P2(t)*Ks^2*siga2^2 + P1(t)*beta_SI*alpa*Mar*siga2 + P1(t)*beta_SI*alpa*Ks*siga2 + P1(t)*beta_SI*M*P2(t)*Mar*Ks*siga2 + P1(t)*beta_SI*M*P2(t)*Mar*siga2^2 + P1(t)*beta_SI*M*P2(t)*Ks^2*siga2 + 2*P1(t)*beta_SI*M*P2(t)*Ks*siga2^2 + P1(t)*beta_SI*P4(t)*Mar*siga2 + P1(t)*beta_SI*P4(t)*Ks*siga2 + P1(t)*beta_SI*P2(t)*Mar*Ks*siga2^2 + P1(t)*beta_SI*P2(t)*Ks^2*siga2^2 - P1(t)*P3(t)*phi*M^2*Mar*beta_SA + P1(t)*P3(t)*phi*M^2*siga2*beta_SA - 2*P1(t)*P3(t)*phi*M*Mar*Ks*beta_SA - P1(t)*P3(t)*phi*M*Mar*siga2*beta_SA + 2*P1(t)*P3(t)*phi*M*Ks*siga2*beta_SA + P1(t)*P3(t)*phi*M*siga2^2*beta_SA - P1(t)*P3(t)*phi*Mar*Ks^2*beta_SA - P1(t)*P3(t)*phi*Mar*Ks*siga2*beta_SA + P1(t)*P3(t)*phi*Ks^2*siga2*beta_SA + P1(t)*P3(t)*phi*Ks*siga2^2*beta_SA + P1(t)*P3(t)*M^2*Mar*beta_SA + P1(t)*P3(t)*M^2*Ks*beta_SA + 2*P1(t)*P3(t)*M*Mar*Ks*beta_SA + P1(t)*P3(t)*M*Mar*siga2*beta_SA + P1(t)*P3(t)*M*Ks^2*beta_SA + 2*P1(t)*P3(t)*M*Ks*siga2*beta_SA + P1(t)*P3(t)*Mar*Ks^2*beta_SA + P1(t)*P3(t)*Mar*Ks*siga2*beta_SA - P1(t)*phi*alpa*M*Mar*beta_SA + P1(t)*phi*alpa*M*siga2*beta_SA - P1(t)*phi*alpa*Mar*Ks*beta_SA + P1(t)*phi*alpa*Ks*siga2*beta_SA - P1(t)*phi*M^2*P2(t)*Mar*Ks*beta_SA - P1(t)*phi*M^2*P2(t)*Mar*siga2*beta_SA + P1(t)*phi*M^2*P2(t)*Ks*siga2*beta_SA + P1(t)*phi*M^2*P2(t)*siga2^2*beta_SA - P1(t)*phi*M^2*Mar*P5(t)*Ks*siga2^2 - P1(t)*phi*M^2*Ks*beta*siga2^2 - P1(t)*phi*M*P4(t)*Mar*beta_SA + P1(t)*phi*M*P4(t)*siga2*beta_SA - P1(t)*phi*M*P2(t)*Mar*Ks^2*beta_SA - 3*P1(t)*phi*M*P2(t)*Mar*Ks*siga2*beta_SA + P1(t)*phi*M*P2(t)*Ks^2*siga2*beta_SA + 3*P1(t)*phi*M*P2(t)*Ks*siga2^2*beta_SA - P1(t)*phi*P4(t)*Mar*Ks*beta_SA + P1(t)*phi*P4(t)*Ks*siga2*beta_SA - P1(t)*phi*P2(t)*Mar*Ks^2*siga2*beta_SA + P1(t)*phi*P2(t)*Ks^2*siga2^2*beta_SA + P1(t)*alpa*M*Mar*beta_SA + P1(t)*alpa*M*Ks*beta_SA + P1(t)*alpa*Mar*Ks*beta_SA + P1(t)*M^2*P2(t)*Mar*Ks*beta_SA + P1(t)*M^2*P2(t)*Mar*siga2*beta_SA + P1(t)*M^2*P2(t)*Ks^2*beta_SA + 2*P1(t)*M^2*P2(t)*Ks*siga2*beta_SA + P1(t)*M^2*Mar*P5(t)*Ks*siga2^2 + P1(t)*M^2*Ks*beta*siga2^2 + P1(t)*M*P4(t)*Mar*beta_SA + P1(t)*M*P4(t)*Ks*beta_SA + P1(t)*M*P2(t)*Mar*Ks^2*beta_SA + 3*P1(t)*M*P2(t)*Mar*Ks*siga2*beta_SA + 2*P1(t)*M*P2(t)*Ks^2*siga2*beta_SA + P1(t)*P4(t)*Mar*Ks*beta_SA + P1(t)*P2(t)*Mar*Ks^2*siga2*beta_SA - beta_SI*P3(t)*phi*M*P2(t)*siga2 - beta_SI*P3(t)*phi*P2(t)*Ks*siga2 - beta_SI*P3(t)*phi*P2(t)*siga2^2 + beta_SI*P3(t)*M*P2(t)*siga2 + beta_SI*P3(t)*P2(t)*Ks*siga2 + beta_SI*P3(t)*P2(t)*siga2^2 - beta_SI*phi*alpa*P2(t)*siga2 - beta_SI*phi*M*P2(t)^2*Ks*siga2 - beta_SI*phi*M*P2(t)^2*siga2^2 - beta_SI*phi*P4(t)*P2(t)*siga2 - beta_SI*phi*P2(t)^2*Ks*siga2^2 + beta_SI*alpa*P2(t)*siga2 + beta_SI*M*P2(t)^2*Ks*siga2 + beta_SI*M*P2(t)^2*siga2^2 + beta_SI*P4(t)*P2(t)*siga2 + beta_SI*P2(t)^2*Ks*siga2^2 + P3(t)^2*M*beta_SA + P3(t)^2*Ks*beta_SA + P3(t)^2*siga2*beta_SA - P3(t)*phi*M^2*Mar*P5(t)*siga2 - P3(t)*phi*M^2*Ks*siga2 - P3(t)*phi*M^2*beta*siga2 - P3(t)*phi*M^2*siga2^2 - P3(t)*phi*M*P2(t)*Mar*beta_SA + P3(t)*phi*M*P2(t)*siga2*beta_SA - P3(t)*phi*M*Mar*P5(t)*Ks*siga2 - P3(t)*phi*M*Mar*P5(t)*siga2^2 - P3(t)*phi*M*Ks*beta*siga2 - P3(t)*phi*M*Ks*siga2^2 - P3(t)*phi*M*beta*siga2^2 - P3(t)*phi*P2(t)*Mar*Ks*beta_SA - P3(t)*phi*P2(t)*Mar*siga2*beta_SA + P3(t)*phi*P2(t)*Ks*siga2*beta_SA + P3(t)*phi*P2(t)*siga2^2*beta_SA + P3(t)*alpa*beta_SA + P3(t)*M^2*P2(t)*beta_SA + P3(t)*M^2*Mar*P5(t)*siga2 + P3(t)*M^2*Ks*siga2 + P3(t)*M^2*beta*siga2 + P3(t)*M^2*siga2^2 + P3(t)*M*P2(t)*Mar*beta_SA + 3*P3(t)*M*P2(t)*Ks*beta_SA + 2*P3(t)*M*P2(t)*siga2*beta_SA + P3(t)*M*Mar*P5(t)*Ks*siga2 + P3(t)*M*Mar*P5(t)*siga2^2 + P3(t)*M*Ks*beta*siga2 + P3(t)*M*Ks*siga2^2 + P3(t)*M*beta*siga2^2 + P3(t)*P4(t)*beta_SA + P3(t)*P2(t)*Mar*Ks*beta_SA + P3(t)*P2(t)*Mar*siga2*beta_SA + P3(t)*P2(t)*Ks^2*beta_SA + 2*P3(t)*P2(t)*Ks*siga2*beta_SA - phi*alpa*M*Mar*P5(t)*siga2 - phi*alpa*M*beta*siga2 - phi*alpa*P2(t)*Mar*beta_SA + phi*alpa*P2(t)*siga2*beta_SA - phi*M^2*P4(t)*siga2 - phi*M^2*P2(t)*Mar*P5(t)*Ks*siga2 - phi*M^2*P2(t)*Mar*P5(t)*siga2^2 - phi*M^2*P2(t)*Ks*beta*siga2 - phi*M^2*P2(t)*Ks*siga2^2 - phi*M^2*P2(t)*beta*siga2^2 - phi*M*P4(t)*Mar*P5(t)*siga2 - phi*M*P4(t)*Ks*siga2 - phi*M*P4(t)*beta*siga2 - phi*M*P4(t)*siga2^2 - phi*M*P2(t)^2*Mar*Ks*beta_SA - phi*M*P2(t)^2*Mar*siga2*beta_SA + phi*M*P2(t)^2*Ks*siga2*beta_SA + phi*M*P2(t)^2*siga2^2*beta_SA - phi*M*P2(t)*Mar*P5(t)*Ks*siga2^2 - phi*M*P2(t)*Ks*beta*siga2^2 - phi*P4(t)*P2(t)*Mar*beta_SA + phi*P4(t)*P2(t)*siga2*beta_SA - phi*P2(t)^2*Mar*Ks*siga2*beta_SA + phi*P2(t)^2*Ks*siga2^2*beta_SA + alpa*M*P2(t)*beta_SA + alpa*M*Mar*P5(t)*siga2 + alpa*M*beta*siga2 + alpa*P2(t)*Mar*beta_SA + alpa*P2(t)*Ks*beta_SA + M^2*P4(t)*siga2 + M^2*P2(t)^2*Ks*beta_SA + M^2*P2(t)^2*siga2*beta_SA + M^2*P2(t)*Mar*P5(t)*Ks*siga2 + M^2*P2(t)*Mar*P5(t)*siga2^2 + M^2*P2(t)*Ks*beta*siga2 + M^2*P2(t)*Ks*siga2^2 + M^2*P2(t)*beta*siga2^2 + M*P4(t)*P2(t)*beta_SA + M*P4(t)*Mar*P5(t)*siga2 + M*P4(t)*Ks*siga2 + M*P4(t)*beta*siga2 + M*P4(t)*siga2^2 + M*P2(t)^2*Mar*Ks*beta_SA + M*P2(t)^2*Mar*siga2*beta_SA + M*P2(t)^2*Ks^2*beta_SA + 2*M*P2(t)^2*Ks*siga2*beta_SA + M*P2(t)*Mar*P5(t)*Ks*siga2^2 + M*P2(t)*Ks*beta*siga2^2 + P4(t)*P2(t)*Mar*beta_SA + P4(t)*P2(t)*Ks*beta_SA + P2(t)^2*Mar*Ks*siga2*beta_SA + P2(t)^2*Ks^2*siga2*beta_SA) / (phi*M*siga2 - M*siga2),
y(t) = P0(t)
