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

Hello,

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...