dy/dx=1+x+y^2+xy^2