xy^2y' = x+1