function fxy=dyydxx(x,y)
dydx1=-0.5*y(1);
dydx2=4-0.3*y(2)-0.1*y(1);
fxy=[dydx1 ; dydx2];