% Function to evaluate the first moment of the normal exponential

function f = gammadf(x);
global A B

f = x .* B^A .* x.^(A-1) .* exp(-B .* x) ./ gamma(A) ;
