function dydt=linder(t,y)
A=[-1 2;3 -1];
dydt=A*y;