Question: pdsolve with elliptic equation

restart;

pde := diff(u(x, y), `$`(x, 2))+4*(diff(u(x, y), `$`(y, 2))) = -x*y;

ibc := {u(0, y) = 0, u(2, y) = 0, u(x, 0) = x*(4-x^2), (D[2](u))(x, 0) = 0};

sol := pdsolve(pde, ibc, numeric);

sol:-plot3d(x = 0 .. 2, y = 0 .. 1, numpoints = 2000, axes = boxed);

 

The final command produce a figure with some sigularity.

how to fix it? I need help, thx you very much.

Please Wait...