%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%
%
% Compute inverse of the cumulative F distribution function
% 
% CALL 
% 
% x = Finv(p,n1,n2);
% 
% 
% INPUT
% 
% - p, probability
% - n1, degree of freedom 1
% - n2, degree of freedom 2
% 
% References:
%   [1] William H. Press, Saul A. Teukolsky, William T. Vetterling, Brian
%   P. Flannery, NUMERICAL RECIPES, Cambridge University Press, 2007
% 
%
% L Ferraioli 06-12-2010
%
% $Id: Finv.m,v 1.1 2010/12/06 19:14:36 luigi Exp $
%
%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%
function x = Finv(p,n1,n2)

u = betaincinv(p,n1/2,n2/2);
x = n2.*u./(n1.*(1-u));

end