function w=fnlab9p4(x,u);
%% y''=(x-y)/(1+y^2)   y(2)=1, y(2.2)=0.8
w=[u(2);(x-u(1))/(1+(u(1)^2)); u(4);-u(1)/(1+(u(1)^2))]';
