param N := 16;

var x{1..N};

minimize f:
sum {k in 1..15} (100*(x[k+1]-x[k]^2)^2+(1-x[k])^2);

data;
var x:=
1	-1.2
2	1
3	-1.2
4	1
5	-1.2
6	1
7	-1.2
8	1
9	-1.2
10	1
11	-1.2
12	1
13	-1.2
14	1
15	-1.2
16	1;

solve;

display x;