function x = binoinv(y,n,p)
%BINOINV Inverse of the binomial cumulative distribution function (cdf).
%   X = BINOINV(Y,N,P) returns the inverse of the binomial cdf with 
%   parameters N and P. Since the binomial distribution is
%   discrete, BINOINV returns the least integer X such that 
%   the binomial cdf evaluated at X, equals or exceeds Y.
%
%   The size of X is the common size of the input arguments. A scalar input  
%   functions as a constant matrix of the same size as the other inputs.    
%
%   Note that X takes the values 0,1,2,...,N.
%
%   See also BINOCDF, BINOFIT, BINOPDF, BINORND, BINOSTAT, ICDF.

%   Copyright 1993-2004 The MathWorks, Inc. 
%   $Revision: 2.12.2.4 $  $Date: 2004/12/06 16:37:04 $

if nargin < 3, 
    error('stats:binoinv:TooFewInputs','Requires three input arguments.'); 
end 

[errorcode y n p] = distchck(3,y,n,p);

if errorcode > 0
    error('stats:binoinv:InputSizeMismatch',...
          'Requires non-scalar arguments to match in size.');
end

k = 1:prod(size(y));
k1 = find(n < 0 | p < 0 | p > 1 | round(n) ~= n | y < 0 | y > 1 | isnan(y)); 
k2 = find(y == 1);

% Initialize X to 0.
if isa(y,'single') || isa(n,'single') || isa(p,'single')
  x = zeros(size(y),'single');
else
   x = zeros(size(y));
end
cumdist = x;
x(k2) = n(k2);
x(k1) = NaN;
k([k1(:); k2(:)]) = [];
if isempty(k), return; end

cumdist(k) = binopdf(0,n(k),p(k));

count = 0;

k = k(cumdist(k) < y(k));
while ~isempty(k)
   x(k) = x(k) + 1;
   count = count + 1;
   cumdist(k) = cumdist(k) + binopdf(count,n(k),p(k));
   k = k(cumdist(k) < y(k));
end