# #************************** # SET UP THE INITIAL DATA * #************************** # Problem : # -------- # Solving symmetric eigenvalue problems as systems of # nonlinear equations. # The problem is, given a symmetric matrix A, to find an orthogonal # matrix Q and diagonal matrix D such that A = Q(T) D Q. # Example B: a tridiagonal matrix with diagonals 2 and off diagonals -1 # Source: An idea by Nick Gould # Least-squares version # SIF input: Nick Gould, Nov 1992. # classification SUR2-AN-V-0 # The dimension of the matrix. #IE N 2 #IE N 100 # other parameter definitions # Define the upper triangular part of the matrix. # Define the eigenvalues # Define the eigenvectors # Introduce the eigen-equations Q(T) D Q - A = 0. param n := 10; param nm1 := -1 + (10); param a1_1 := 2.0; param jm1 := -1 + (10); param jm2 := -2 + (10); param a1_2 := -1.0; param a2_2 := 2.0; param a1_3 := 0.0; param a2_3 := -1.0; param a3_3 := 2.0; param a1_4 := 0.0; param a2_4 := 0.0; param a3_4 := -1.0; param a4_4 := 2.0; param a1_5 := 0.0; param a2_5 := 0.0; param a3_5 := 0.0; param a4_5 := -1.0; param a5_5 := 2.0; param a1_6 := 0.0; param a2_6 := 0.0; param a3_6 := 0.0; param a4_6 := 0.0; param a5_6 := -1.0; param a6_6 := 2.0; param a1_7 := 0.0; param a2_7 := 0.0; param a3_7 := 0.0; param a4_7 := 0.0; param a5_7 := 0.0; param a6_7 := -1.0; param a7_7 := 2.0; param a1_8 := 0.0; param a2_8 := 0.0; param a3_8 := 0.0; param a4_8 := 0.0; param a5_8 := 0.0; param a6_8 := 0.0; param a7_8 := -1.0; param a8_8 := 2.0; param a1_9 := 0.0; param a2_9 := 0.0; param a3_9 := 0.0; param a4_9 := 0.0; param a5_9 := 0.0; param a6_9 := 0.0; param a7_9 := 0.0; param a8_9 := -1.0; param a9_9 := 2.0; param a1_10 := 0.0; param a2_10 := 0.0; param a3_10 := 0.0; param a4_10 := 0.0; param a5_10 := 0.0; param a6_10 := 0.0; param a7_10 := 0.0; param a8_10 := 0.0; param a9_10 := -1.0; param a10_10 := 2.0; var d1 := 1.0; var q1_1 := 1.0; var q2_1; var q3_1; var q4_1; var q5_1; var q6_1; var q7_1; var q8_1; var q9_1; var q10_1; var d2 := 1.0; var q1_2; var q2_2 := 1.0; var q3_2; var q4_2; var q5_2; var q6_2; var q7_2; var q8_2; var q9_2; var q10_2; var d3 := 1.0; var q1_3; var q2_3; var q3_3 := 1.0; var q4_3; var q5_3; var q6_3; var q7_3; var q8_3; var q9_3; var q10_3; var d4 := 1.0; var q1_4; var q2_4; var q3_4; var q4_4 := 1.0; var q5_4; var q6_4; var q7_4; var q8_4; var q9_4; var q10_4; var d5 := 1.0; var q1_5; var q2_5; var q3_5; var q4_5; var q5_5 := 1.0; var q6_5; var q7_5; var q8_5; var q9_5; var q10_5; var d6 := 1.0; var q1_6; var q2_6; var q3_6; var q4_6; var q5_6; var q6_6 := 1.0; var q7_6; var q8_6; var q9_6; var q10_6; var d7 := 1.0; var q1_7; var q2_7; var q3_7; var q4_7; var q5_7; var q6_7; var q7_7 := 1.0; var q8_7; var q9_7; var q10_7; var d8 := 1.0; var q1_8; var q2_8; var q3_8; var q4_8; var q5_8; var q6_8; var q7_8; var q8_8 := 1.0; var q9_8; var q10_8; var d9 := 1.0; var q1_9; var q2_9; var q3_9; var q4_9; var q5_9; var q6_9; var q7_9; var q8_9; var q9_9 := 1.0; var q10_9; var d10 := 1.0; var q1_10; var q2_10; var q3_10; var q4_10; var q5_10; var q6_10; var q7_10; var q8_10; var q9_10; var q10_10 := 1.0; minimize obj: (q1_1 * q1_1 * d1 + q2_1 * q2_1 * d2 + q3_1 * q3_1 * d3 + q4_1 * q4_1 * d4 + q5_1 * q5_1 * d5 + q6_1 * q6_1 * d6 + q7_1 * q7_1 * d7 + q8_1 * q8_1 * d8 + q9_1 * q9_1 * d9 + q10_1 * q10_1 * d10 - 2.0)*(q1_1 * q1_1 * d1 + q2_1 * q2_1 * d2 + q3_1 * q3_1 * d3 + q4_1 * q4_1 * d4 + q5_1 * q5_1 * d5 + q6_1 * q6_1 * d6 + q7_1 * q7_1 * d7 + q8_1 * q8_1 * d8 + q9_1 * q9_1 * d9 + q10_1 * q10_1 * d10 - 2.0) + (q1_1 * q1_1 + q2_1 * q2_1 + q3_1 * q3_1 + q4_1 * q4_1 + q5_1 * q5_1 + q6_1 * q6_1 + q7_1 * q7_1 + q8_1 * q8_1 + q9_1 * q9_1 + q10_1 * q10_1 - 1.0)*(q1_1 * q1_1 + q2_1 * q2_1 + q3_1 * q3_1 + q4_1 * q4_1 + q5_1 * q5_1 + q6_1 * q6_1 + q7_1 * q7_1 + q8_1 * q8_1 + q9_1 * q9_1 + q10_1 * q10_1 - 1.0) + (q1_1 * q1_2 * d1 + q2_1 * q2_2 * d2 + q3_1 * q3_2 * d3 + q4_1 * q4_2 * d4 + q5_1 * q5_2 * d5 + q6_1 * q6_2 * d6 + q7_1 * q7_2 * d7 + q8_1 * q8_2 * d8 + q9_1 * q9_2 * d9 + q10_1 * q10_2 * d10 + 1.0)*(q1_1 * q1_2 * d1 + q2_1 * q2_2 * d2 + q3_1 * q3_2 * d3 + q4_1 * q4_2 * d4 + q5_1 * q5_2 * d5 + q6_1 * q6_2 * d6 + q7_1 * q7_2 * d7 + q8_1 * q8_2 * d8 + q9_1 * q9_2 * d9 + q10_1 * q10_2 * d10 + 1.0) + (q1_1 * q1_2 + q2_1 * q2_2 + q3_1 * q3_2 + q4_1 * q4_2 + q5_1 * q5_2 + q6_1 * q6_2 + q7_1 * q7_2 + q8_1 * q8_2 + q9_1 * q9_2 + q10_1 * q10_2)*(q1_1 * q1_2 + q2_1 * q2_2 + q3_1 * q3_2 + q4_1 * q4_2 + q5_1 * q5_2 + q6_1 * q6_2 + q7_1 * q7_2 + q8_1 * q8_2 + q9_1 * q9_2 + q10_1 * q10_2) + (q1_2 * q1_2 * d1 + q2_2 * q2_2 * d2 + q3_2 * q3_2 * d3 + q4_2 * q4_2 * d4 + q5_2 * q5_2 * d5 + q6_2 * q6_2 * d6 + q7_2 * q7_2 * d7 + q8_2 * q8_2 * d8 + q9_2 * q9_2 * d9 + q10_2 * q10_2 * d10 - 2.0)*(q1_2 * q1_2 * d1 + q2_2 * q2_2 * d2 + q3_2 * q3_2 * d3 + q4_2 * q4_2 * d4 + q5_2 * q5_2 * d5 + q6_2 * q6_2 * d6 + q7_2 * q7_2 * d7 + q8_2 * q8_2 * d8 + q9_2 * q9_2 * d9 + q10_2 * q10_2 * d10 - 2.0) + (q1_2 * q1_2 + q2_2 * q2_2 + q3_2 * q3_2 + q4_2 * q4_2 + q5_2 * q5_2 + q6_2 * q6_2 + q7_2 * q7_2 + q8_2 * q8_2 + q9_2 * q9_2 + q10_2 * q10_2 - 1.0)*(q1_2 * q1_2 + q2_2 * q2_2 + q3_2 * q3_2 + q4_2 * q4_2 + q5_2 * q5_2 + q6_2 * q6_2 + q7_2 * q7_2 + q8_2 * q8_2 + q9_2 * q9_2 + q10_2 * q10_2 - 1.0) + (q1_1 * q1_3 * d1 + q2_1 * q2_3 * d2 + q3_1 * q3_3 * d3 + q4_1 * q4_3 * d4 + q5_1 * q5_3 * d5 + q6_1 * q6_3 * d6 + q7_1 * q7_3 * d7 + q8_1 * q8_3 * d8 + q9_1 * q9_3 * d9 + q10_1 * q10_3 * d10)*(q1_1 * q1_3 * d1 + q2_1 * q2_3 * d2 + q3_1 * q3_3 * d3 + q4_1 * q4_3 * d4 + q5_1 * q5_3 * d5 + q6_1 * q6_3 * d6 + q7_1 * q7_3 * d7 + q8_1 * q8_3 * d8 + q9_1 * q9_3 * d9 + q10_1 * q10_3 * d10) + (q1_1 * q1_3 + q2_1 * q2_3 + q3_1 * q3_3 + q4_1 * q4_3 + q5_1 * q5_3 + q6_1 * q6_3 + q7_1 * q7_3 + q8_1 * q8_3 + q9_1 * q9_3 + q10_1 * q10_3)*(q1_1 * q1_3 + q2_1 * q2_3 + q3_1 * q3_3 + q4_1 * q4_3 + q5_1 * q5_3 + q6_1 * q6_3 + q7_1 * q7_3 + q8_1 * q8_3 + q9_1 * q9_3 + q10_1 * q10_3) + (q1_2 * q1_3 * d1 + q2_2 * q2_3 * d2 + q3_2 * q3_3 * d3 + q4_2 * q4_3 * d4 + q5_2 * q5_3 * d5 + q6_2 * q6_3 * d6 + q7_2 * q7_3 * d7 + q8_2 * q8_3 * d8 + q9_2 * q9_3 * d9 + q10_2 * q10_3 * d10 + 1.0)*(q1_2 * q1_3 * d1 + q2_2 * q2_3 * d2 + q3_2 * q3_3 * d3 + q4_2 * q4_3 * d4 + q5_2 * q5_3 * d5 + q6_2 * q6_3 * d6 + q7_2 * q7_3 * d7 + q8_2 * q8_3 * d8 + q9_2 * q9_3 * d9 + q10_2 * q10_3 * d10 + 1.0) + (q1_2 * q1_3 + q2_2 * q2_3 + q3_2 * q3_3 + q4_2 * q4_3 + q5_2 * q5_3 + q6_2 * q6_3 + q7_2 * q7_3 + q8_2 * q8_3 + q9_2 * q9_3 + q10_2 * q10_3)*(q1_2 * q1_3 + q2_2 * q2_3 + q3_2 * q3_3 + q4_2 * q4_3 + q5_2 * q5_3 + q6_2 * q6_3 + q7_2 * q7_3 + q8_2 * q8_3 + q9_2 * q9_3 + q10_2 * q10_3) + (q1_3 * q1_3 * d1 + q2_3 * q2_3 * d2 + q3_3 * q3_3 * d3 + q4_3 * q4_3 * d4 + q5_3 * q5_3 * d5 + q6_3 * q6_3 * d6 + q7_3 * q7_3 * d7 + q8_3 * q8_3 * d8 + q9_3 * q9_3 * d9 + q10_3 * q10_3 * d10 - 2.0)*(q1_3 * q1_3 * d1 + q2_3 * q2_3 * d2 + q3_3 * q3_3 * d3 + q4_3 * q4_3 * d4 + q5_3 * q5_3 * d5 + q6_3 * q6_3 * d6 + q7_3 * q7_3 * d7 + q8_3 * q8_3 * d8 + q9_3 * q9_3 * d9 + q10_3 * q10_3 * d10 - 2.0) + (q1_3 * q1_3 + q2_3 * q2_3 + q3_3 * q3_3 + q4_3 * q4_3 + q5_3 * q5_3 + q6_3 * q6_3 + q7_3 * q7_3 + q8_3 * q8_3 + q9_3 * q9_3 + q10_3 * q10_3 - 1.0)*(q1_3 * q1_3 + q2_3 * q2_3 + q3_3 * q3_3 + q4_3 * q4_3 + q5_3 * q5_3 + q6_3 * q6_3 + q7_3 * q7_3 + q8_3 * q8_3 + q9_3 * q9_3 + q10_3 * q10_3 - 1.0) + (q1_1 * q1_4 * d1 + q2_1 * q2_4 * d2 + q3_1 * q3_4 * d3 + q4_1 * q4_4 * d4 + q5_1 * q5_4 * d5 + q6_1 * q6_4 * d6 + q7_1 * q7_4 * d7 + q8_1 * q8_4 * d8 + q9_1 * q9_4 * d9 + q10_1 * q10_4 * d10)*(q1_1 * q1_4 * d1 + q2_1 * q2_4 * d2 + q3_1 * q3_4 * d3 + q4_1 * q4_4 * d4 + q5_1 * q5_4 * d5 + q6_1 * q6_4 * d6 + q7_1 * q7_4 * d7 + q8_1 * q8_4 * d8 + q9_1 * q9_4 * d9 + q10_1 * q10_4 * d10) + (q1_1 * q1_4 + q2_1 * q2_4 + q3_1 * q3_4 + q4_1 * q4_4 + q5_1 * q5_4 + q6_1 * q6_4 + q7_1 * q7_4 + q8_1 * q8_4 + q9_1 * q9_4 + q10_1 * q10_4)*(q1_1 * q1_4 + q2_1 * q2_4 + q3_1 * q3_4 + q4_1 * q4_4 + q5_1 * q5_4 + q6_1 * q6_4 + q7_1 * q7_4 + q8_1 * q8_4 + q9_1 * q9_4 + q10_1 * q10_4) + (q1_2 * q1_4 * d1 + q2_2 * q2_4 * d2 + q3_2 * q3_4 * d3 + q4_2 * q4_4 * d4 + q5_2 * q5_4 * d5 + q6_2 * q6_4 * d6 + q7_2 * q7_4 * d7 + q8_2 * q8_4 * d8 + q9_2 * q9_4 * d9 + q10_2 * q10_4 * d10)*(q1_2 * q1_4 * d1 + q2_2 * q2_4 * d2 + q3_2 * q3_4 * d3 + q4_2 * q4_4 * d4 + q5_2 * q5_4 * d5 + q6_2 * q6_4 * d6 + q7_2 * q7_4 * d7 + q8_2 * q8_4 * d8 + q9_2 * q9_4 * d9 + q10_2 * q10_4 * d10) + (q1_2 * q1_4 + q2_2 * q2_4 + q3_2 * q3_4 + q4_2 * q4_4 + q5_2 * q5_4 + q6_2 * q6_4 + q7_2 * q7_4 + q8_2 * q8_4 + q9_2 * q9_4 + q10_2 * q10_4)*(q1_2 * q1_4 + q2_2 * q2_4 + q3_2 * q3_4 + q4_2 * q4_4 + q5_2 * q5_4 + q6_2 * q6_4 + q7_2 * q7_4 + q8_2 * q8_4 + q9_2 * q9_4 + q10_2 * q10_4) + (q1_3 * q1_4 * d1 + q2_3 * q2_4 * d2 + q3_3 * q3_4 * d3 + q4_3 * q4_4 * d4 + q5_3 * q5_4 * d5 + q6_3 * q6_4 * d6 + q7_3 * q7_4 * d7 + q8_3 * q8_4 * d8 + q9_3 * q9_4 * d9 + q10_3 * q10_4 * d10 + 1.0)*(q1_3 * q1_4 * d1 + q2_3 * q2_4 * d2 + q3_3 * q3_4 * d3 + q4_3 * q4_4 * d4 + q5_3 * q5_4 * d5 + q6_3 * q6_4 * d6 + q7_3 * q7_4 * d7 + q8_3 * q8_4 * d8 + q9_3 * q9_4 * d9 + q10_3 * q10_4 * d10 + 1.0) + (q1_3 * q1_4 + q2_3 * q2_4 + q3_3 * q3_4 + q4_3 * q4_4 + q5_3 * q5_4 + q6_3 * q6_4 + q7_3 * q7_4 + q8_3 * q8_4 + q9_3 * q9_4 + q10_3 * q10_4)*(q1_3 * q1_4 + q2_3 * q2_4 + q3_3 * q3_4 + q4_3 * q4_4 + q5_3 * q5_4 + q6_3 * q6_4 + q7_3 * q7_4 + q8_3 * q8_4 + q9_3 * q9_4 + q10_3 * q10_4) + (q1_4 * q1_4 * d1 + q2_4 * q2_4 * d2 + q3_4 * q3_4 * d3 + q4_4 * q4_4 * d4 + q5_4 * q5_4 * d5 + q6_4 * q6_4 * d6 + q7_4 * q7_4 * d7 + q8_4 * q8_4 * d8 + q9_4 * q9_4 * d9 + q10_4 * q10_4 * d10 - 2.0)*(q1_4 * q1_4 * d1 + q2_4 * q2_4 * d2 + q3_4 * q3_4 * d3 + q4_4 * q4_4 * d4 + q5_4 * q5_4 * d5 + q6_4 * q6_4 * d6 + q7_4 * q7_4 * d7 + q8_4 * q8_4 * d8 + q9_4 * q9_4 * d9 + q10_4 * q10_4 * d10 - 2.0) + (q1_4 * q1_4 + q2_4 * q2_4 + q3_4 * q3_4 + q4_4 * q4_4 + q5_4 * q5_4 + q6_4 * q6_4 + q7_4 * q7_4 + q8_4 * q8_4 + q9_4 * q9_4 + q10_4 * q10_4 - 1.0)*(q1_4 * q1_4 + q2_4 * q2_4 + q3_4 * q3_4 + q4_4 * q4_4 + q5_4 * q5_4 + q6_4 * q6_4 + q7_4 * q7_4 + q8_4 * q8_4 + q9_4 * q9_4 + q10_4 * q10_4 - 1.0) + (q1_1 * q1_5 * d1 + q2_1 * q2_5 * d2 + q3_1 * q3_5 * d3 + q4_1 * q4_5 * d4 + q5_1 * q5_5 * d5 + q6_1 * q6_5 * d6 + q7_1 * q7_5 * d7 + q8_1 * q8_5 * d8 + q9_1 * q9_5 * d9 + q10_1 * q10_5 * d10)*(q1_1 * q1_5 * d1 + q2_1 * q2_5 * d2 + q3_1 * q3_5 * d3 + q4_1 * q4_5 * d4 + q5_1 * q5_5 * d5 + q6_1 * q6_5 * d6 + q7_1 * q7_5 * d7 + q8_1 * q8_5 * d8 + q9_1 * q9_5 * d9 + q10_1 * q10_5 * d10) + (q1_1 * q1_5 + q2_1 * q2_5 + q3_1 * q3_5 + q4_1 * q4_5 + q5_1 * q5_5 + q6_1 * q6_5 + q7_1 * q7_5 + q8_1 * q8_5 + q9_1 * q9_5 + q10_1 * q10_5)*(q1_1 * q1_5 + q2_1 * q2_5 + q3_1 * q3_5 + q4_1 * q4_5 + q5_1 * q5_5 + q6_1 * q6_5 + q7_1 * q7_5 + q8_1 * q8_5 + q9_1 * q9_5 + q10_1 * q10_5) + (q1_2 * q1_5 * d1 + q2_2 * q2_5 * d2 + q3_2 * q3_5 * d3 + q4_2 * q4_5 * d4 + q5_2 * q5_5 * d5 + q6_2 * q6_5 * d6 + q7_2 * q7_5 * d7 + q8_2 * q8_5 * d8 + q9_2 * q9_5 * d9 + q10_2 * q10_5 * d10)*(q1_2 * q1_5 * d1 + q2_2 * q2_5 * d2 + q3_2 * q3_5 * d3 + q4_2 * q4_5 * d4 + q5_2 * q5_5 * d5 + q6_2 * q6_5 * d6 + q7_2 * q7_5 * d7 + q8_2 * q8_5 * d8 + q9_2 * q9_5 * d9 + q10_2 * q10_5 * d10) + (q1_2 * q1_5 + q2_2 * q2_5 + q3_2 * q3_5 + q4_2 * q4_5 + q5_2 * q5_5 + q6_2 * q6_5 + q7_2 * q7_5 + q8_2 * q8_5 + q9_2 * q9_5 + q10_2 * q10_5)*(q1_2 * q1_5 + q2_2 * q2_5 + q3_2 * q3_5 + q4_2 * q4_5 + q5_2 * q5_5 + q6_2 * q6_5 + q7_2 * q7_5 + q8_2 * q8_5 + q9_2 * q9_5 + q10_2 * q10_5) + (q1_3 * q1_5 * d1 + q2_3 * q2_5 * d2 + q3_3 * q3_5 * d3 + q4_3 * q4_5 * d4 + q5_3 * q5_5 * d5 + q6_3 * q6_5 * d6 + q7_3 * q7_5 * d7 + q8_3 * q8_5 * d8 + q9_3 * q9_5 * d9 + q10_3 * q10_5 * d10)*(q1_3 * q1_5 * d1 + q2_3 * q2_5 * d2 + q3_3 * q3_5 * d3 + q4_3 * q4_5 * d4 + q5_3 * q5_5 * d5 + q6_3 * q6_5 * d6 + q7_3 * q7_5 * d7 + q8_3 * q8_5 * d8 + q9_3 * q9_5 * d9 + q10_3 * q10_5 * d10) + (q1_3 * q1_5 + q2_3 * q2_5 + q3_3 * q3_5 + q4_3 * q4_5 + q5_3 * q5_5 + q6_3 * q6_5 + q7_3 * q7_5 + q8_3 * q8_5 + q9_3 * q9_5 + q10_3 * q10_5)*(q1_3 * q1_5 + q2_3 * q2_5 + q3_3 * q3_5 + q4_3 * q4_5 + q5_3 * q5_5 + q6_3 * q6_5 + q7_3 * q7_5 + q8_3 * q8_5 + q9_3 * q9_5 + q10_3 * q10_5) + (q1_4 * q1_5 * d1 + q2_4 * q2_5 * d2 + q3_4 * q3_5 * d3 + q4_4 * q4_5 * d4 + q5_4 * q5_5 * d5 + q6_4 * q6_5 * d6 + q7_4 * q7_5 * d7 + q8_4 * q8_5 * d8 + q9_4 * q9_5 * d9 + q10_4 * q10_5 * d10 + 1.0)*(q1_4 * q1_5 * d1 + q2_4 * q2_5 * d2 + q3_4 * q3_5 * d3 + q4_4 * q4_5 * d4 + q5_4 * q5_5 * d5 + q6_4 * q6_5 * d6 + q7_4 * q7_5 * d7 + q8_4 * q8_5 * d8 + q9_4 * q9_5 * d9 + q10_4 * q10_5 * d10 + 1.0) + (q1_4 * q1_5 + q2_4 * q2_5 + q3_4 * q3_5 + q4_4 * q4_5 + q5_4 * q5_5 + q6_4 * q6_5 + q7_4 * q7_5 + q8_4 * q8_5 + q9_4 * q9_5 + q10_4 * q10_5)*(q1_4 * q1_5 + q2_4 * q2_5 + q3_4 * q3_5 + q4_4 * q4_5 + q5_4 * q5_5 + q6_4 * q6_5 + q7_4 * q7_5 + q8_4 * q8_5 + q9_4 * q9_5 + q10_4 * q10_5) + (q1_5 * q1_5 * d1 + q2_5 * q2_5 * d2 + q3_5 * q3_5 * d3 + q4_5 * q4_5 * d4 + q5_5 * q5_5 * d5 + q6_5 * q6_5 * d6 + q7_5 * q7_5 * d7 + q8_5 * q8_5 * d8 + q9_5 * q9_5 * d9 + q10_5 * q10_5 * d10 - 2.0)*(q1_5 * q1_5 * d1 + q2_5 * q2_5 * d2 + q3_5 * q3_5 * d3 + q4_5 * q4_5 * d4 + q5_5 * q5_5 * d5 + q6_5 * q6_5 * d6 + q7_5 * q7_5 * d7 + q8_5 * q8_5 * d8 + q9_5 * q9_5 * d9 + q10_5 * q10_5 * d10 - 2.0) + (q1_5 * q1_5 + q2_5 * q2_5 + q3_5 * q3_5 + q4_5 * q4_5 + q5_5 * q5_5 + q6_5 * q6_5 + q7_5 * q7_5 + q8_5 * q8_5 + q9_5 * q9_5 + q10_5 * q10_5 - 1.0)*(q1_5 * q1_5 + q2_5 * q2_5 + q3_5 * q3_5 + q4_5 * q4_5 + q5_5 * q5_5 + q6_5 * q6_5 + q7_5 * q7_5 + q8_5 * q8_5 + q9_5 * q9_5 + q10_5 * q10_5 - 1.0) + (q1_1 * q1_6 * d1 + q2_1 * q2_6 * d2 + q3_1 * q3_6 * d3 + q4_1 * q4_6 * d4 + q5_1 * q5_6 * d5 + q6_1 * q6_6 * d6 + q7_1 * q7_6 * d7 + q8_1 * q8_6 * d8 + q9_1 * q9_6 * d9 + q10_1 * q10_6 * d10)*(q1_1 * q1_6 * d1 + q2_1 * q2_6 * d2 + q3_1 * q3_6 * d3 + q4_1 * q4_6 * d4 + q5_1 * q5_6 * d5 + q6_1 * q6_6 * d6 + q7_1 * q7_6 * d7 + q8_1 * q8_6 * d8 + q9_1 * q9_6 * d9 + q10_1 * q10_6 * d10) + (q1_1 * q1_6 + q2_1 * q2_6 + q3_1 * q3_6 + q4_1 * q4_6 + q5_1 * q5_6 + q6_1 * q6_6 + q7_1 * q7_6 + q8_1 * q8_6 + q9_1 * q9_6 + q10_1 * q10_6)*(q1_1 * q1_6 + q2_1 * q2_6 + q3_1 * q3_6 + q4_1 * q4_6 + q5_1 * q5_6 + q6_1 * q6_6 + q7_1 * q7_6 + q8_1 * q8_6 + q9_1 * q9_6 + q10_1 * q10_6) + (q1_2 * q1_6 * d1 + q2_2 * q2_6 * d2 + q3_2 * q3_6 * d3 + q4_2 * q4_6 * d4 + q5_2 * q5_6 * d5 + q6_2 * q6_6 * d6 + q7_2 * q7_6 * d7 + q8_2 * q8_6 * d8 + q9_2 * q9_6 * d9 + q10_2 * q10_6 * d10)*(q1_2 * q1_6 * d1 + q2_2 * q2_6 * d2 + q3_2 * q3_6 * d3 + q4_2 * q4_6 * d4 + q5_2 * q5_6 * d5 + q6_2 * q6_6 * d6 + q7_2 * q7_6 * d7 + q8_2 * q8_6 * d8 + q9_2 * q9_6 * d9 + q10_2 * q10_6 * d10) + (q1_2 * q1_6 + q2_2 * q2_6 + q3_2 * q3_6 + q4_2 * q4_6 + q5_2 * q5_6 + q6_2 * q6_6 + q7_2 * q7_6 + q8_2 * q8_6 + q9_2 * q9_6 + q10_2 * q10_6)*(q1_2 * q1_6 + q2_2 * q2_6 + q3_2 * q3_6 + q4_2 * q4_6 + q5_2 * q5_6 + q6_2 * q6_6 + q7_2 * q7_6 + q8_2 * q8_6 + q9_2 * q9_6 + q10_2 * q10_6) + (q1_3 * q1_6 * d1 + q2_3 * q2_6 * d2 + q3_3 * q3_6 * d3 + q4_3 * q4_6 * d4 + q5_3 * q5_6 * d5 + q6_3 * q6_6 * d6 + q7_3 * q7_6 * d7 + q8_3 * q8_6 * d8 + q9_3 * q9_6 * d9 + q10_3 * q10_6 * d10)*(q1_3 * q1_6 * d1 + q2_3 * q2_6 * d2 + q3_3 * q3_6 * d3 + q4_3 * q4_6 * d4 + q5_3 * q5_6 * d5 + q6_3 * q6_6 * d6 + q7_3 * q7_6 * d7 + q8_3 * q8_6 * d8 + q9_3 * q9_6 * d9 + q10_3 * q10_6 * d10) + (q1_3 * q1_6 + q2_3 * q2_6 + q3_3 * q3_6 + q4_3 * q4_6 + q5_3 * q5_6 + q6_3 * q6_6 + q7_3 * q7_6 + q8_3 * q8_6 + q9_3 * q9_6 + q10_3 * q10_6)*(q1_3 * q1_6 + q2_3 * q2_6 + q3_3 * q3_6 + q4_3 * q4_6 + q5_3 * q5_6 + q6_3 * q6_6 + q7_3 * q7_6 + q8_3 * q8_6 + q9_3 * q9_6 + q10_3 * q10_6) + (q1_4 * q1_6 * d1 + q2_4 * q2_6 * d2 + q3_4 * q3_6 * d3 + q4_4 * q4_6 * d4 + q5_4 * q5_6 * d5 + q6_4 * q6_6 * d6 + q7_4 * q7_6 * d7 + q8_4 * q8_6 * d8 + q9_4 * q9_6 * d9 + q10_4 * q10_6 * d10)*(q1_4 * q1_6 * d1 + q2_4 * q2_6 * d2 + q3_4 * q3_6 * d3 + q4_4 * q4_6 * d4 + q5_4 * q5_6 * d5 + q6_4 * q6_6 * d6 + q7_4 * q7_6 * d7 + q8_4 * q8_6 * d8 + q9_4 * q9_6 * d9 + q10_4 * q10_6 * d10) + (q1_4 * q1_6 + q2_4 * q2_6 + q3_4 * q3_6 + q4_4 * q4_6 + q5_4 * q5_6 + q6_4 * q6_6 + q7_4 * q7_6 + q8_4 * q8_6 + q9_4 * q9_6 + q10_4 * q10_6)*(q1_4 * q1_6 + q2_4 * q2_6 + q3_4 * q3_6 + q4_4 * q4_6 + q5_4 * q5_6 + q6_4 * q6_6 + q7_4 * q7_6 + q8_4 * q8_6 + q9_4 * q9_6 + q10_4 * q10_6) + (q1_5 * q1_6 * d1 + q2_5 * q2_6 * d2 + q3_5 * q3_6 * d3 + q4_5 * q4_6 * d4 + q5_5 * q5_6 * d5 + q6_5 * q6_6 * d6 + q7_5 * q7_6 * d7 + q8_5 * q8_6 * d8 + q9_5 * q9_6 * d9 + q10_5 * q10_6 * d10 + 1.0)*(q1_5 * q1_6 * d1 + q2_5 * q2_6 * d2 + q3_5 * q3_6 * d3 + q4_5 * q4_6 * d4 + q5_5 * q5_6 * d5 + q6_5 * q6_6 * d6 + q7_5 * q7_6 * d7 + q8_5 * q8_6 * d8 + q9_5 * q9_6 * d9 + q10_5 * q10_6 * d10 + 1.0) + (q1_5 * q1_6 + q2_5 * q2_6 + q3_5 * q3_6 + q4_5 * q4_6 + q5_5 * q5_6 + q6_5 * q6_6 + q7_5 * q7_6 + q8_5 * q8_6 + q9_5 * q9_6 + q10_5 * q10_6)*(q1_5 * q1_6 + q2_5 * q2_6 + q3_5 * q3_6 + q4_5 * q4_6 + q5_5 * q5_6 + q6_5 * q6_6 + q7_5 * q7_6 + q8_5 * q8_6 + q9_5 * q9_6 + q10_5 * q10_6) + (q1_6 * q1_6 * d1 + q2_6 * q2_6 * d2 + q3_6 * q3_6 * d3 + q4_6 * q4_6 * d4 + q5_6 * q5_6 * d5 + q6_6 * q6_6 * d6 + q7_6 * q7_6 * d7 + q8_6 * q8_6 * d8 + q9_6 * q9_6 * d9 + q10_6 * q10_6 * d10 - 2.0)*(q1_6 * q1_6 * d1 + q2_6 * q2_6 * d2 + q3_6 * q3_6 * d3 + q4_6 * q4_6 * d4 + q5_6 * q5_6 * d5 + q6_6 * q6_6 * d6 + q7_6 * q7_6 * d7 + q8_6 * q8_6 * d8 + q9_6 * q9_6 * d9 + q10_6 * q10_6 * d10 - 2.0) + (q1_6 * q1_6 + q2_6 * q2_6 + q3_6 * q3_6 + q4_6 * q4_6 + q5_6 * q5_6 + q6_6 * q6_6 + q7_6 * q7_6 + q8_6 * q8_6 + q9_6 * q9_6 + q10_6 * q10_6 - 1.0)*(q1_6 * q1_6 + q2_6 * q2_6 + q3_6 * q3_6 + q4_6 * q4_6 + q5_6 * q5_6 + q6_6 * q6_6 + q7_6 * q7_6 + q8_6 * q8_6 + q9_6 * q9_6 + q10_6 * q10_6 - 1.0) + (q1_1 * q1_7 * d1 + q2_1 * q2_7 * d2 + q3_1 * q3_7 * d3 + q4_1 * q4_7 * d4 + q5_1 * q5_7 * d5 + q6_1 * q6_7 * d6 + q7_1 * q7_7 * d7 + q8_1 * q8_7 * d8 + q9_1 * q9_7 * d9 + q10_1 * q10_7 * d10)*(q1_1 * q1_7 * d1 + q2_1 * q2_7 * d2 + q3_1 * q3_7 * d3 + q4_1 * q4_7 * d4 + q5_1 * q5_7 * d5 + q6_1 * q6_7 * d6 + q7_1 * q7_7 * d7 + q8_1 * q8_7 * d8 + q9_1 * q9_7 * d9 + q10_1 * q10_7 * d10) + (q1_1 * q1_7 + q2_1 * q2_7 + q3_1 * q3_7 + q4_1 * q4_7 + q5_1 * q5_7 + q6_1 * q6_7 + q7_1 * q7_7 + q8_1 * q8_7 + q9_1 * q9_7 + q10_1 * q10_7)*(q1_1 * q1_7 + q2_1 * q2_7 + q3_1 * q3_7 + q4_1 * q4_7 + q5_1 * q5_7 + q6_1 * q6_7 + q7_1 * q7_7 + q8_1 * q8_7 + q9_1 * q9_7 + q10_1 * q10_7) + (q1_2 * q1_7 * d1 + q2_2 * q2_7 * d2 + q3_2 * q3_7 * d3 + q4_2 * q4_7 * d4 + q5_2 * q5_7 * d5 + q6_2 * q6_7 * d6 + q7_2 * q7_7 * d7 + q8_2 * q8_7 * d8 + q9_2 * q9_7 * d9 + q10_2 * q10_7 * d10)*(q1_2 * q1_7 * d1 + q2_2 * q2_7 * d2 + q3_2 * q3_7 * d3 + q4_2 * q4_7 * d4 + q5_2 * q5_7 * d5 + q6_2 * q6_7 * d6 + q7_2 * q7_7 * d7 + q8_2 * q8_7 * d8 + q9_2 * q9_7 * d9 + q10_2 * q10_7 * d10) + (q1_2 * q1_7 + q2_2 * q2_7 + q3_2 * q3_7 + q4_2 * q4_7 + q5_2 * q5_7 + q6_2 * q6_7 + q7_2 * q7_7 + q8_2 * q8_7 + q9_2 * q9_7 + q10_2 * q10_7)*(q1_2 * q1_7 + q2_2 * q2_7 + q3_2 * q3_7 + q4_2 * q4_7 + q5_2 * q5_7 + q6_2 * q6_7 + q7_2 * q7_7 + q8_2 * q8_7 + q9_2 * q9_7 + q10_2 * q10_7) + (q1_3 * q1_7 * d1 + q2_3 * q2_7 * d2 + q3_3 * q3_7 * d3 + q4_3 * q4_7 * d4 + q5_3 * q5_7 * d5 + q6_3 * q6_7 * d6 + q7_3 * q7_7 * d7 + q8_3 * q8_7 * d8 + q9_3 * q9_7 * d9 + q10_3 * q10_7 * d10)*(q1_3 * q1_7 * d1 + q2_3 * q2_7 * d2 + q3_3 * q3_7 * d3 + q4_3 * q4_7 * d4 + q5_3 * q5_7 * d5 + q6_3 * q6_7 * d6 + q7_3 * q7_7 * d7 + q8_3 * q8_7 * d8 + q9_3 * q9_7 * d9 + q10_3 * q10_7 * d10) + (q1_3 * q1_7 + q2_3 * q2_7 + q3_3 * q3_7 + q4_3 * q4_7 + q5_3 * q5_7 + q6_3 * q6_7 + q7_3 * q7_7 + q8_3 * q8_7 + q9_3 * q9_7 + q10_3 * q10_7)*(q1_3 * q1_7 + q2_3 * q2_7 + q3_3 * q3_7 + q4_3 * q4_7 + q5_3 * q5_7 + q6_3 * q6_7 + q7_3 * q7_7 + q8_3 * q8_7 + q9_3 * q9_7 + q10_3 * q10_7) + (q1_4 * q1_7 * d1 + q2_4 * q2_7 * d2 + q3_4 * q3_7 * d3 + q4_4 * q4_7 * d4 + q5_4 * q5_7 * d5 + q6_4 * q6_7 * d6 + q7_4 * q7_7 * d7 + q8_4 * q8_7 * d8 + q9_4 * q9_7 * d9 + q10_4 * q10_7 * d10)*(q1_4 * q1_7 * d1 + q2_4 * q2_7 * d2 + q3_4 * q3_7 * d3 + q4_4 * q4_7 * d4 + q5_4 * q5_7 * d5 + q6_4 * q6_7 * d6 + q7_4 * q7_7 * d7 + q8_4 * q8_7 * d8 + q9_4 * q9_7 * d9 + q10_4 * q10_7 * d10) + (q1_4 * q1_7 + q2_4 * q2_7 + q3_4 * q3_7 + q4_4 * q4_7 + q5_4 * q5_7 + q6_4 * q6_7 + q7_4 * q7_7 + q8_4 * q8_7 + q9_4 * q9_7 + q10_4 * q10_7)*(q1_4 * q1_7 + q2_4 * q2_7 + q3_4 * q3_7 + q4_4 * q4_7 + q5_4 * q5_7 + q6_4 * q6_7 + q7_4 * q7_7 + q8_4 * q8_7 + q9_4 * q9_7 + q10_4 * q10_7) + (q1_5 * q1_7 * d1 + q2_5 * q2_7 * d2 + q3_5 * q3_7 * d3 + q4_5 * q4_7 * d4 + q5_5 * q5_7 * d5 + q6_5 * q6_7 * d6 + q7_5 * q7_7 * d7 + q8_5 * q8_7 * d8 + q9_5 * q9_7 * d9 + q10_5 * q10_7 * d10)*(q1_5 * q1_7 * d1 + q2_5 * q2_7 * d2 + q3_5 * q3_7 * d3 + q4_5 * q4_7 * d4 + q5_5 * q5_7 * d5 + q6_5 * q6_7 * d6 + q7_5 * q7_7 * d7 + q8_5 * q8_7 * d8 + q9_5 * q9_7 * d9 + q10_5 * q10_7 * d10) + (q1_5 * q1_7 + q2_5 * q2_7 + q3_5 * q3_7 + q4_5 * q4_7 + q5_5 * q5_7 + q6_5 * q6_7 + q7_5 * q7_7 + q8_5 * q8_7 + q9_5 * q9_7 + q10_5 * q10_7)*(q1_5 * q1_7 + q2_5 * q2_7 + q3_5 * q3_7 + q4_5 * q4_7 + q5_5 * q5_7 + q6_5 * q6_7 + q7_5 * q7_7 + q8_5 * q8_7 + q9_5 * q9_7 + q10_5 * q10_7) + (q1_6 * q1_7 * d1 + q2_6 * q2_7 * d2 + q3_6 * q3_7 * d3 + q4_6 * q4_7 * d4 + q5_6 * q5_7 * d5 + q6_6 * q6_7 * d6 + q7_6 * q7_7 * d7 + q8_6 * q8_7 * d8 + q9_6 * q9_7 * d9 + q10_6 * q10_7 * d10 + 1.0)*(q1_6 * q1_7 * d1 + q2_6 * q2_7 * d2 + q3_6 * q3_7 * d3 + q4_6 * q4_7 * d4 + q5_6 * q5_7 * d5 + q6_6 * q6_7 * d6 + q7_6 * q7_7 * d7 + q8_6 * q8_7 * d8 + q9_6 * q9_7 * d9 + q10_6 * q10_7 * d10 + 1.0) + (q1_6 * q1_7 + q2_6 * q2_7 + q3_6 * q3_7 + q4_6 * q4_7 + q5_6 * q5_7 + q6_6 * q6_7 + q7_6 * q7_7 + q8_6 * q8_7 + q9_6 * q9_7 + q10_6 * q10_7)*(q1_6 * q1_7 + q2_6 * q2_7 + q3_6 * q3_7 + q4_6 * q4_7 + q5_6 * q5_7 + q6_6 * q6_7 + q7_6 * q7_7 + q8_6 * q8_7 + q9_6 * q9_7 + q10_6 * q10_7) + (q1_7 * q1_7 * d1 + q2_7 * q2_7 * d2 + q3_7 * q3_7 * d3 + q4_7 * q4_7 * d4 + q5_7 * q5_7 * d5 + q6_7 * q6_7 * d6 + q7_7 * q7_7 * d7 + q8_7 * q8_7 * d8 + q9_7 * q9_7 * d9 + q10_7 * q10_7 * d10 - 2.0)*(q1_7 * q1_7 * d1 + q2_7 * q2_7 * d2 + q3_7 * q3_7 * d3 + q4_7 * q4_7 * d4 + q5_7 * q5_7 * d5 + q6_7 * q6_7 * d6 + q7_7 * q7_7 * d7 + q8_7 * q8_7 * d8 + q9_7 * q9_7 * d9 + q10_7 * q10_7 * d10 - 2.0) + (q1_7 * q1_7 + q2_7 * q2_7 + q3_7 * q3_7 + q4_7 * q4_7 + q5_7 * q5_7 + q6_7 * q6_7 + q7_7 * q7_7 + q8_7 * q8_7 + q9_7 * q9_7 + q10_7 * q10_7 - 1.0)*(q1_7 * q1_7 + q2_7 * q2_7 + q3_7 * q3_7 + q4_7 * q4_7 + q5_7 * q5_7 + q6_7 * q6_7 + q7_7 * q7_7 + q8_7 * q8_7 + q9_7 * q9_7 + q10_7 * q10_7 - 1.0) + (q1_1 * q1_8 * d1 + q2_1 * q2_8 * d2 + q3_1 * q3_8 * d3 + q4_1 * q4_8 * d4 + q5_1 * q5_8 * d5 + q6_1 * q6_8 * d6 + q7_1 * q7_8 * d7 + q8_1 * q8_8 * d8 + q9_1 * q9_8 * d9 + q10_1 * q10_8 * d10)*(q1_1 * q1_8 * d1 + q2_1 * q2_8 * d2 + q3_1 * q3_8 * d3 + q4_1 * q4_8 * d4 + q5_1 * q5_8 * d5 + q6_1 * q6_8 * d6 + q7_1 * q7_8 * d7 + q8_1 * q8_8 * d8 + q9_1 * q9_8 * d9 + q10_1 * q10_8 * d10) + (q1_1 * q1_8 + q2_1 * q2_8 + q3_1 * q3_8 + q4_1 * q4_8 + q5_1 * q5_8 + q6_1 * q6_8 + q7_1 * q7_8 + q8_1 * q8_8 + q9_1 * q9_8 + q10_1 * q10_8)*(q1_1 * q1_8 + q2_1 * q2_8 + q3_1 * q3_8 + q4_1 * q4_8 + q5_1 * q5_8 + q6_1 * q6_8 + q7_1 * q7_8 + q8_1 * q8_8 + q9_1 * q9_8 + q10_1 * q10_8) + (q1_2 * q1_8 * d1 + q2_2 * q2_8 * d2 + q3_2 * q3_8 * d3 + q4_2 * q4_8 * d4 + q5_2 * q5_8 * d5 + q6_2 * q6_8 * d6 + q7_2 * q7_8 * d7 + q8_2 * q8_8 * d8 + q9_2 * q9_8 * d9 + q10_2 * q10_8 * d10)*(q1_2 * q1_8 * d1 + q2_2 * q2_8 * d2 + q3_2 * q3_8 * d3 + q4_2 * q4_8 * d4 + q5_2 * q5_8 * d5 + q6_2 * q6_8 * d6 + q7_2 * q7_8 * d7 + q8_2 * q8_8 * d8 + q9_2 * q9_8 * d9 + q10_2 * q10_8 * d10) + (q1_2 * q1_8 + q2_2 * q2_8 + q3_2 * q3_8 + q4_2 * q4_8 + q5_2 * q5_8 + q6_2 * q6_8 + q7_2 * q7_8 + q8_2 * q8_8 + q9_2 * q9_8 + q10_2 * q10_8)*(q1_2 * q1_8 + q2_2 * q2_8 + q3_2 * q3_8 + q4_2 * q4_8 + q5_2 * q5_8 + q6_2 * q6_8 + q7_2 * q7_8 + q8_2 * q8_8 + q9_2 * q9_8 + q10_2 * q10_8) + (q1_3 * q1_8 * d1 + q2_3 * q2_8 * d2 + q3_3 * q3_8 * d3 + q4_3 * q4_8 * d4 + q5_3 * q5_8 * d5 + q6_3 * q6_8 * d6 + q7_3 * q7_8 * d7 + q8_3 * q8_8 * d8 + q9_3 * q9_8 * d9 + q10_3 * q10_8 * d10)*(q1_3 * q1_8 * d1 + q2_3 * q2_8 * d2 + q3_3 * q3_8 * d3 + q4_3 * q4_8 * d4 + q5_3 * q5_8 * d5 + q6_3 * q6_8 * d6 + q7_3 * q7_8 * d7 + q8_3 * q8_8 * d8 + q9_3 * q9_8 * d9 + q10_3 * q10_8 * d10) + (q1_3 * q1_8 + q2_3 * q2_8 + q3_3 * q3_8 + q4_3 * q4_8 + q5_3 * q5_8 + q6_3 * q6_8 + q7_3 * q7_8 + q8_3 * q8_8 + q9_3 * q9_8 + q10_3 * q10_8)*(q1_3 * q1_8 + q2_3 * q2_8 + q3_3 * q3_8 + q4_3 * q4_8 + q5_3 * q5_8 + q6_3 * q6_8 + q7_3 * q7_8 + q8_3 * q8_8 + q9_3 * q9_8 + q10_3 * q10_8) + (q1_4 * q1_8 * d1 + q2_4 * q2_8 * d2 + q3_4 * q3_8 * d3 + q4_4 * q4_8 * d4 + q5_4 * q5_8 * d5 + q6_4 * q6_8 * d6 + q7_4 * q7_8 * d7 + q8_4 * q8_8 * d8 + q9_4 * q9_8 * d9 + q10_4 * q10_8 * d10)*(q1_4 * q1_8 * d1 + q2_4 * q2_8 * d2 + q3_4 * q3_8 * d3 + q4_4 * q4_8 * d4 + q5_4 * q5_8 * d5 + q6_4 * q6_8 * d6 + q7_4 * q7_8 * d7 + q8_4 * q8_8 * d8 + q9_4 * q9_8 * d9 + q10_4 * q10_8 * d10) + (q1_4 * q1_8 + q2_4 * q2_8 + q3_4 * q3_8 + q4_4 * q4_8 + q5_4 * q5_8 + q6_4 * q6_8 + q7_4 * q7_8 + q8_4 * q8_8 + q9_4 * q9_8 + q10_4 * q10_8)*(q1_4 * q1_8 + q2_4 * q2_8 + q3_4 * q3_8 + q4_4 * q4_8 + q5_4 * q5_8 + q6_4 * q6_8 + q7_4 * q7_8 + q8_4 * q8_8 + q9_4 * q9_8 + q10_4 * q10_8) + (q1_5 * q1_8 * d1 + q2_5 * q2_8 * d2 + q3_5 * q3_8 * d3 + q4_5 * q4_8 * d4 + q5_5 * q5_8 * d5 + q6_5 * q6_8 * d6 + q7_5 * q7_8 * d7 + q8_5 * q8_8 * d8 + q9_5 * q9_8 * d9 + q10_5 * q10_8 * d10)*(q1_5 * q1_8 * d1 + q2_5 * q2_8 * d2 + q3_5 * q3_8 * d3 + q4_5 * q4_8 * d4 + q5_5 * q5_8 * d5 + q6_5 * q6_8 * d6 + q7_5 * q7_8 * d7 + q8_5 * q8_8 * d8 + q9_5 * q9_8 * d9 + q10_5 * q10_8 * d10) + (q1_5 * q1_8 + q2_5 * q2_8 + q3_5 * q3_8 + q4_5 * q4_8 + q5_5 * q5_8 + q6_5 * q6_8 + q7_5 * q7_8 + q8_5 * q8_8 + q9_5 * q9_8 + q10_5 * q10_8)*(q1_5 * q1_8 + q2_5 * q2_8 + q3_5 * q3_8 + q4_5 * q4_8 + q5_5 * q5_8 + q6_5 * q6_8 + q7_5 * q7_8 + q8_5 * q8_8 + q9_5 * q9_8 + q10_5 * q10_8) + (q1_6 * q1_8 * d1 + q2_6 * q2_8 * d2 + q3_6 * q3_8 * d3 + q4_6 * q4_8 * d4 + q5_6 * q5_8 * d5 + q6_6 * q6_8 * d6 + q7_6 * q7_8 * d7 + q8_6 * q8_8 * d8 + q9_6 * q9_8 * d9 + q10_6 * q10_8 * d10)*(q1_6 * q1_8 * d1 + q2_6 * q2_8 * d2 + q3_6 * q3_8 * d3 + q4_6 * q4_8 * d4 + q5_6 * q5_8 * d5 + q6_6 * q6_8 * d6 + q7_6 * q7_8 * d7 + q8_6 * q8_8 * d8 + q9_6 * q9_8 * d9 + q10_6 * q10_8 * d10) + (q1_6 * q1_8 + q2_6 * q2_8 + q3_6 * q3_8 + q4_6 * q4_8 + q5_6 * q5_8 + q6_6 * q6_8 + q7_6 * q7_8 + q8_6 * q8_8 + q9_6 * q9_8 + q10_6 * q10_8)*(q1_6 * q1_8 + q2_6 * q2_8 + q3_6 * q3_8 + q4_6 * q4_8 + q5_6 * q5_8 + q6_6 * q6_8 + q7_6 * q7_8 + q8_6 * q8_8 + q9_6 * q9_8 + q10_6 * q10_8) + (q1_7 * q1_8 * d1 + q2_7 * q2_8 * d2 + q3_7 * q3_8 * d3 + q4_7 * q4_8 * d4 + q5_7 * q5_8 * d5 + q6_7 * q6_8 * d6 + q7_7 * q7_8 * d7 + q8_7 * q8_8 * d8 + q9_7 * q9_8 * d9 + q10_7 * q10_8 * d10 + 1.0)*(q1_7 * q1_8 * d1 + q2_7 * q2_8 * d2 + q3_7 * q3_8 * d3 + q4_7 * q4_8 * d4 + q5_7 * q5_8 * d5 + q6_7 * q6_8 * d6 + q7_7 * q7_8 * d7 + q8_7 * q8_8 * d8 + q9_7 * q9_8 * d9 + q10_7 * q10_8 * d10 + 1.0) + (q1_7 * q1_8 + q2_7 * q2_8 + q3_7 * q3_8 + q4_7 * q4_8 + q5_7 * q5_8 + q6_7 * q6_8 + q7_7 * q7_8 + q8_7 * q8_8 + q9_7 * q9_8 + q10_7 * q10_8)*(q1_7 * q1_8 + q2_7 * q2_8 + q3_7 * q3_8 + q4_7 * q4_8 + q5_7 * q5_8 + q6_7 * q6_8 + q7_7 * q7_8 + q8_7 * q8_8 + q9_7 * q9_8 + q10_7 * q10_8) + (q1_8 * q1_8 * d1 + q2_8 * q2_8 * d2 + q3_8 * q3_8 * d3 + q4_8 * q4_8 * d4 + q5_8 * q5_8 * d5 + q6_8 * q6_8 * d6 + q7_8 * q7_8 * d7 + q8_8 * q8_8 * d8 + q9_8 * q9_8 * d9 + q10_8 * q10_8 * d10 - 2.0)*(q1_8 * q1_8 * d1 + q2_8 * q2_8 * d2 + q3_8 * q3_8 * d3 + q4_8 * q4_8 * d4 + q5_8 * q5_8 * d5 + q6_8 * q6_8 * d6 + q7_8 * q7_8 * d7 + q8_8 * q8_8 * d8 + q9_8 * q9_8 * d9 + q10_8 * q10_8 * d10 - 2.0) + (q1_8 * q1_8 + q2_8 * q2_8 + q3_8 * q3_8 + q4_8 * q4_8 + q5_8 * q5_8 + q6_8 * q6_8 + q7_8 * q7_8 + q8_8 * q8_8 + q9_8 * q9_8 + q10_8 * q10_8 - 1.0)*(q1_8 * q1_8 + q2_8 * q2_8 + q3_8 * q3_8 + q4_8 * q4_8 + q5_8 * q5_8 + q6_8 * q6_8 + q7_8 * q7_8 + q8_8 * q8_8 + q9_8 * q9_8 + q10_8 * q10_8 - 1.0) + (q1_1 * q1_9 * d1 + q2_1 * q2_9 * d2 + q3_1 * q3_9 * d3 + q4_1 * q4_9 * d4 + q5_1 * q5_9 * d5 + q6_1 * q6_9 * d6 + q7_1 * q7_9 * d7 + q8_1 * q8_9 * d8 + q9_1 * q9_9 * d9 + q10_1 * q10_9 * d10)*(q1_1 * q1_9 * d1 + q2_1 * q2_9 * d2 + q3_1 * q3_9 * d3 + q4_1 * q4_9 * d4 + q5_1 * q5_9 * d5 + q6_1 * q6_9 * d6 + q7_1 * q7_9 * d7 + q8_1 * q8_9 * d8 + q9_1 * q9_9 * d9 + q10_1 * q10_9 * d10) + (q1_1 * q1_9 + q2_1 * q2_9 + q3_1 * q3_9 + q4_1 * q4_9 + q5_1 * q5_9 + q6_1 * q6_9 + q7_1 * q7_9 + q8_1 * q8_9 + q9_1 * q9_9 + q10_1 * q10_9)*(q1_1 * q1_9 + q2_1 * q2_9 + q3_1 * q3_9 + q4_1 * q4_9 + q5_1 * q5_9 + q6_1 * q6_9 + q7_1 * q7_9 + q8_1 * q8_9 + q9_1 * q9_9 + q10_1 * q10_9) + (q1_2 * q1_9 * d1 + q2_2 * q2_9 * d2 + q3_2 * q3_9 * d3 + q4_2 * q4_9 * d4 + q5_2 * q5_9 * d5 + q6_2 * q6_9 * d6 + q7_2 * q7_9 * d7 + q8_2 * q8_9 * d8 + q9_2 * q9_9 * d9 + q10_2 * q10_9 * d10)*(q1_2 * q1_9 * d1 + q2_2 * q2_9 * d2 + q3_2 * q3_9 * d3 + q4_2 * q4_9 * d4 + q5_2 * q5_9 * d5 + q6_2 * q6_9 * d6 + q7_2 * q7_9 * d7 + q8_2 * q8_9 * d8 + q9_2 * q9_9 * d9 + q10_2 * q10_9 * d10) + (q1_2 * q1_9 + q2_2 * q2_9 + q3_2 * q3_9 + q4_2 * q4_9 + q5_2 * q5_9 + q6_2 * q6_9 + q7_2 * q7_9 + q8_2 * q8_9 + q9_2 * q9_9 + q10_2 * q10_9)*(q1_2 * q1_9 + q2_2 * q2_9 + q3_2 * q3_9 + q4_2 * q4_9 + q5_2 * q5_9 + q6_2 * q6_9 + q7_2 * q7_9 + q8_2 * q8_9 + q9_2 * q9_9 + q10_2 * q10_9) + (q1_3 * q1_9 * d1 + q2_3 * q2_9 * d2 + q3_3 * q3_9 * d3 + q4_3 * q4_9 * d4 + q5_3 * q5_9 * d5 + q6_3 * q6_9 * d6 + q7_3 * q7_9 * d7 + q8_3 * q8_9 * d8 + q9_3 * q9_9 * d9 + q10_3 * q10_9 * d10)*(q1_3 * q1_9 * d1 + q2_3 * q2_9 * d2 + q3_3 * q3_9 * d3 + q4_3 * q4_9 * d4 + q5_3 * q5_9 * d5 + q6_3 * q6_9 * d6 + q7_3 * q7_9 * d7 + q8_3 * q8_9 * d8 + q9_3 * q9_9 * d9 + q10_3 * q10_9 * d10) + (q1_3 * q1_9 + q2_3 * q2_9 + q3_3 * q3_9 + q4_3 * q4_9 + q5_3 * q5_9 + q6_3 * q6_9 + q7_3 * q7_9 + q8_3 * q8_9 + q9_3 * q9_9 + q10_3 * q10_9)*(q1_3 * q1_9 + q2_3 * q2_9 + q3_3 * q3_9 + q4_3 * q4_9 + q5_3 * q5_9 + q6_3 * q6_9 + q7_3 * q7_9 + q8_3 * q8_9 + q9_3 * q9_9 + q10_3 * q10_9) + (q1_4 * q1_9 * d1 + q2_4 * q2_9 * d2 + q3_4 * q3_9 * d3 + q4_4 * q4_9 * d4 + q5_4 * q5_9 * d5 + q6_4 * q6_9 * d6 + q7_4 * q7_9 * d7 + q8_4 * q8_9 * d8 + q9_4 * q9_9 * d9 + q10_4 * q10_9 * d10)*(q1_4 * q1_9 * d1 + q2_4 * q2_9 * d2 + q3_4 * q3_9 * d3 + q4_4 * q4_9 * d4 + q5_4 * q5_9 * d5 + q6_4 * q6_9 * d6 + q7_4 * q7_9 * d7 + q8_4 * q8_9 * d8 + q9_4 * q9_9 * d9 + q10_4 * q10_9 * d10) + (q1_4 * q1_9 + q2_4 * q2_9 + q3_4 * q3_9 + q4_4 * q4_9 + q5_4 * q5_9 + q6_4 * q6_9 + q7_4 * q7_9 + q8_4 * q8_9 + q9_4 * q9_9 + q10_4 * q10_9)*(q1_4 * q1_9 + q2_4 * q2_9 + q3_4 * q3_9 + q4_4 * q4_9 + q5_4 * q5_9 + q6_4 * q6_9 + q7_4 * q7_9 + q8_4 * q8_9 + q9_4 * q9_9 + q10_4 * q10_9) + (q1_5 * q1_9 * d1 + q2_5 * q2_9 * d2 + q3_5 * q3_9 * d3 + q4_5 * q4_9 * d4 + q5_5 * q5_9 * d5 + q6_5 * q6_9 * d6 + q7_5 * q7_9 * d7 + q8_5 * q8_9 * d8 + q9_5 * q9_9 * d9 + q10_5 * q10_9 * d10)*(q1_5 * q1_9 * d1 + q2_5 * q2_9 * d2 + q3_5 * q3_9 * d3 + q4_5 * q4_9 * d4 + q5_5 * q5_9 * d5 + q6_5 * q6_9 * d6 + q7_5 * q7_9 * d7 + q8_5 * q8_9 * d8 + q9_5 * q9_9 * d9 + q10_5 * q10_9 * d10) + (q1_5 * q1_9 + q2_5 * q2_9 + q3_5 * q3_9 + q4_5 * q4_9 + q5_5 * q5_9 + q6_5 * q6_9 + q7_5 * q7_9 + q8_5 * q8_9 + q9_5 * q9_9 + q10_5 * q10_9)*(q1_5 * q1_9 + q2_5 * q2_9 + q3_5 * q3_9 + q4_5 * q4_9 + q5_5 * q5_9 + q6_5 * q6_9 + q7_5 * q7_9 + q8_5 * q8_9 + q9_5 * q9_9 + q10_5 * q10_9) + (q1_6 * q1_9 * d1 + q2_6 * q2_9 * d2 + q3_6 * q3_9 * d3 + q4_6 * q4_9 * d4 + q5_6 * q5_9 * d5 + q6_6 * q6_9 * d6 + q7_6 * q7_9 * d7 + q8_6 * q8_9 * d8 + q9_6 * q9_9 * d9 + q10_6 * q10_9 * d10)*(q1_6 * q1_9 * d1 + q2_6 * q2_9 * d2 + q3_6 * q3_9 * d3 + q4_6 * q4_9 * d4 + q5_6 * q5_9 * d5 + q6_6 * q6_9 * d6 + q7_6 * q7_9 * d7 + q8_6 * q8_9 * d8 + q9_6 * q9_9 * d9 + q10_6 * q10_9 * d10) + (q1_6 * q1_9 + q2_6 * q2_9 + q3_6 * q3_9 + q4_6 * q4_9 + q5_6 * q5_9 + q6_6 * q6_9 + q7_6 * q7_9 + q8_6 * q8_9 + q9_6 * q9_9 + q10_6 * q10_9)*(q1_6 * q1_9 + q2_6 * q2_9 + q3_6 * q3_9 + q4_6 * q4_9 + q5_6 * q5_9 + q6_6 * q6_9 + q7_6 * q7_9 + q8_6 * q8_9 + q9_6 * q9_9 + q10_6 * q10_9) + (q1_7 * q1_9 * d1 + q2_7 * q2_9 * d2 + q3_7 * q3_9 * d3 + q4_7 * q4_9 * d4 + q5_7 * q5_9 * d5 + q6_7 * q6_9 * d6 + q7_7 * q7_9 * d7 + q8_7 * q8_9 * d8 + q9_7 * q9_9 * d9 + q10_7 * q10_9 * d10)*(q1_7 * q1_9 * d1 + q2_7 * q2_9 * d2 + q3_7 * q3_9 * d3 + q4_7 * q4_9 * d4 + q5_7 * q5_9 * d5 + q6_7 * q6_9 * d6 + q7_7 * q7_9 * d7 + q8_7 * q8_9 * d8 + q9_7 * q9_9 * d9 + q10_7 * q10_9 * d10) + (q1_7 * q1_9 + q2_7 * q2_9 + q3_7 * q3_9 + q4_7 * q4_9 + q5_7 * q5_9 + q6_7 * q6_9 + q7_7 * q7_9 + q8_7 * q8_9 + q9_7 * q9_9 + q10_7 * q10_9)*(q1_7 * q1_9 + q2_7 * q2_9 + q3_7 * q3_9 + q4_7 * q4_9 + q5_7 * q5_9 + q6_7 * q6_9 + q7_7 * q7_9 + q8_7 * q8_9 + q9_7 * q9_9 + q10_7 * q10_9) + (q1_8 * q1_9 * d1 + q2_8 * q2_9 * d2 + q3_8 * q3_9 * d3 + q4_8 * q4_9 * d4 + q5_8 * q5_9 * d5 + q6_8 * q6_9 * d6 + q7_8 * q7_9 * d7 + q8_8 * q8_9 * d8 + q9_8 * q9_9 * d9 + q10_8 * q10_9 * d10 + 1.0)*(q1_8 * q1_9 * d1 + q2_8 * q2_9 * d2 + q3_8 * q3_9 * d3 + q4_8 * q4_9 * d4 + q5_8 * q5_9 * d5 + q6_8 * q6_9 * d6 + q7_8 * q7_9 * d7 + q8_8 * q8_9 * d8 + q9_8 * q9_9 * d9 + q10_8 * q10_9 * d10 + 1.0) + (q1_8 * q1_9 + q2_8 * q2_9 + q3_8 * q3_9 + q4_8 * q4_9 + q5_8 * q5_9 + q6_8 * q6_9 + q7_8 * q7_9 + q8_8 * q8_9 + q9_8 * q9_9 + q10_8 * q10_9)*(q1_8 * q1_9 + q2_8 * q2_9 + q3_8 * q3_9 + q4_8 * q4_9 + q5_8 * q5_9 + q6_8 * q6_9 + q7_8 * q7_9 + q8_8 * q8_9 + q9_8 * q9_9 + q10_8 * q10_9) + (q1_9 * q1_9 * d1 + q2_9 * q2_9 * d2 + q3_9 * q3_9 * d3 + q4_9 * q4_9 * d4 + q5_9 * q5_9 * d5 + q6_9 * q6_9 * d6 + q7_9 * q7_9 * d7 + q8_9 * q8_9 * d8 + q9_9 * q9_9 * d9 + q10_9 * q10_9 * d10 - 2.0)*(q1_9 * q1_9 * d1 + q2_9 * q2_9 * d2 + q3_9 * q3_9 * d3 + q4_9 * q4_9 * d4 + q5_9 * q5_9 * d5 + q6_9 * q6_9 * d6 + q7_9 * q7_9 * d7 + q8_9 * q8_9 * d8 + q9_9 * q9_9 * d9 + q10_9 * q10_9 * d10 - 2.0) + (q1_9 * q1_9 + q2_9 * q2_9 + q3_9 * q3_9 + q4_9 * q4_9 + q5_9 * q5_9 + q6_9 * q6_9 + q7_9 * q7_9 + q8_9 * q8_9 + q9_9 * q9_9 + q10_9 * q10_9 - 1.0)*(q1_9 * q1_9 + q2_9 * q2_9 + q3_9 * q3_9 + q4_9 * q4_9 + q5_9 * q5_9 + q6_9 * q6_9 + q7_9 * q7_9 + q8_9 * q8_9 + q9_9 * q9_9 + q10_9 * q10_9 - 1.0) + (q1_1 * q1_10 * d1 + q2_1 * q2_10 * d2 + q3_1 * q3_10 * d3 + q4_1 * q4_10 * d4 + q5_1 * q5_10 * d5 + q6_1 * q6_10 * d6 + q7_1 * q7_10 * d7 + q8_1 * q8_10 * d8 + q9_1 * q9_10 * d9 + q10_1 * q10_10 * d10)*(q1_1 * q1_10 * d1 + q2_1 * q2_10 * d2 + q3_1 * q3_10 * d3 + q4_1 * q4_10 * d4 + q5_1 * q5_10 * d5 + q6_1 * q6_10 * d6 + q7_1 * q7_10 * d7 + q8_1 * q8_10 * d8 + q9_1 * q9_10 * d9 + q10_1 * q10_10 * d10) + (q1_1 * q1_10 + q2_1 * q2_10 + q3_1 * q3_10 + q4_1 * q4_10 + q5_1 * q5_10 + q6_1 * q6_10 + q7_1 * q7_10 + q8_1 * q8_10 + q9_1 * q9_10 + q10_1 * q10_10)*(q1_1 * q1_10 + q2_1 * q2_10 + q3_1 * q3_10 + q4_1 * q4_10 + q5_1 * q5_10 + q6_1 * q6_10 + q7_1 * q7_10 + q8_1 * q8_10 + q9_1 * q9_10 + q10_1 * q10_10) + (q1_2 * q1_10 * d1 + q2_2 * q2_10 * d2 + q3_2 * q3_10 * d3 + q4_2 * q4_10 * d4 + q5_2 * q5_10 * d5 + q6_2 * q6_10 * d6 + q7_2 * q7_10 * d7 + q8_2 * q8_10 * d8 + q9_2 * q9_10 * d9 + q10_2 * q10_10 * d10)*(q1_2 * q1_10 * d1 + q2_2 * q2_10 * d2 + q3_2 * q3_10 * d3 + q4_2 * q4_10 * d4 + q5_2 * q5_10 * d5 + q6_2 * q6_10 * d6 + q7_2 * q7_10 * d7 + q8_2 * q8_10 * d8 + q9_2 * q9_10 * d9 + q10_2 * q10_10 * d10) + (q1_2 * q1_10 + q2_2 * q2_10 + q3_2 * q3_10 + q4_2 * q4_10 + q5_2 * q5_10 + q6_2 * q6_10 + q7_2 * q7_10 + q8_2 * q8_10 + q9_2 * q9_10 + q10_2 * q10_10)*(q1_2 * q1_10 + q2_2 * q2_10 + q3_2 * q3_10 + q4_2 * q4_10 + q5_2 * q5_10 + q6_2 * q6_10 + q7_2 * q7_10 + q8_2 * q8_10 + q9_2 * q9_10 + q10_2 * q10_10) + (q1_3 * q1_10 * d1 + q2_3 * q2_10 * d2 + q3_3 * q3_10 * d3 + q4_3 * q4_10 * d4 + q5_3 * q5_10 * d5 + q6_3 * q6_10 * d6 + q7_3 * q7_10 * d7 + q8_3 * q8_10 * d8 + q9_3 * q9_10 * d9 + q10_3 * q10_10 * d10)*(q1_3 * q1_10 * d1 + q2_3 * q2_10 * d2 + q3_3 * q3_10 * d3 + q4_3 * q4_10 * d4 + q5_3 * q5_10 * d5 + q6_3 * q6_10 * d6 + q7_3 * q7_10 * d7 + q8_3 * q8_10 * d8 + q9_3 * q9_10 * d9 + q10_3 * q10_10 * d10) + (q1_3 * q1_10 + q2_3 * q2_10 + q3_3 * q3_10 + q4_3 * q4_10 + q5_3 * q5_10 + q6_3 * q6_10 + q7_3 * q7_10 + q8_3 * q8_10 + q9_3 * q9_10 + q10_3 * q10_10)*(q1_3 * q1_10 + q2_3 * q2_10 + q3_3 * q3_10 + q4_3 * q4_10 + q5_3 * q5_10 + q6_3 * q6_10 + q7_3 * q7_10 + q8_3 * q8_10 + q9_3 * q9_10 + q10_3 * q10_10) + (q1_4 * q1_10 * d1 + q2_4 * q2_10 * d2 + q3_4 * q3_10 * d3 + q4_4 * q4_10 * d4 + q5_4 * q5_10 * d5 + q6_4 * q6_10 * d6 + q7_4 * q7_10 * d7 + q8_4 * q8_10 * d8 + q9_4 * q9_10 * d9 + q10_4 * q10_10 * d10)*(q1_4 * q1_10 * d1 + q2_4 * q2_10 * d2 + q3_4 * q3_10 * d3 + q4_4 * q4_10 * d4 + q5_4 * q5_10 * d5 + q6_4 * q6_10 * d6 + q7_4 * q7_10 * d7 + q8_4 * q8_10 * d8 + q9_4 * q9_10 * d9 + q10_4 * q10_10 * d10) + (q1_4 * q1_10 + q2_4 * q2_10 + q3_4 * q3_10 + q4_4 * q4_10 + q5_4 * q5_10 + q6_4 * q6_10 + q7_4 * q7_10 + q8_4 * q8_10 + q9_4 * q9_10 + q10_4 * q10_10)*(q1_4 * q1_10 + q2_4 * q2_10 + q3_4 * q3_10 + q4_4 * q4_10 + q5_4 * q5_10 + q6_4 * q6_10 + q7_4 * q7_10 + q8_4 * q8_10 + q9_4 * q9_10 + q10_4 * q10_10) + (q1_5 * q1_10 * d1 + q2_5 * q2_10 * d2 + q3_5 * q3_10 * d3 + q4_5 * q4_10 * d4 + q5_5 * q5_10 * d5 + q6_5 * q6_10 * d6 + q7_5 * q7_10 * d7 + q8_5 * q8_10 * d8 + q9_5 * q9_10 * d9 + q10_5 * q10_10 * d10)*(q1_5 * q1_10 * d1 + q2_5 * q2_10 * d2 + q3_5 * q3_10 * d3 + q4_5 * q4_10 * d4 + q5_5 * q5_10 * d5 + q6_5 * q6_10 * d6 + q7_5 * q7_10 * d7 + q8_5 * q8_10 * d8 + q9_5 * q9_10 * d9 + q10_5 * q10_10 * d10) + (q1_5 * q1_10 + q2_5 * q2_10 + q3_5 * q3_10 + q4_5 * q4_10 + q5_5 * q5_10 + q6_5 * q6_10 + q7_5 * q7_10 + q8_5 * q8_10 + q9_5 * q9_10 + q10_5 * q10_10)*(q1_5 * q1_10 + q2_5 * q2_10 + q3_5 * q3_10 + q4_5 * q4_10 + q5_5 * q5_10 + q6_5 * q6_10 + q7_5 * q7_10 + q8_5 * q8_10 + q9_5 * q9_10 + q10_5 * q10_10) + (q1_6 * q1_10 * d1 + q2_6 * q2_10 * d2 + q3_6 * q3_10 * d3 + q4_6 * q4_10 * d4 + q5_6 * q5_10 * d5 + q6_6 * q6_10 * d6 + q7_6 * q7_10 * d7 + q8_6 * q8_10 * d8 + q9_6 * q9_10 * d9 + q10_6 * q10_10 * d10)*(q1_6 * q1_10 * d1 + q2_6 * q2_10 * d2 + q3_6 * q3_10 * d3 + q4_6 * q4_10 * d4 + q5_6 * q5_10 * d5 + q6_6 * q6_10 * d6 + q7_6 * q7_10 * d7 + q8_6 * q8_10 * d8 + q9_6 * q9_10 * d9 + q10_6 * q10_10 * d10) + (q1_6 * q1_10 + q2_6 * q2_10 + q3_6 * q3_10 + q4_6 * q4_10 + q5_6 * q5_10 + q6_6 * q6_10 + q7_6 * q7_10 + q8_6 * q8_10 + q9_6 * q9_10 + q10_6 * q10_10)*(q1_6 * q1_10 + q2_6 * q2_10 + q3_6 * q3_10 + q4_6 * q4_10 + q5_6 * q5_10 + q6_6 * q6_10 + q7_6 * q7_10 + q8_6 * q8_10 + q9_6 * q9_10 + q10_6 * q10_10) + (q1_7 * q1_10 * d1 + q2_7 * q2_10 * d2 + q3_7 * q3_10 * d3 + q4_7 * q4_10 * d4 + q5_7 * q5_10 * d5 + q6_7 * q6_10 * d6 + q7_7 * q7_10 * d7 + q8_7 * q8_10 * d8 + q9_7 * q9_10 * d9 + q10_7 * q10_10 * d10)*(q1_7 * q1_10 * d1 + q2_7 * q2_10 * d2 + q3_7 * q3_10 * d3 + q4_7 * q4_10 * d4 + q5_7 * q5_10 * d5 + q6_7 * q6_10 * d6 + q7_7 * q7_10 * d7 + q8_7 * q8_10 * d8 + q9_7 * q9_10 * d9 + q10_7 * q10_10 * d10) + (q1_7 * q1_10 + q2_7 * q2_10 + q3_7 * q3_10 + q4_7 * q4_10 + q5_7 * q5_10 + q6_7 * q6_10 + q7_7 * q7_10 + q8_7 * q8_10 + q9_7 * q9_10 + q10_7 * q10_10)*(q1_7 * q1_10 + q2_7 * q2_10 + q3_7 * q3_10 + q4_7 * q4_10 + q5_7 * q5_10 + q6_7 * q6_10 + q7_7 * q7_10 + q8_7 * q8_10 + q9_7 * q9_10 + q10_7 * q10_10) + (q1_8 * q1_10 * d1 + q2_8 * q2_10 * d2 + q3_8 * q3_10 * d3 + q4_8 * q4_10 * d4 + q5_8 * q5_10 * d5 + q6_8 * q6_10 * d6 + q7_8 * q7_10 * d7 + q8_8 * q8_10 * d8 + q9_8 * q9_10 * d9 + q10_8 * q10_10 * d10)*(q1_8 * q1_10 * d1 + q2_8 * q2_10 * d2 + q3_8 * q3_10 * d3 + q4_8 * q4_10 * d4 + q5_8 * q5_10 * d5 + q6_8 * q6_10 * d6 + q7_8 * q7_10 * d7 + q8_8 * q8_10 * d8 + q9_8 * q9_10 * d9 + q10_8 * q10_10 * d10) + (q1_8 * q1_10 + q2_8 * q2_10 + q3_8 * q3_10 + q4_8 * q4_10 + q5_8 * q5_10 + q6_8 * q6_10 + q7_8 * q7_10 + q8_8 * q8_10 + q9_8 * q9_10 + q10_8 * q10_10)*(q1_8 * q1_10 + q2_8 * q2_10 + q3_8 * q3_10 + q4_8 * q4_10 + q5_8 * q5_10 + q6_8 * q6_10 + q7_8 * q7_10 + q8_8 * q8_10 + q9_8 * q9_10 + q10_8 * q10_10) + (q1_9 * q1_10 * d1 + q2_9 * q2_10 * d2 + q3_9 * q3_10 * d3 + q4_9 * q4_10 * d4 + q5_9 * q5_10 * d5 + q6_9 * q6_10 * d6 + q7_9 * q7_10 * d7 + q8_9 * q8_10 * d8 + q9_9 * q9_10 * d9 + q10_9 * q10_10 * d10 + 1.0)*(q1_9 * q1_10 * d1 + q2_9 * q2_10 * d2 + q3_9 * q3_10 * d3 + q4_9 * q4_10 * d4 + q5_9 * q5_10 * d5 + q6_9 * q6_10 * d6 + q7_9 * q7_10 * d7 + q8_9 * q8_10 * d8 + q9_9 * q9_10 * d9 + q10_9 * q10_10 * d10 + 1.0) + (q1_9 * q1_10 + q2_9 * q2_10 + q3_9 * q3_10 + q4_9 * q4_10 + q5_9 * q5_10 + q6_9 * q6_10 + q7_9 * q7_10 + q8_9 * q8_10 + q9_9 * q9_10 + q10_9 * q10_10)*(q1_9 * q1_10 + q2_9 * q2_10 + q3_9 * q3_10 + q4_9 * q4_10 + q5_9 * q5_10 + q6_9 * q6_10 + q7_9 * q7_10 + q8_9 * q8_10 + q9_9 * q9_10 + q10_9 * q10_10) + (q1_10 * q1_10 * d1 + q2_10 * q2_10 * d2 + q3_10 * q3_10 * d3 + q4_10 * q4_10 * d4 + q5_10 * q5_10 * d5 + q6_10 * q6_10 * d6 + q7_10 * q7_10 * d7 + q8_10 * q8_10 * d8 + q9_10 * q9_10 * d9 + q10_10 * q10_10 * d10 - 2.0)*(q1_10 * q1_10 * d1 + q2_10 * q2_10 * d2 + q3_10 * q3_10 * d3 + q4_10 * q4_10 * d4 + q5_10 * q5_10 * d5 + q6_10 * q6_10 * d6 + q7_10 * q7_10 * d7 + q8_10 * q8_10 * d8 + q9_10 * q9_10 * d9 + q10_10 * q10_10 * d10 - 2.0) + (q1_10 * q1_10 + q2_10 * q2_10 + q3_10 * q3_10 + q4_10 * q4_10 + q5_10 * q5_10 + q6_10 * q6_10 + q7_10 * q7_10 + q8_10 * q8_10 + q9_10 * q9_10 + q10_10 * q10_10 - 1.0)*(q1_10 * q1_10 + q2_10 * q2_10 + q3_10 * q3_10 + q4_10 * q4_10 + q5_10 * q5_10 + q6_10 * q6_10 + q7_10 * q7_10 + q8_10 * q8_10 + q9_10 * q9_10 + q10_10 * q10_10 - 1.0); solve; display d1; display q1_1; display q2_1; display q3_1; display q4_1; display q5_1; display q6_1; display q7_1; display q8_1; display q9_1; display q10_1; display d2; display q1_2; display q2_2; display q3_2; display q4_2; display q5_2; display q6_2; display q7_2; display q8_2; display q9_2; display q10_2; display d3; display q1_3; display q2_3; display q3_3; display q4_3; display q5_3; display q6_3; display q7_3; display q8_3; display q9_3; display q10_3; display d4; display q1_4; display q2_4; display q3_4; display q4_4; display q5_4; display q6_4; display q7_4; display q8_4; display q9_4; display q10_4; display d5; display q1_5; display q2_5; display q3_5; display q4_5; display q5_5; display q6_5; display q7_5; display q8_5; display q9_5; display q10_5; display d6; display q1_6; display q2_6; display q3_6; display q4_6; display q5_6; display q6_6; display q7_6; display q8_6; display q9_6; display q10_6; display d7; display q1_7; display q2_7; display q3_7; display q4_7; display q5_7; display q6_7; display q7_7; display q8_7; display q9_7; display q10_7; display d8; display q1_8; display q2_8; display q3_8; display q4_8; display q5_8; display q6_8; display q7_8; display q8_8; display q9_8; display q10_8; display d9; display q1_9; display q2_9; display q3_9; display q4_9; display q5_9; display q6_9; display q7_9; display q8_9; display q9_9; display q10_9; display d10; display q1_10; display q2_10; display q3_10; display q4_10; display q5_10; display q6_10; display q7_10; display q8_10; display q9_10; display q10_10; display obj;