Based on topography and gravitational field data, model variations of the thickness of the crust of Mars and Venus are calculated by the method of Love numbers. The method takes into account the adjustment of the planetary subsurface to the loads on the surface and in the subsurface. Numerical modeling is carried out using spherical harmonic expansion of topography and gravitational field data up to the 90th degree and order for Mars and up to the 70th for Venus. The relief of the crust-mantle boundary suggests partial Airy isostatic compensation. The model of the crust of Mars is consistent with the range of values of the crustal thickness beneath the InSight mission station in the southwestern part of the Elysium Planitia, obtained from the results of the seismic experiment.