function f = t_pdf(t,k);
% f = t_pdf(t,k);
% student's t PDF

f = (exp (-(k + 1) * log (1 + t .^ 2 / k)/2) / (sqrt (k) * beta (k/2, 1/2)));

