results is NaN help

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

[Image Analyst]

http://www.mathworks.com/matlabcentral/answers/13205#answer_18099

Answer 1

[Torsten]

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

Best wishes

Torsten.

[Rajeev Ray]

Thank You very much..