% Function to evaluate the second moment of the Chi-Square distribution

function f = chis(x);
global NU

f = x.^2 .* 0.5^(NU/2) .* x.^(NU/2-1) .* exp(-0.5 .* x) / gamma(NU/2)  ;
