Added Oct 6, 2019.
Taken from https://www.mapleprimes.com/posts/211274-Integral-Transforms-revamped-And-PDE
Solve \[ {x}^{2}{\frac {\partial ^{2}}{\partial {x}^{2}}}u \left ( x,y \right ) + x{\frac {\partial }{\partial x}}u \left ( x,y \right ) +{\frac { \partial ^{2}}{\partial {y}^{2}}}u \left ( x,y \right ) =0 \] With boundary conditions \begin {align*} u \left ( x,1 \right ) &=\begin {cases} 1& 0\leq x \hspace {5 pt} \text {and} \hspace {5pt} x<1\\ 0&1<x\\ \end {cases}\\ u(x, 0) &= 0 \end {align*}
Mathematica ✗
ClearAll["Global`*"]; pde = x^2*D[u[x, y], {x, 2}] + x*D[u[x, y], x] + D[u[x, y], {y, 2}] ==0; bc = {u[x, 0] == 0, u[x, 1] == Piecewise[{{1, 0 <= x < 1}, {0, 1 < x}}]}; sol = AbsoluteTiming[TimeConstrained[DSolve[{pde, bc}, u[x, y], {x, y}], 60*10]];
Failed
Maple ✓
restart; pde := x^2*diff(u(x, y), x, x) + x*diff(u(x, y), x) + diff(u(x, y), y, y) = 0; iv := u(x, 0) = 0, u(x, 1) = piecewise(0 <= x and x < 1, 1, 1 < x, 0); cpu_time := timelimit(60*10,CodeTools[Usage](assign('sol',pdsolve([pde,iv],u(x,y))),output='realtime'));
\[u \left ( x,y \right ) ={\it invmellin} \left ( {\frac {\sin \left ( sy \right ) }{\sin \left ( s \right ) s}},s,x \right ) \]