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