function A = model(k)

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

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