function Psharp = deblur(Pblur,S)  
% MATLAB function for deblurring via solving a sparse linear system
[m,n] = size(Pblur);
v = buildpsfmatrix(m,n,S)\makeimgvector(Pblur);
Psharp = reshape(v,n,m)';