From Mathematica DSolve helps pages.
Boundary value problem for the Tricomi equation.
Solve for \(u(x,y)\)
With boundary conditions
Mathematica ✓
ClearAll["Global`*"]; pde = D[u[x, y], {x, 2}] + y*D[u[x, y], {y, 2}] == 0; bc = {u[x, 0] == 0, Derivative[0, 1][u][x, 0] == x^2}; sol = AbsoluteTiming[TimeConstrained[DSolve[{pde, bc}, u[x, y], {x, y}], 60*10]];
Maple ✓
restart; pde := diff(u(x,y),x$2)+ y*diff(u(x,y),y$2)=0; bc := u(x,0)=0,(D[2](u))(x,0)=x^2; cpu_time := timelimit(60*10,CodeTools[Usage](assign('sol',pdsolve([pde,bc],u(x,y))),output='realtime'));
___________________________________________________________________________________________________________________________________________________________________________________________________________________________________________________________________________________________________________________________________________________________________________________________