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