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);