function b = bin(n,k);
% function b = bin(n,k);   Binomialkoeffizient

if (n == k | k == 0)
  b = 1;
elseif (0 < k & k < n)
  b = fac(n)/fac(n-k)/fac(k);
else
  error(['erster Eingabewert muss groesser',...
 	 'als der zweite sein']);
end

