Wir betrachten numerische Lösungsverfahren für singuläre Randwertprobleme der Form
z'(t)=M(t)z(t)/t+f(t,z(t)), t∈ (0,1],
Baz(0)+Bbz(1)=β,
(1)
Speziell wird dabei auch auf ein von uns entwickeltes Matlab package eingegangen. Unsere Lösungsroutine basiert auf Kollokationsverfahren, wobei zur Fehlerschätzung ein neu entwickelter Schätzer des globalen Fehlers herangezogen wird. Diese Fehlerschätzung erweist sich als robust gegenüber der Singularität und ermöglicht eine Gittersteuerung, die sich an der Glattheit der Lösung orientiert und nicht durch das unbeschränkte Richtungsfeld beeinträchtigt wird. Die Konstruktion dieses Fehlerschätzers wurde inspiriert durch die Idee der Iterierten Defektkorrektur (IDeC), die unter anderem auf Stetter ([2]) zurückgeht. Eine detaillierte Beschreibung der Algorithmen findet sich in [1].

References

[1] Auzinger, W., O. Koch and E. Weinmüller, Efficient Collocation Schemes for Singular Boundary Value Problems. Submitted to Numer. Algorithms.

[2] Stetter, H.J., The Defect Correction Principle and Discretization Methods, Numer. Math., 27:425-443, 1978.