Dear Maple users

I have a question about applying pdsolve MAPLE for solving two dimensional heat equations:

My codes have been provided but it shows to me this error:

Error, (in pdsolve/numeric/process_PDEs) can only numerically solve PDE with two independent variables, got {t, x, y}

If kindly is possible, please help me in this case.

With kind regards,

Emran Tohidi.

> restart;

> with(plots);

print(??); # input placeholder

> with(PDEtools);

print(??); # input placeholder

> declare(u(x, y, t));

print(`output redirected...`); # input placeholder

u(x, y, t) will now be displayed as u

> S := 1/100; tR := 0 .. 1; xR := 0 .. 1; yR := 0 .. 1; NF := 30; NP := 100;

print(??); # input placeholder

> N := 3; L1 := [red, blue, green]; L2 := [0, 1/2, 1]; Ops := spacestep = S, timestep = S;

print(??); # input placeholder

> Op1 := frames = NF, numpoints = NP;

print(??); # input placeholder

> PDE1 := diff(u(x, y, t), t)-(diff(u(x, y, t), `$`(x, 2)))-(diff(u(x, y, t), `$`(y, 2))) = 0;

print(??); # input placeholder

> IC := {u(x, y, 0) = exp(x+y)}; BC := {u(0, y, t) = exp(2*t+y), u(1, y, t) = exp(2*t+y+1), u(x, 0, t) = exp(2*t+x), u(x, 1, t) = exp(2*t+x+1)};

print(??); # input placeholder

> Sol := pdsolve(PDE1, `union`(IC, BC), numeric, u(x, t), Ops);

Error, (in pdsolve/numeric/process_PDEs) can only numerically solve PDE with two independent variables, got {t, x, y}