Theoretical modeling of lasers due to vibrational-rotational transitions in diatomic molecules taking anharmonicity and rotational nonequilibrium into account