Numerical solution of initial and singularly perturbed two-point boundary value problems using adaptive spline function approximation.
We consider Sturm-Liouville differential operators on a finite interval with discontinuous potentials having one jump. As the main result we obtain a procedure of recovering the location of the discontinuity and the height of the jump. Using our result, we apply a generalized Rundell-Sacks algorithm of Rafler and Böckmann for a more effective reconstruction of the potential and present some numerical examples.
In this paper we present an iterative algorithm for the construction of Mathieu functions of any order in the form of Fourier series (practically, polynomials), and also the corresponding Quick-BASIC program for realization of this algorithm with numerical values of the parameter.