Question

matlab大神来帮我看下这个求解微分方程组的程序是怎么回事啊？总是计算到一个时间点就报错，不收敛了

具体的微分方程组如上图所示，其中I是拟合后的关于时间t的函数，具体的m文件程序如下所示：
function dy=wffcz1(t,y)
a1=1/15;A1=160445.089209;a2=1/19;A2=147921.11406183671;
R1=94.25;C=6.897e-11;L1=11.75e-6;
%options=odeset('wffcz1',1e-13);dy=zeros(2,1);dy(1)=(25079999999999999898935095034332140556544638976*t^9 - 3116000000000000083537608917399754195861504*t^8 + 164799999999999990698756673337884147712*t^7 - 4823999999999999755901658355204096*t^6 + 83320000000000003537039785984*t^5 - 770099999999999994757120*t^4 + 3330000000000000000*t^3 - 94230000000000*t^2 + 1885000000*t + 1847/50-((y(1)-R1*((y(1)/A1)^(1/a1)-y(2)))/A2)^(1/a2)-(y(1)/A1)^(1/a1))/C;dy(2)=((((y(1)-R1*((y(1)/A1)^(1/a1)-y(2)))/A2)^(1/a2)-y(2))*R1)/L1;

在matlab工作空间输入调用计算语句如下：
options=odeset('reltol',1e-12,'abstol',1e-12);
[t,y]=ode15s('wffcz1',[0:1e-10:0.000026],[eps,eps],options);

但在计算过程中总是出现如下的提示：
> [t,y]=ode15s('wffcz1',[0:1e-10:0.000026],[eps,eps],options);Warning: Failure at t=2.943201e-006. Unable to meet integration tolerances without reducingthe step size below the smallest value allowed (6.776264e-021) at time t. > In ode15s at 819
哪位帮我看下是怎么回事啊，谢谢了啊，问题解决后另当感谢！

Answer 1

一般这种情况是由于方程发散造成的，检查方程、初值和t的范围吧

追问

初值肯定都是0，检查方程怎么个检查法啊？这个我不是特别懂