function A=lampomatr(n,r)
% 25.11.03
% Muodostetaan lämpöyhtälön differenssimenetelmämatriisi,
%jolla iteroimalla saadaan sarake kerrallaan
% uusia lämpösarakkeita.
%   Syötteet:
% r = menetelmän parametri, r=c^2*k/h^2.
%     (Stabiilisuusehto: r <= 1/2 .)
% n = osavälien lukumäärä. (pisteitä n+1, sisäpisteitä n-1)
%    Tulos:
% Palautetaan (n-1)x(n-1)-matriisi
%
D=diag((1-2*r)*ones(1,n-1));
Yla=diag(r*ones(1,n-2),1);
Ala=diag(r*ones(1,n-2),-1);
A=Ala+D+Yla;
