We report the implementation of the complete set of the lowest-order relativistic corrections of the order of 2 where is the fine structure constant for calculating vibrational states of diatomic molecular systems within ...