NDSolve and Plot report Replace Issue!

21 Views Asked by At

NDSolve and Plot

Thanks Anyone know what's the issue? NDSolve[{x'[t] == x[t]*(1 - x[t]) - 9*x[t]*y[t]/(10*x[t] + 1), y'[t] == 0.3*y[t]*(1 - 10*y[t]/(10*x[t] + 1)), x[0] == 0.15, y[0] == 0.25}, {x, y}, {t, 0, 150}] Plot[Evaluate[x[t] /. %%], {t, 0, 150}]

1

There are 1 best solutions below

0
On

It is a very long time I did not use Mathematica. Looking at the documentation, writing

 s=NDSolve[{x'[t] == x[t]*(1 - x[t]) - 9*x[t]*y[t]/(10*x[t] + 1),
            y'[t] == 0.3*y[t]*(1 - 10*y[t]/(10*x[t] + 1)), x[0] == 0.15, y[0] == 0.25},
            {x, y}, {t, 0, 150}] 

 Plot[Evaluate[x[t] /. s], {t, 0, 150}, PlotRange -> All]
 Plot[Evaluate[y[t] /. s], {t, 0, 150}, PlotRange -> All]

would (in principle !) not make problems.

Thanks for letting me know.