q1, q2, q3, q4, q5, q6, q7, q8 = symbols('q1 q2 q3 q4 q5 q6 q7 q8 ', real=True)
f1 = q1 - qb - q2 - q3
f2 = q7 - qc + q2 + q4 -q5
f3 = q8 - qd - q4 - q6
f4 = q5 + q3 + q6 - qf
f5 = k2*q2**2 + k5*q5**2 + k3*q3**2
f6 = k4*q4**2 + k5*q5**2 -k6*q6**2
f7 = -30 + k1*q1**2 + k2*q2**2 - k7*q7**2
f8 = -40 + k7*q7**2 - k4*q4**2 - k8*q8**2
f_eqs = [f1,f2,f3,f4,f5,f6,f7,f8]