function A = sparsemodel(n)

% Model matrix for iterative solvers: 
% 5-point finite difference Laplacian on a square (n-by-n grid)
% Sparse format

A = delsq(numgrid('S',n));
