function yprime=f(t,y);
   yprime=[1 1;4 -1]*y;
