% Function to evaluate the area under the gamma distribution

% Passes back a single argument containing f(x) for the normal PDF function

function f = gammad(x);
global A B

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