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