Maple 2020.
restart; PDE := diff(u(x,y), y$2 ) + diff(u(x,y), x$2) = 0; BC:= u(x,0)=0, u(x,100)=100, u(0,y)=0, u(10,y)=0; sol:=pdsolve(PDE,[BC] ,numeric); Error, (in pdsolve/numeric) unable to handle elliptic PDEs
Compare to
restart; PDE := diff(u(x,y), y$2 ) + diff(u(x,y), x$2) = 0; BC:= u(x,0)=0, u(x,100)=100, u(0,y)=0, u(10,y)=0; sol:=pdsolve([PDE,BC]);