Consider the following toy example:
sol = NDSolve[{
x'[t] == {-1, -2} x[t],
x[0] == {2, 2}
},
x, {t, 0, 4}
];
Plot[x[t] /. sol, {t, 0, 4}, PlotRange -> All]
As seen above, both plots are drawn in the same colour.
The usual trick of adding Evaluate
does not work here, I guess because x[t] /. sol
evaluates to a list with a single element:
and it's the InterpolatingFunction
that later evaluates to a list.
Manually specifying the PlotStyle
also doesn't seem to work.
How can I ensure that the different plots are assigned different colours/styles?