2.1.59 \(u_x+x u_y=y\) with \(u(0,y)=y^2\) Problem 3.9(a) Lokenath Debnath

problem number 59

Added June 3, 2019.

Problem 3.9(a) nonlinear pde’s by Lokenath Debnath, 3rd edition.

Solve for \(u(x,y)\) \[ u_x+x u_y=y \] With \(u(0,y)=y^2\).

Mathematica

ClearAll["Global`*"]; 
pde =  D[u[x, y], x] +x*D[u[x, y], y]== y; 
ic  = u[0,y]==y^2; 
sol =  AbsoluteTiming[TimeConstrained[DSolve[{pde,ic} ,u[x, y], {x, y}], 60*10]];
 

\[\left \{\left \{u(x,y)\to \frac {x^4}{4}-\frac {x^3}{3}-x^2 y+x y+y^2\right \}\right \}\]

Maple

restart; 
pde := diff(u(x,y),x) + x*diff(u(x,y),y)= y; 
ic  := u(0,y)=y^2; 
cpu_time := timelimit(60*10,CodeTools[Usage](assign('sol',pdsolve([pde,ic],u(x,y))),output='realtime'));
 

\[u \left ( x,y \right ) =-1/3\,{x}^{3}+yx+1/4\,{x}^{4}-{x}^{2}y+{y}^{2}\]

____________________________________________________________________________________