function[P] = pascal(n)

% function[P] = pascal(n)
% 
% Diese Funktion berechnet das n-te Pascalische Dreieck, womit man 
% zeilenweise die Koeffizienten der Polynome (x+y)^(n-1) erhaelt. 
% Die Rekursionsformel ist
%    pascal(row, col) = pascal(row-1, col-1) + pascal(row-1, col)
% Dabei ist pascal(1,1) = 1, und die fehlenden Eintraege sind als 
% Null zu verstehen.
% Geliefert wird die entsprechende Dreiecksmatrix.

n = fix(n);
if (n < 0) 
    error('Der Eingabeparamter muss eine natuerliche Zahl sein !!');
    return;
end
if (n == 0)
    P = 1;
    return
end

P = zeros(n,n+1);
P(1,2) = 1;
for row = 2:n
    for col = 2:n+1
        P(row, col) = P(row-1, col-1) + P(row-1, col);
    end
end

P(:,1) = [];
