matlab二阶微分方程组解functionDy=dfun(t,y)D2y(1,1)=(0.2023*y(2)-960-23.85/(0.2023-0.2023*y(1)).^1.4-0.1*Dy(1,1))/96;Dy(2,1)=(-0.2023*Dy(1,1)-0.1*y(2)+2)/(2.4e7*y(1));t=linspace(0,0.15,100);ode113(#39;dfun#39;

问题：matlab二阶微分方程组解functionDy=dfun(t,y)D2y(1,1)=(0.2023*y(2)-960-23.85/(0.2023-0.2023*y(1)).^1.4-0.1*Dy(1,1))/96;Dy(2,1)=(-0.2023*Dy(1,1)-0.1*y(2)+2)/(2.4e7*y(1));t=linspace(0,0.15,100);ode113(#39;dfun#39;,t,[0,16e6]);这有什么问题

答案：↓↓↓

孙作玉的回答：

网友采纳　　把原方程写出来看看,你这种写法比较非主流

孙作玉的回答：

网友采纳　　你的初值有问题，再检查检查　　clc　　f=@(t,x)([x(2);-.0027/96*x(3)+10-9.5e-6/96/(.0088-.0072*x(1))^1.4-.1/96*x(2);(2-.0176*x(2)-.1*x(3))/x(1)/2.4e7]);　　[t,x]=ode45(f,[0.15],[1116e6]);　　subplot(1,2,1),plot(t,x(:,3))　　xlabel('t'),ylabel('p1')　　subplot(1,2,2),plot(t,x(:,1))　　xlabel('t'),ylabel('x')