%Pascal/MATLAB-Kurs WS 2001, Gerhard Opfer
%Der fuer die Darstellung von Flaechen ueber dem R^2 konzipierte Befehl
%surf (von surface=Oberflaeche) kann gut zur Visualisierung wesentlicher 
%Merkmale von Matrizen verwendet werden, insbesondere kann Bandstruktur 
%gut erkannt werden.

n=12;
J=(eye(n)); %eye(n)=(nxn)-Einheitsmatrix
figure(1);
surf(J); view(-45,30); %Blickwinkel: Azimut und Hoehe
zlabel('Einheitsmatrix')
figure(2)
I=[zeros(n-1,1),eye(n-1);zeros(1,n)]; %Superdiagonale ist 1
K=[[zeros(1,n-1);eye(n-1)],zeros(n,1)]; %Subdiagonale ist 1 
surf(0.8*I+J-2*K); view(-45,30);   %Blickwinkel
xlabel('3-Bandmatrix')
figure(3) 
H=invhilb(round(n/2)); %inverse Hilbertmatrix, Eintraege ganzzahlig,
                       %Vorzeichen schachbrettartig,
surf(H); view(-37.5,0);
xlabel('Inverse Hilbertmatrix, frontal');
