function UpperTailP=binomutp(n,k,p);
% format UpperTailp=binomutp(n,k,p);
% Input n=number of cases
%       k=number of successes
%       p=expected proportion
% Output: UpperTailP Upper-Tail probabilities.
% Reference: Hollander & Wolfe p. 567.
% written by Eugene.Gallagher@umb.edu
if k>n
    error('k must be <= n');
end
i=k:n;
UpperTailP=sum(binopdf(i,n,p));