diff(diff(y(x),x),x)+a^2*sin(y(x))-b*sin(x)