Noncentral F cumulative distribution function
p = ncfcdf(x,nu1,nu2,delta)
p = ncfcdf(x,nu1,nu2,delta,'upper')
p = ncfcdf(x,nu1,nu2,delta) computes the noncentral F cdf at each value in x using the corresponding numerator degrees of freedom in nu1, denominator degrees of freedom in nu2, and positive noncentrality parameters in delta. nu1, nu2, and delta can be vectors, matrices, or multidimensional arrays that have the same size, which is also the size of p. A scalar input for x, nu1, nu2, or delta is expanded to a constant array with the same dimensions as the other inputs.
p = ncfcdf(x,nu1,nu2,delta,'upper') returns the complement of the noncentral F cdf at each value in x, using an algorithm that more accurately computes the extreme upper tail probabilities.
The noncentral F cdf is
where I(x|a,b) is the incomplete beta function with parameters a and b.
 Johnson, N., and S. Kotz. Distributions in Statistics: Continuous Univariate Distributions-2. Hoboken, NJ: John Wiley & Sons, Inc., 1970, pp. 189–200.