%script to run diffuse
%This script runs stadard diffusion on a vertical scarp


%You should change these below:
dx = 1;   %usually m

Hinitial = [1 1 1 1 1 1 1 1 1 1 1 -1 -1 -1 -1 -1 -1 -1 -1 -1 -1 -1];

k = 1   %usually m^2/kyr
finaltime = 50;   %usually kyr

pt = 100 ; %every how many time steps you want an output profile drawn

%make sure that Harb and Hinitial have the same length (number of elements or there will be a problem.  Or just comment this line out with a precending %:
Harb = [0 0 0 0 0 0 0 0 0 0 0 0 0 0 0 0 0 0 0 0 0 0];



%don't mess with this:
diffuse(dx, Hinitial, k, finaltime, pt, Harb)