We numerically construct static and spherically symmetric electrically charged black hole solutions in Einstein-Born-Infeld gravity with massive dilaton. The numerical solutions show that the dilaton potential allows many more black hole causal structures than the massless dilaton. We find that depending on the black hole mass and charge and the dilaton mass the black holes can have either one, two, or three horizons. The extremal solutions are also found out. As an interesting peculiarity we note that there are extremal black holes with an inner horizon and with triply degenerated horizon.

We investigate numerically models of the static spherically symmetric boson-fermion stars in the scalar-tensor theory of gravity with massive dilaton field. The proper mathematical model of such stars is interpreted as a nonlinear two-parametric eigenvalue problem with unknown internal boundary. To solve this problem the Continuous Analogue of Newton Method is used.

We develop the general theory of stars in Saa's model of gravity with propagating torsion and study the basic stationary state of neutron star. Our numerical results show that the torsion force decreases the role of the gravity in the star configuration leading to significant changes in the neutron star masses depending on the equation of state of star matter. The inconsistency of the Saa's model with Roll-Krotkov-Dicke and Braginsky-Panov experiments is discussed.