y(1)' = y(2);
y(2)' = [-y(1)-y(2)^2]/(y(1)+C0);
y(3)' = y(3)/y(4)*y(2)*C1;
y(4)' = y(2)*FF;
C0, C1,C2和C3是 Const.,FF是C2*y(3)^0.5if y(3) > 100; FF是C3*y(3)^0.5如果y(3)<=100。
这种情况下怎么设置y(4)'？
我试过y(4)' = y(2)*((y(3) > 100)*C2*y(3)^0.5 + (y(3)<=100)*C3*y(3)^0.5);了，没用。有什么建议么？
先感谢您。
上面的方程被简化了。
好吧，在使用 ode15s 和其他求解器进行了几次测试后，我意识到矩阵接近奇异矩阵。所以身体状况是个问题。