function y=humps(x)
y=1./(((x-1)-.3).^2 + .01) +1./(((x-1)-.9).^2 + .04) -6;