results is NaN help

2 vues (au cours des 30 derniers jours)
Rajeev Ray
Rajeev Ray le 7 Déc 2017
Commenté : Rajeev Ray le 8 Déc 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

Connectez-vous pour répondre à cette question.

Réponses (1)

Torsten
Torsten le 7 Déc 2017
You divide by t^0.5. So don't start with t=0, but with t=a small value.
Best wishes
Torsten.
  1 commentaire
Rajeev Ray
Rajeev Ray le 8 Déc 2017
Thank You very much..

Connectez-vous pour commenter.

Catégories

En savoir plus sur Numerical Integration and Differential Equations dans Help Center et File Exchange

Tags

Community Treasure Hunt

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

Start Hunting!

Translated by