param xinit{1..2};
var x{i in 1..2} := xinit[i];

minimize f:
	(x[2]-x[1]^2)^2/0.01+(x[1]-1)^2;

data;
param xinit:= 1 -1.2 2 1.0;

solve; display f; display x;