gama.m

关于GP算法,不过好像有些问题.大家一起解决吧

function f=gama(x)
t=0:0.001:100;
y=t.^(x-1).*exp(-t);
f=trapz(y)*0.001;