# ParametricNDSolve many functions and plotting

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

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