function y = objective(x)
y = (x(1)-1)^2 + 5*(x(2)-2)^2 + (x(3)-3)^2;