function A=lapm(m)
% Laplace (Poisson) "5 point. star discr."
%
%        m -- Jakopisteiden lkm. yhdessä suunnassa (reunat mukaanlukien)
%  Tulos A -- (m-2)^2 x (m-2)^2-matriisi sparse-talletusmuodossa 
n=m-2;  % Sisäsolmut yhdessä suunnassa 
N=n^2; % Sisäsolmujen lkm = tuntemattomien (u(i,j) lkm)
e=ones(N,1);
d=~rem(1:N,n)==0;du=[1,d];du=du(1:N);
A=-spdiags([e,d',-4*e,du',e],[-n,-1,0,1,n],N,N);
