swh:1:snp:d91631a580497fff421f83ef3e706159a0d67602
Tip revision: ca0cbd288f03c7a29050b3d6698a96b45ccfa4b2 authored by Marko Laine on 15 August 2021, 17:59:37 UTC
cleanup
cleanup
Tip revision: ca0cbd2
tr.m
function y=tr(m,n,df)
%TR Random numbers from the t distribution
% TR(m,n,df) m,n shape of the result, df dergees of freedom
% Marko Laine <marko.laine@fmi.fi>
% $Revision: 1.6 $ $Date: 2012/09/27 11:47:40 $
if ~isunix & exist('distribs') == 3 % Windows and we have the mex code
y = distribs('tr',m,n,df);
else
% alternatively
z = norr(m,n);
x = chir(m,n,df);
y = z.*sqrt(df)./sqrt(x);
end