Simulation of the Heat equation in the half-space ℝ_{+} with Robin-type boundary conditions; the problem is the following:
starting from a step-like datum. We used finite difference method: a Crank-Nicolson scheme for the time derivative and a centered-three-points scheme for the Laplacian. (Scilab code here)