A pure two-body problem has seven integrals including the Kepler energy, the Laplace vector, and the angular momentum vector. However, only five of them are independent. When the five independent integrals are preserved, the two other dependent integrals are naturally preserved from a theoretical viewpoint; but they may not be either from a numerical computational viewpoint. Because of this, we use seven scale factors to adjust the integrated positions and velocities so that the adjusted solutions strictly satisfy the seven constraints. Noticing the existence of the two dependent integrals, we adopt the Newton iterative method combined with the singular value decomposition to calculate these factors. This correction scheme can be applied to perturbed two-body and N-body problems in the solar system. In this case, the seven quantities of each planet slowly vary with time. More accurate values can be given to the seven slowly-varying quantities by integrating the integral invariant relations of these quantities and the equations of motion. They should be satisfied with the adjusted solutions. Numerical tests show that the new method can significantly reduce the rapid growth of numerical errors of all orbital elements.

For most hierarchical triple stars, the classical double two-body model of zeroth-order cannot describe the motions of the components under the current observational accuracy. In this paper, Marchal's first-order analytical solution is implemented and a more efficient simplified version is applied to real triple stars. The results show that, for most triple stars, the proposed first-order model is preferable to the zeroth-order model both in fitting observational data and in predicting component positions.

This work is devoted to improving empirical mass-luminosity relations and mass-metallicity-luminosity relation for low mass stars. For these stars, observational data in the mass-luminosity plane or the mass-metallicity-luminosity space subject to non-negligible errors in all coordinates with different dimensions. Thus a reasonable weight assigning scheme is needed for obtaining more reliable results. Such a scheme is developed, with which each data point can have its own due contribution. Previous studies have shown that there exists a plateau feature in the mass-luminosity relation. Taking into account the constraints from the observational luminosity function, we find by fitting the observational data using our weight assigning scheme that the plateau spans from 0.28 to 0.50 solar mass. Three-piecewise continuous improved mass-luminosity relations in K, J, H and V bands, respectively, are obtained. The visual mass-metallicity-luminosity relation is also improved based on our K band mass-luminosity relation and the available observational metallicity data.