param xinit{1..2}; var x{i in 1..2} := xinit[i]; minimize f: (x[1]^2)^2 + (10*x[1]/(x[1]+0.1)+2*x[2]^2)^2; data; param xinit:= 1 3 2 1; solve; display f; display x;