Added June 3, 2019.
Problem 3.8(d) nonlinear pde’s by Lokenath Debnath, 3rd edition.
Solve for \(u(x,y,z)\)
Mathematica ✓
ClearAll["Global`*"]; pde = y*z*D[u[x, y,z], x] - x*z*D[u[x, y,z], y]+x*y*(x^2+y^2)*D[u[x, y,z], z]== 0; sol = AbsoluteTiming[TimeConstrained[DSolve[pde ,u[x, y,z], {x, y,z}], 60*10]];
Maple ✓
restart; pde :=y*z*diff(u(x,y,z),x) - x*z*diff(u(x,y,z),y)+x*y*(x^2+y^2)*diff(u(x,y,z),z)= 0; cpu_time := timelimit(60*10,CodeTools[Usage](assign('sol',pdsolve(pde,u(x,y,z))),output='realtime'));
___________________________________________________________________________________________________________________________________________________________________________________________________________________________________________________________________________________________________________________________________________________________________________________________