function y = t_distribution(x,v)
%
A=gamma((v+1)/2);
A=A/gamma(v/2);
A=A/(sqrt(v*pi));
B=-(v+1)/2;
y=A*(1+(x.^2)/v).^B;