dt = 0.01;
dx = 0.1;
L = 1;
tmax = 0.04;

function lbVals = alph(t)
  lbVals = zeros(size(t));
end

function rbVals = bet(t)
  rbVals = zeros(size(t));
end

function iVals = ic(x)
  iVals = (x <= 1/5) .* (-x);
  iVals += (x > 1/5 & x <= 7/10) .* (x - 2/5);
  iVals += (x > 7/10) .* (1 - x);
end

soln = fdHeatSolver(dt, dx, L, tmax, @alph, @bet, @ic);
for j = 1:size(soln)(1)
  plot(0:dx:L, soln(j,:))
  axis([0 1 -1 1])
  pause
end