Lunar numerical theory and determination of parameters k2, δM from analysis of LLR data

Analysis of lunar ranging data of 1970-2001 was carried out to estimate values of the lunar phase angle delta(M) = 0.0446 +/- 0.0001 and Love's number k(2) = 0.0205 +/- 0.0001. The Moon's dynamic theory has been constructed by simultaneous numerical integration of equations of motion of the major and minor planets and the Moon with equations of lunar librations. Partial derivatives with respect to estimated parameters were obtained following the method of arbitrary constant varying.