Computation of the M matrix for fourth-order problems
This paper is concerned with the numerical computation of the Titchmarsh-Weyl M matrix. We show how an algorithm may be developed which relates solutions of an initial-value problem to an approximation of the M matrix. In some special cases we compare the results from our algorithm with results obtained using special functions.