• Keine Ergebnisse gefunden

Numerik II (Numerik gewöhnlicher Differentialgleichungen)

N/A
N/A
Protected

Academic year: 2021

Aktie "Numerik II (Numerik gewöhnlicher Differentialgleichungen)"

Copied!
2
0
0

Wird geladen.... (Jetzt Volltext ansehen)

Volltext

(1)

BTU COTTBUS

LEHRSTUHLNUMERISCHE MATHEMATIK UND

WISSENSCHAFTLICHESRECHNEN Prof. Dr. G. Bader, Dipl.-Math. Friedemann Kemm

Numerik II (Numerik gewöhnlicher Differentialgleichungen)

Serie 5

www.math.tu-cottbus.de/˜kemm/lehre/ode/index.html

Hausaufgaben zur Abgabe am 14.05.04

1. Bearbeiten Sie die Teile a) und b) der Programmieraufgabe, für die wir Ihnen das Material auf der Internetseite der Vorlesung hinterlegt haben. Die genaue Beschreibung der Aufgabe finden Sie in der Datei.README.

2. Man zeige für eine allgemeine RK-Formel 2. Ordnung, daß der lokale Abschneidefehlerτh(t) die Darstellung

τh(t) =τ(t)h2+O(h3)

erlaubt. Spezialisieren Sie Ihr Ergebnis auf die Heunschen Formeln 2. Ordnung.

3. Unter der Autonomisierung einer AWA y0 = f(t,y),y(t0) =y(0) versteht man das erweiterte System mit ˆy= (t,y)

ˆ

y0= fˆ(y) =ˆ

1

f(t,y)

, y(tˆ 0) = t0

y0

welches offensichtlich dieselbe Lösung besitzt. Ein Einschrittverfahren heißtinvariant gegen Au- tonomisierung, wenn die zu f bzw. ˆf gehörenden Verfahrensfunktionen F und ˆF dieselben Er- gebnisse produzieren:

Fˆ(y; h) =ˆ

1

F(t,y; h)

Zeigen Sie: Ein explizites konsistentes Runge-Kutta-Verfahren der Stufe s ist genau dann invari- ant gegen Autonomisierung, wenn gilt

ci=

i−1

j=1

ai j, 1≤is

Hausaufgabe zur Abgabe am 21.05.04

Damit Sie etwas mehr Zeit haben, stellen wir Ihnen die Programmieraufgabe für das nächste Übungs- blatt jetzt schon zur Verfügung:

1. Bearbeiten Sie Teil c) der Programmieraufgabe, für die wir Ihnen das Material auf der Inter- netseite der Vorlesung hinterlegt haben. Die genaue Beschreibung der Aufgabe finden Sie in der Datei.README.

(2)

Bonusaufgabe

1. Bearbeiten Sie Teil d) der Programmieraufgabe, für die wir Ihnen das Material auf der Inter- netseite der Vorlesung hinterlegt haben. Die genaue Beschreibung der Aufgabe finden Sie in der Datei.README.

Referenzen

ÄHNLICHE DOKUMENTE

Test the stability of the solution to (2) with respect to changes in the initial conditi- ons by changing x 0 = 80, y 0 = 30 by a unit amount in each direction (four different

Hint: The midpoint rule is only stable if |hλ n−1 (A)| < 2, where λ n−1 (A) is the lar- gest eigenvalue of A.. Implement the trapezoidal method on the linear system of

for solving ordinary differential equations and systems of ordinary differential equations using step size control, where fun is a function handle, y0 an initial condition, t =

Universität Konstanz SS 10 Fachbereich Mathematik und StatistikS. Discuss the behavior of the sequence {y n } as n

Further, represent the update using the stability function.. Further, use the Newton method to solve the non- linear problem arising from the trapezoidal rule. Plot all results into