We have studied two probability potentials of a collinear oscillator atom chain and developed a Langevin dynamics approach for calculation. In the case of the harmonic chains, results of the Monte Carlo simulations are compared with the analytical solutions to verify the validity of this approach. In the case of 4-times coupled oscillator chains, the results of numerical simulations are used to the calculation of Hamiltonian. Then the system's energy distribution and the Maxwell-Boltzmann distribution are compared, and found to be in agreement with each other.