Noncentral chi-square random numbers
R = ncx2rnd(V,DELTA)
R = ncx2rnd(V,DELTA,m,n,...)
R = ncx2rnd(V,DELTA,[m,n,...])
R = ncx2rnd(V,DELTA) returns a matrix of random numbers chosen from the noncentral chi-square distribution using the corresponding degrees of freedom in V and positive noncentrality parameters in DELTA. V and DELTA can be vectors, matrices, or multidimensional arrays that have the same size, which is also the size of R. A scalar input for V or DELTA is expanded to a constant array with the same dimensions as the other input.
R = ncx2rnd(V,DELTA,m,n,...) or R = ncx2rnd(V,DELTA,[m,n,...]) generates an m-by-n-by-... array. The V, DELTA parameters can each be scalars or arrays of the same size as R.
ncx2rnd(4,2,6,3) ans = 6.8552 5.9650 11.2961 5.2631 4.2640 5.9495 9.1939 6.7162 3.8315 10.3100 4.4828 7.1653 2.1142 1.9826 4.6400 3.8852 5.3999 0.9282
 Evans, M., N. Hastings, and B. Peacock. Statistical Distributions. 2nd ed., Hoboken, NJ: John Wiley & Sons, Inc., 1993, pp. 50–52.
 Johnson, N., and S. Kotz. Distributions in Statistics: Continuous Univariate Distributions-2. Hoboken, NJ: John Wiley & Sons, Inc., 1970, pp. 130–148.