function y=zkc(x,k)
% function y=zkc(x,k)
% zobecnene kombinacni cislo
% y=x*(x-1)*...*(x-k+1)/k!

y=[];
for kk=k
 if kk==0
  y1=1;
 else
  y1=prod((x-(0:kk-1))./(1:kk));
 end
 y=[y,y1];
end
