function yprime=linearde(t,y,flag,a,b,c,d)
yprime=[a*y(1)+b*y(2); c*y(1)+d*y(2)];
