function  p = dhypg(k,n,K,N)
%DHYPGEO  The hypergeometric probability function
%
%        p = dhypg(k,n,K,N)

%       Anders Holtsberg, 18-11-93
%       Copyright (c) Anders Holtsberg

if any(any( n>N | K>N | K<0 ));
   error('Incompatible input arguments');
end

z = k<0 | n-k<0 | k>K | n-k>N-K;
i = find(~z);
if length(k)>1, k = k(i); end
if length(K)>1, K = K(i); end
if length(n)>1, n = n(i); end
if length(N)>1, N = N(i); end
pp = bincoef(k,K) .* bincoef(n-k,N-K) ./ bincoef(n,N);
p = z*0;
p(i) = pp(:);