Question: error in solving diffusion equation


I have a problem with maple solving diffusion equation. I solved the equation as follows

Es := 0.117108e12;
Ef := 0.78125e11;
l := 0.150e-6;
s := 0.500000e-3;
f := 0.5898334197e-6;
o := 0.9e-5;
d := 0.10e-17;
cb := 0.1e7/(19.9);
c := l*f/(d*cb);
PDE := diff(u(x, t), t)-(diff(u(x, t), x, x)) = 0;
IBC1 := {u(x, 0) = 0, (D[1](u))(0, t) = 0, (D[1](u))(1, t) = c};
S1 := pdsolve(PDE, IBC1, numeric, time = t);
S1:-plot(t = 0.6e-2);

i am plotting cincentration along the thickness of the material. i want to know why the concentration is negative. It should not be negative because constant flux 'c' is present at x=1. and flux at x=0 is 0, and also the initial concentration is 0.



Please Wait...