import classes

atome1 = classes.Atome('C', 1, [2, 3, 4, 5])
atome2 = classes.Atome('H', 2, [1])
atome3 = classes.Atome('H', 3, [1])
atome4 = classes.Atome('H', 4, [1])
atome5 = classes.Atome('C', 5, [1, 6, 6, 7])
atome6 = classes.Atome('C', 6, [5, 5, 8, 9])
atome7 = classes.Atome('H', 7, [5])
atome8 = classes.Atome('H', 8, [6])
atome9 = classes.Atome('C', 9, [6, 10, 11, 12])
atome10 = classes.Atome('H', 10, [9])
atome11 = classes.Atome('H', 11, [9])
atome12 = classes.Atome('H', 12, [9])

tableau_coeffs_HO = classes.Tableau_de_coefficients([atome1, atome2, atome3, atome4, atome5, atome6, atome7, atome8, atome9, atome10, atome11, atome12], [-0.158227, 0, -0.1922, 0.1922, 0.54271, 0.54271, 0, 0, -0.15828, 0, 0.15028, -0.15028])

tableau_coeffs_BV = classes.Tableau_de_coefficients([atome1, atome2, atome3, atome4, atome5, atome6, atome7, atome8, atome9, atome10, atome11, atome12], [0.04378, 0, -0.15027, 0.15027, 0.58242, -0.58242, 0, 0, -0.04379, 0, 0.15028, -0.15028])

HO = classes.Orbitale_moleculaire(-9.569977, tableau_coeffs_HO, 'HO')

BV = classes.Orbitale_moleculaire(1.276215, tableau_coeffs_BV, 'BV')

molecule = classes.Molecule('(Z)-but-2-ène', HO, BV)
