Abstract
Spherical splines are used to define approximate solutions to strongly elliptic pseudodiffererential equations on the unit sphere. These equations arise from geodesy. The approximate solutions are found by using Galerkin method. We prove optimal convergence (in Sobolev norms) of the approximate solution by spherical splines to the exact solution. Our numerical results underlie the theoretical result.