x = linspace(0,1);
y = 1./((x-0.3).^2+0.02)+1./((x-0.9).^2+0.04)-6;
plot(x,y)
grid
xlabel('x')
ylabel('y')
title('humps function')
