0
$\begingroup$

I don't understand why I don't get the plots :

dzeta = 0.1
alpha = 0.1
beta = 0.1
kappa = 0.1
t0 = 0.5
nc = 1/2
chic = 1/2

sol = ParametricNDSolve[{c'[t] == 
    alpha*c[t]^(2/3)*q[t] - beta*(1 - q[t])*c[t], 
   chi'[t] == kappa*c[t]*q[t] - dzeta*chi[t], 
   q'[t] + q[t]/t0 == chi[t]/((chi[t] + chic)), q[0] == 1/2, 
   chi[t /; t <= 0] == 1/2*(Tanh[-t*10]), c[t /; t <= 0] == C1}, {c, 
   chi, q}, {t, 0, 14}, {C1}]


Plot[Evaluate[c[1][t] /. sol], {t, 0, 14}]

And the plot is empty !! I've looked at ParametricNDSolve documentation and I didn't see what could be wrong...

thx

$\endgroup$
  • 2
    $\begingroup$ Why do you use Condition (/;) in your initial condiions? If you change to chi[0] == 0 and c[0] == C1 then it will work. $\endgroup$ – Carl Woll May 23 at 15:02

Your Answer

By clicking “Post Your Answer”, you agree to our terms of service, privacy policy and cookie policy

Browse other questions tagged or ask your own question.