dsolve(diff(y(x),x)+y(x)^2+(x*y(x)-1)*f(x) = 0,y(x))