dsolve(diff(y(x),x) = -(-1/x-_F1(y(x)+1/x))/x,y(x))