added vector of initial conditions,
code loops accross these, solving each IVP individually.
I tried vectorizing this but it was a mess of 3-D matrices.
Added colored phase-plane plot.
%todo flesh-out chart titles, add color legend to phase portrait