var x{1..4};

param best_val_found := 30.49654415;
param eps := 1; 		# = max(1, 1% x best_val_found)

subject to f:
	x[3]-1 +
	x[1]^2+
	x[2]^2 + (x[3]+x[4])^2 +
	sin(x[3])^2 + x[1]^2*x[2]^2 + x[4]-3 +
	sin(x[3])^2 +
	(x[4]-1)^2 +
	(x[2]^2)^2+
	(x[3]^2 + (x[4]+x[1])^2)^2 +
	(x[1]-4 + sin(x[4])^2 + x[2]^2*x[3]^2)^2 +
	sin(x[4])^4 <= best_val_found + eps;
subject to cons1:
	x[2] >= 1;
subject to cons2:
	-1D+10 <= x[3] <= 1;
subject to cons3:
	x[4] = 2;
subject to cons4:
	x[1]^2+x[2]^2 -1 <= 0;

solve;
#display f;
display x;