设dz=(xy(x+y)-f(x)y)dx+(f'(x)+x^2y)dy
∂z/∂y=f'(x)+x^2y
z=f'(x)y+x^2y^2/2+g(x)
∂z/∂x=f''(x)y+xy^2+g'(x)
由:f''(x)y+xy^2+g'(x)=xy(x+y)-f(x)y
f''(x)y+g'(x)=x^2y-f(x)y (要解出f(x),除非g'(x)=0)
f''(x)+f(x)=x^2
解得:f(x)=C1sinx+C2cosx+x^2-2
f(0)=0得:C2=2
f'(x)=C1cosx-C2sinx+x^2-2
f'(0)=1得:C1=3
f(x)=3sinx+2cosx+x^2-2