Taken from Maple pdsolve help pages, probem 5
Nonlinear first order PDE
Solve for \(f(x,y,z)\)
Mathematica ✓
ClearAll["Global`*"]; pde = D[f[x, y, z], x] + D[f[x, y, z], y]^2 == f[x, y, z] + z; sol = AbsoluteTiming[TimeConstrained[DSolve[pde, f[x, y, z], {x, y, z}], 60*10]];
Maple ✓
restart; pde := diff(f(x,y,z),x) + (diff(f(x,y,z),y))^2 = f(x,y,z)+z; cpu_time := timelimit(60*10,CodeTools[Usage](assign('sol',pdsolve(pde,f(x,y,z),'build')),output='realtime'));
___________________________________________________________________________________________________________________________________________________________________________________________________________________________________________________________________________________________________________________________________________________________________________________________