If you are attempting to solve ODEs one through five in sequence, one way to do so is GraphicsGrid[{Table[s = NDSolve[eqns[[n]], vars[[n]], {t, 0, 50}]; Plot[Evaluate[vars[[n]] /. s], {t, 0, 50}, PlotRange -> All], {n, 5}]}, ImageSize -> Large] [![enter image description here][1]][1] [1]: https://i.sstatic.net/j6gsu.png