Question: How to solve a system of ODEs first and then plot later?


I am solving a large system of ODEs, using the following command,

> Sol := dsolve({seq(ode[j], j = 0 .. 21), seq(v[j](0) = 0, j = 1 .. 21), v[0](0) = 1}, [seq(v[j](t), j = 0 .. 21)]);



and then plot the quantities I want by something like

> plots:-odeplot(Sol, [t, v[3](t)+v[5](t)], t = 0 .. 1.5);

My problem is that, I do not know a priori which quantity I want to plot, and plotting using above method requires solving the ODEs each time separately, which takes a long time.


So I was curious if there is a scheme that I can solve my system for once and for all, and then plot any quantities that I would like to see.

Please Wait...