function fxy=dydx3rdorder(x,y)
fxy1 = y(2);
fxy2 = y(3);
fxy3 = -3*y(3) - 5*y(2) -7*y(1) +6;
fxy = [fxy1; fxy2; fxy3];