To solve differential equation (x+y)(p+q)2+(x−y)(p−q)2=1
Let us put x+y=X2 and x−y=Y2
then
p=∂x∂z=∂X∂z∂x∂X+∂Y∂z∂x∂Y=2X1∂X∂z+2Y1∂Y∂z
q=∂y∂z=∂X∂z∂y∂X+∂Y∂z∂y∂Y=2X1∂X∂z−2Y1∂Y∂z
Then
p+q=X1∂X∂z
and
p−q=Y1∂Y∂z
Now substitute it in the equation, then we get
(∂X∂z)2+(∂Y∂z)2=1
Complete integral of the equation is
z=aX+bY+c where a2+b2=1
Put value of X and Y, we get
z=a(x+y)+1−a2(x−y)+c
where a and c are arbitrary constants.