dsolve(diff(diff(y(x),x),x)+a^2*sin(y(x))-b*f(x)=0,y(x))