results is NaN help

3 visualizaciones (últimos 30 días)
Rajeev Ray
Rajeev Ray el 7 de Dic. de 2017
Comentada: Rajeev Ray el 8 de Dic. de 2017
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

Iniciar sesión para responder a esta pregunta.

Respuestas (1)

Torsten
Torsten el 7 de Dic. de 2017
You divide by t^0.5. So don't start with t=0, but with t=a small value.
Best wishes
Torsten.
  1 comentario
Rajeev Ray
Rajeev Ray el 8 de Dic. de 2017
Thank You very much..

Iniciar sesión para comentar.

Categorías

Más información sobre Numerical Integration and Differential Equations en Help Center y File Exchange.

Etiquetas

Community Treasure Hunt

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

Start Hunting!

Translated by