Taken from Mathematica DSolve help pages
Another example of nonlinear Clairaut equation
Solve for \(u(x,y)\)
Mathematica ✓
ClearAll["Global`*"]; pde = u[x, y] == x*D[u[x, y], x] + y*D[u[x, y], y] + Sin[D[u[x, y], x] + D[u[x, y], y]]; sol = AbsoluteTiming[TimeConstrained[DSolve[pde, u[x, y], {x, y}], 60*10]];
Maple ✓
restart; pde := u(x,y)= x*diff(u(x,y),x)+y*diff(u(x,y),y)+sin(diff(u(x,y),x)+diff(u(x,y),y)); cpu_time := timelimit(60*10,CodeTools[Usage](assign('sol',pdsolve(pde,u(x,y))),output='realtime'));
___________________________________________________________________________________________________________________________________________________________________________________________________________________________________________________________________________________________________________________________________________________________________________________________