Hello,

I have found numerous ways of plotting PDEs, but I am trying to ask Maple to calculate the simple analytical solution of one.

Now, Maple is very happy to solve the following with one initial condition:

restart;
with(PDEtools);

pde := diff(u(z, t), t)+c*(diff(u(z, t), z)) = A;
IBC := (u(z, 0) = f(z));

sys:=[pde,IBC];

ans:=pdsolve(sys); 


However, when we take, 
IBC := (u(z, 0) = f(z), D[1](u)(0, t) = 0, D[1](u)(-h, t) = 0);