results is NaN help

2017 12 7
1 个回答
更新时间：2017 12 8
10 次查看（30 天）

0 个投票

function dydt = odefcn(t,y,A,B,C) dydt = zeros(2,1); dydt(1) = -A*y(1)-((B*y(1)^2)/(t^0.5))+C*y(2); dydt(2) = ((0.5*B*y(1)^2)/(t^0.5))-C*y(2); end
A = 0.000001; B = 0.05; C = 0.15; tspan = [0 10]; y0 = [0 0.0001]; [t,y] = ode45(@(t,y) odefcn(t,y,A,B,C), tspan, y0); plot(t,y(:,1),'-o',t,y(:,2),'-.')
It is giving results NaN NaN please help

请先登录，再回答此问题。

回答（1 个）

Torsten
Torsten 2017-12-7

0 个投票

You divide by t^0.5. So don't start with t=0, but with t=a small value.
Best wishes
Torsten.

类别

在 帮助中心 和 File Exchange 中查找有关 Programming 的更多信息

标签

Community Treasure Hunt

Find the treasures in MATLAB Central and discover how the community can help you!

Start Hunting!

Translated by