The following is the code for the 2nd Fick's law with initial/boundary conditions. Unfortunately this code does not show me the result required for next process. It would be useful for us to proceed this work. Thanks in advance. ( Currently we are using Maple 10.) > restart; > DiffusionCoefficientST := 0.5e-2; > ExperienceConstant := 5; > Temperature := 273+25; > PDE := diff(C(x, t), t) = > DiffusionCoefficientST*(exp(1))(-ExperienceConstant(1/Temperature-1/296))*(­diff(C(x, > t), x, x)); > IBCondition := {C(x, 0) = 2, C(0, t) = 0, ((D[1])(C))(5, t) = 0}; > pds := pdsolve(PDE, IBCondition, numeric);

Please Wait...