% Diskretisoidaan ja luodaan ohjeen mukaisesti diagonaalit 

n=3;
B=4*diag(ones(1,n))-diag(ones(1,n-1),1)-diag(ones(1,n-1),-1);
O=zeros(n,n);
mI=-eye(n,n);

A=[B mI O;mI B O;O mI B];

A

% Kirjoitin tämän nyt kuitenkin m-fileksi, mutta eipä saiitä
% kummoinen tullut.  Matlabin helpeistä on toki apua komentoja
% tutkittaessa...
