Added June 3, 2019.
Problem 3.14(g) nonlinear pde’s by Lokenath Debnath, 3rd edition.
Solve for \(u(x,y)\) \[ x u_x+y u_y=x y \] With \(u=\frac {x^2}{2}\) at \(y=x\).
Mathematica ✓
ClearAll["Global`*"]; pde = x*D[u[x, y], x] +y*D[u[x, y], y]== x*y; ic=u[x,x]==x^2/2; sol = AbsoluteTiming[TimeConstrained[DSolve[pde ,u[x, y], {x, y}], 60*10]];
\[\left \{\left \{u(x,y)\to c_1\left (\frac {y}{x}\right )+\frac {x y}{2}\right \}\right \}\] Mathematica does not support this Cauchy data I.C.
Maple ✓
restart; pde :=x*diff(u(x,y),x) + y*diff(u(x,y),y)= x*y; ic := u(x,x)=x^2/2; cpu_time := timelimit(60*10,CodeTools[Usage](assign('sol',pdsolve(pde,u(x,y)) ),output='realtime'));
\[u \left ( x,y \right ) =1/2\,yx+{\it \_F1} \left ( {\frac {y}{x}} \right ) \] Maple does not support this Cauchy data I.C.
____________________________________________________________________________________