ExplicitLMM2:=proc(f,y0,y1,dt,n)
local y,Y,i;  
  y[0]:=y0; y[1]:=y1;
  Y:=[[0,y[0]],[dt,y[1]]];
  for i from 1 to n-1 do
    y[i+1]:=-4*y[i]+5*y[i-1]+dt*(4*f(y[i])+2*f(y[i-1]));
    Y:=[op(Y),[dt*(i+1),y[i+1]]];
  od;
  return Y;
end;