constants
	r1 = 4.9;
	r2 = 3.0;
	c1 = 5.0;
	c2 = 1.0;
variables
	x in [-oo,oo];
	y in [-oo,oo];
constraints
	(x-c1)^2+(y-c2)^2 <= r1^2;
	(x-(c1-r1+r2))^2+(y-c2)^2 >= r2^2;
	(x-(c1+r2))^2+(y-c2)^2 >= (r1-r2)^2;
	y >= c2;
end