How to do this using Python?

Problem #3 The structure of a carbon dioxide molecule may be idealized as three masses connected by two springs, where the masses are the carbon and oxygen atoms, and the springs represent the chemical bonds between the atoms. The equations of motion for the three atoms, i.e. masses, are given by dx1 = -631 mo 212 = -kx + kx2, at2 dt2 mc-4122 = kx1 2kx2 + ku3, most = kx2 k23. a) By substituting the oscillatory solutions xi(t) = Ajeiwt (or xi(t) = Al cos(wt)], 22(t) = Azeiwt for x2(t) = A, cos(wt)], x3(t) = Azeiwt (or x3(t) = A3 cos(wt)], into the differential equations by hand, determine the system of algebraic equations for the natural frequencies of vibration w of the system of atoms from these equations of motion. b) For k = 1,420 kg/sa, mo = 2.6568 x 10-26 kg, and mc = 1.9926 x 10-26 kg, use built-in functions in Python, Matlab, or Mathematica to compute the natural frequencies. Problem #3 The structure of a carbon dioxide molecule may be idealized as three masses connected by two springs, where the masses are the carbon and oxygen atoms, and the springs represent the chemical bonds between the atoms. The equations of motion for the three atoms, i.e. masses, are given by dx1 = -631 mo 212 = -kx + kx2, at2 dt2 mc-4122 = kx1 2kx2 + ku3, most = kx2 k23. a) By substituting the oscillatory solutions xi(t) = Ajeiwt (or xi(t) = Al cos(wt)], 22(t) = Azeiwt for x2(t) = A, cos(wt)], x3(t) = Azeiwt (or x3(t) = A3 cos(wt)], into the differential equations by hand, determine the system of algebraic equations for the natural frequencies of vibration w of the system of atoms from these equations of motion. b) For k = 1,420 kg/sa, mo = 2.6568 x 10-26 kg, and mc = 1.9926 x 10-26 kg, use built-in functions in Python, Matlab, or Mathematica to compute the natural frequencies