@Markiyan Hirnyk So in conclusion there is no way to save memory on this since interpolate has been removed for IVP's?
btw this was related to an earlier post http://www.mapleprimes.com/questions/149602-Solving-First-Order-Linear-Coupled-DE-Matrix-Form where I was trying to solve a set of first order, linear ODE's in matrix form, I would love to be able to solve all 961 equations for more accuracy, and the minimun is aroudn 132, I can od about 170 before I get a memory error.
diff~(X,t)=mat.X; #Your system
t0:=1: #Initial t
res:=dsolve(eqs union inits,numeric,parameters=[seq(x0[i],i=1..(L+1)^2)],range=1..15);