For example, in Maple:

> PDE := diff(u(x,t),t)=1/10*diff(u(x,t),x,x);

> IBC := {u(x,0)=1, u(0,t)=0, D[1] (u)(1,t)=0};

> pds := pdsolve(PDE,IBC,numeric);