function [A,Asp]=lampomatr(n,r)
% 4.12.08
% Muodostetaan lämpöyhtälön differenssimenetelmämatriisi,
%jolla iteroimalla saadaan sarake kerrallaan
% uusia lämpösarakkeita (kuvassa rivejä)
%   Syötteet:
% r = menetelmän parametri, r=c^2*k/h^2.
%     (Stabiilisuusehto: r <= 1/2 .)
% n = Pisteiden lkm reunat mukaanlukien, sisäpisteitä n-2
%    Tulos:
% Palautetaan (n-2)x(n-2)-matriisi
%
D=diag((1-2*r)*ones(1,n-2));
Yla=diag(r*ones(1,n-3),1);
Ala=diag(r*ones(1,n-3),-1);
A=Ala+D+Yla;
e=ones(n-2,1);
Asp=spdiags([r*e (1-2*r)*e r*e],-1:1,n-2,n-2);

