function [M,z]=trimat(a,b,c,n);
% function [M,z] = trimat(a,b,c,n)
%ertellt eine Tridiagonalmatrix der Größe n x n 
% mit einträgen a auf der Hauptdiagonale, b auf der linken und c auf der
% rechten ersten Nebendiagonale. Zusätzlich wird die Determinante dieser
% Matrix berechnet.
%
% Input:  a   Skalar
%         b   Skalar
%         c   Skalar
%         n   natürliche Zahl
% Output: M   Matrix
%         z   Skalar - Determinante von M
M=a*eye(n); % a mal die Einheitsmatrix
for i=1:(n-1)
    M(i+1,i)=b;
    M(i,i+1)=c;
end
z=det(M);