为什么我这个程序出来y值全是NAN？

function dx= fun4(t,x) %建立函数 g=10; m=0.03; r=sqrt(x(1)^2+x(3)^2); r1=-x(1)/r; r2=-x(3)/r; dx=[x(2);g*m*r1/r^2;x(4);g*m*r2/r^2]; end x0=[0;0;0;0]; %调用函数 [t,y]=ode45(@fun4,[0,20],x0) plot(y(:,1),y(:,3))

初始值为全0的时候， r=sqrt(x(1)^2+x(3)^2)的取值为0了，r1=-x(1)/r; r2=-x(3)/r均为0/0，即为nan，nan与任何数运算返回值均为nan。解决办法是修改初始条件为非0值。

x0=[1;-1;2;-2]; %调用函数 
[t,y]=ode45(@fun4,[0,20],x0) 
plot(y(:,1),y(:,3))

您好，我是有问必答小助手，您的问题已经有小伙伴解答了，您看下是否解决，可以追评进行沟通哦~

如果有您比较满意的答案 / 帮您提供解决思路的答案，可以点击【采纳】按钮，给回答的小伙伴一些鼓励哦～～

ps:问答VIP仅需29元，即可享受5次/月 有问必答服务，了解详情>>>https://vip.csdn.net/askvip?utm_source=1146287632