result = true;

for n = [3:100]
  ## Check the sparse matrix against the full one.
  expected = advection_matrix(n,0,1);
  actual = advection_matrix_sparse(n);

  if (expected != actual)
     result = false;
  end
end


unit_test_equals("Sparse advection matrices agree with full ones", ...
		 true, ...
		 result);
