function F=imeu(x)
global h
global y0
F=x-h.*[x(2); -3*sin(x(1))-0.2*x(2)]'-y0;