function fxy=dydx2var(x,y)
fxy1 = -0.5*y(1);
fxy2 = 4-0.3*y(2)-0.1*y(1);
fxy = [fxy1; fxy2];