Input isn't a differential equation error

Questions
,
1

I have been working up a pretty simple code and I keep running up against an error saying input isn’t a differential equation. Any help would be much appreciated.

Clear["Global'*"];
ka = .2; kd = .02; km = 150; ns = 10; sens = 3; ro2 = 50; a = .2;
odeca = ka nut[t] cella[t] - 
   kd (1 - (1 + (nut[t]/ns)^(sens))) cella[t];
odecd = kd (1 - (1 + (nut[t]/ns)^(sens))) cella[t];
odeo2 = ro2 (1200 - o2[t]) - km cella[t];
oden = -ka nut[t] cella[t];
vars = {cella[t], celld[t], o2[t], nut[t]};
solution = 
  NDSolve[{cella'[t] == odeca, celld'[t] == odecd, o2'[t] == odeo2, 
    nut'[t] == oden, cella[0] == 10, celld[0] == 1, o2[0] == 1200, 
    n[0] == 200}, vars, {t, 0, 31}];
vars;
vars1 = vars /. solution /. t -> 31; vars2 = Flatten[vars1]
2

If you look at your code, the call to NDSolve you give is:

NDSolve[{
  Derivative[1][cella][t] == 0.2 cella[t] nut[t] + 0.00002 cella[t] nut[t]^3, 
  Derivative[1][celld][t] == -0.00002 cella[t] nut[t]^3, 
  Derivative[1][o2][t] == -150 cella[t] + 50 (1200 - o2[t]), 
  Derivative[1][nut][t] == -0.2 cella[t] nut[t], 
cella[0] == 10, celld[0] == 1, o2[0] == 1200, n[0] == 200}, 
{cella[t], celld[t], o2[t], nut[t]}, {t, 0, 31}]

The first argument to NDSolve is a list of 8 equations. I would take a look at every equation individually and ask myself whether they make sense.
The first four are differential equations. The last four are initial conditions. The very last equation you give is incorrect. Correct the very last equation and NDSolve should evaluate.

3

In[25]:= Clear[“Global’*”];
ka = .2; kd = .02; km = 150; ns = 10; sens = 3; ro2 = 50; a = .2;
solution = NDSolve[{
cella’[t] ==
ka nut[t] cella[t] - kd (1 - (1 + (nut[t]/ns)^(sens))) cella[t],
celld’[t] == kd (1 - (1 + (nut[t]/ns)^(sens))) cella[t],
o2’[t] == ro2 (1200 - o2[t]) - km cella[t],
nut’[t] == -ka nut[t] cella[t],
cella[0] == 10, celld[0] == 1, o2[0] == 1200, nut[0] == 200},
{cella[t], celld[t], o2[t], nut[t]}, {t, 0, 31}]

Out[27]= {{cella[t] -> InterpolatingFunction[{{0., 31.}}, <>][t], 
  celld[t] -> InterpolatingFunction[{{0., 31.}}, <>][t], 
  o2[t] -> InterpolatingFunction[{{0., 31.}}, <>][t], 
  nut[t] -> InterpolatingFunction[{{0., 31.}}, <>][t]}}

I think there were two problems. Your initial condition for nut was give as n[0]. Also your intermediate variables, odeca, etc. were not explicit functions of time and problem variables.

4

Thanks, can’t believe it was only a typo but thanks for loaning me a second pair of eyes to look at it.