The numerical modeling of a cesium cycle in the upper atmosphere by an L-stable method of second-order 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.