% Random chi-square variates

s = 23322;

nu = 6 ;

n = 5000;

z = zeros(n,1);

i = 1;

while i <= n
   
   out = chi_sqr(nu,s);
   
   z(i) = out(1);
   
   s = out(2);
   
   i = i + 1 ;
   
end