CONFIG
 TrackType: 7;
END;
INPUT
 % add sphere 8
 variable_group v41,v42,v43,v51,v52,v53,v61,v62,v63,v71,v72,v73,v81,v82,v83;
 function f1,f2,f3,f4,f5,f6,f7,f8,f9,f10,f11;
 constant a,v11,v12,v13,v21,v22,v23,v31,v32,v33;

 % sqrt(3)
 a = 3^(1/2);

 % 3 fixed locations
 v11 = 0;
 v12 = 0;
 v13 = 0;

 v21 = 1;
 v22 = 0;
 v23 = 0;

 v31 = 1/2;
 v32 = a/2;
 v33 = 0;

 % D_{1,4} = (v41 - v11)^2 + (v42 - v12)^2 + (v43 - v13)^2 - 1;
 f1 = v41^2 + v42^2 + v43^2 - 1;

 % D_{2,4} - D_{1,4} = (v41 - v21)^2 + (v42 - v22)^2 + (v43 - v23)^2 - (v41 - v11)^2 - (v42 - v12)^2 - (v43 - v13)^2;
 f2 = 1 - 2*(v21*v41 + v22*v42 + v23*v43);

 % D_{1,5} = (v51 - v11)^2 + (v52 - v12)^2 + (v53 - v13)^2 - 1;
 f3 = v51^2 + v52^2 + v53^2 - 1;

 % D_{3,5} - D_{1,5} = (v51 - v31)^2 + (v52 - v32)^2 + (v53 - v33)^2 - (v51 - v11)^2 - (v52 - v12)^2 - (v53 - v13)^2;
 f4 = 1 - 2*(v31*v51 + v32*v52 + v33*v53);

 % D_{4,5} = (v51 - v41)^2 + (v52 - v42)^2 + (v53 - v43)^2 - 1;
 f5 = (v51 - v41)^2 + (v52 - v42)^2 + (v53 - v43)^2 - 1;

 % D_{2,6} = (v61 - v21)^2 + (v62 - v22)^2 + (v63 - v23)^2 - 1;
 f6 = (v61 - v21)^2 + (v62 - v22)^2 + (v63 - v23)^2 - 1;

 % D_{3,6} - D_{2,6} = (v61 - v31)^2 + (v62 - v32)^2 + (v63 - v33)^2 - (v61 - v21)^2 - (v62 - v22)^2 - (v63 - v23)^2;
 f7 = v61*v21 + v62*v22 + v63*v23 - v61*v31 - v62*v32 - v63*v33;

 % D_{4,6} = (v61 - v41)^2 + (v62 - v42)^2 + (v63 - v43)^2 - 1;
 f8 = (v61 - v41)^2 + (v62 - v42)^2 + (v63 - v43)^2 - 1;

 % D_{3,7} = (v71 - v31)^2 + (v72 - v32)^2 + (v73 - v33)^2 - 1;
 f9 = (v71 - v31)^2 + (v72 - v32)^2 + (v73 - v33)^2 - 1;

 % D_{5,7} = (v71 - v51)^2 + (v72 - v52)^2 + (v73 - v53)^2 - 1;
 f10= (v71 - v51)^2 + (v72 - v52)^2 + (v73 - v53)^2 - 1;

 % D_{6,7} = (v71 - v61)^2 + (v72 - v62)^2 + (v73 - v63)^2 - 1;
 f11= (v71 - v61)^2 + (v72 - v62)^2 + (v73 - v63)^2 - 1;

 % add distance constraint on 4 and 8, 5 and 8, and 6 and 8
 function g1,g2,g3;

 % D_{4,8} = (v81 - v41)^2 + (v82 - v42)^2 + (v83 - v43)^2 - 1;
 g1 = (v81 - v41)^2 + (v82 - v42)^2 + (v83 - v43)^2 - 1;

 % D_{5,8} = (v81 - v51)^2 + (v82 - v52)^2 + (v83 - v53)^2 - 1;
 g2 = (v81 - v51)^2 + (v82 - v52)^2 + (v83 - v53)^2 - 1;

 % D_{6,8} = (v81 - v61)^2 + (v82 - v62)^2 + (v83 - v63)^2 - 1;
 g3 = (v81 - v61)^2 + (v82 - v62)^2 + (v83 - v63)^2 - 1;

END;