A Collocation/Quadrature-Based Sturm-Liouville Problem Solver

Loading...
Thumbnail Image

Files

TR_99-1.pdf (434.85 KB)
No. of downloads: 1739

Publication or External Link

Date

1999

Advisor

Citation

DRUM DOI

Abstract

We present a computational method for solving a class of boundary-value problemsin Sturm-Liouville form. The algorithms are based on global polynomialcollocation methods and produce discrete representationsof the eigenfunctions. Error control is performed by evaluating theeigenvalue problem residuals generated when the eigenfunctions are interpolatedto a finer discretization grid; eigenfunctions thatproduce residuals exceeding an infinity-norm bound are discarded.Because the computational approach involves the generationof quadrature weights and discrete differentiation operations, our computationalmethods provide a convenient framework for solving boundary-value problemsby eigenfunction expansion and other projection methods.

Notes

Rights