function A = model(n)

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

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