求助各位大神,我用fsolve求解了一个四元非线性方程组,方程比较复杂,得到了结果不符合条件。
代码
syms a b c d e g s r t
e=0.2;
g=0.08;
s=1.5;
r=0.1;
t=2;
eq1=0.014696*(1-s*(1+r))*exp(t*(-(1-a-b-c-d-e)-s*a*(1+r)-s*d))+0.985304*exp(t*(-(1-a-b-c-d-e)))+0.020183903*(s-s^2*(1+r)-(1+r))*exp(t*(-s*(1-a-b-c-d-e)-a*(1+r)-s^2*a*(1+r)-s^2*d))+0.001513038*(s-(1+r))*exp(t*(-s*(1-a-b-c-d-e)-c/0.004-a*(1+r)))+0.581751*(s-(1+r))*exp(t*(-s*(1-a-b-c-d-e)-a*(1+r)))+0.026544385*(s^2-s*(1+r)-s^3*(1+r))*exp(t*(-s^2*(1-a-b-d-e)-a*s*(1+r)-a*s^3*(1+r)-d*s^3))+0.000703601*(s^2-s*(1+r))*exp(t*(-s^2*(1-a-b-d-g)-a*s*(1+r)-c*s/0.004))+0.336909*(s^2-s*(1+r))*exp(t*(-s^2*(1-a-b-d-g)-a*s*(1+r)))+0.219802093*(s^3-s^2*(1+r))*exp(t*(-s^3*(1-a-b-d-g)-a*s^2*(1+r)-d*s^3))-0.132650563*s^3*(1+r)*exp(t*(-s^3*