Rovnice a hraniční podmínky popisující slapové deformace a gravitační potenciál předpjatých elastických sféricky symetrických těles tvoří sadu okrajových úloh, které použitím pseudospektrálních schémat na čebyševovské síti přejdou na soustavu lineárních algebraických rovnic. Metoda je použita na výpočet slapových Loveových čísel modelů Marsu a Venuše