function f = gammafn2(x,y)
f = -gamma(x)*(y/(1+y^2));