Numerical simulation of the cesium cycle in the upper atmosphere by means of L -stable method of the second order of accuracy

An algorithm of right-hand side and Jacobian formation of differential equations of chemical kinetics is described. Numerical simulation of the cesium cycle in the upper atmosphere is conducted by means of the L-stable method of the second order of accuracy with the control accuracy. The results of the computation are presented.