• 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!
1
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 9

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

Hausaufgaben zur Abgabe am 11.06.04

1. Bearbeiten Sie die 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. Leiten Sie die Stabilitätsfunktionen für folgende Verfahren her:

a) implizite Mittelpunktsregel b) Hammer-Hollingsworth c) SDIRK p=3

d) Gauss p=4

3. Berechnen Sie die Stabilitätsbereiche für folgende Verfahren:

a) yn+1 = yn+h2(f(tn,yn) +f(tn+1,yn+1)) Trapezregel b) yn+1 = yn+h f(tn,yn) +h2 f0(tn,yn)

Taylor-Formel 2.Ordnung c) yn+1 = yn+h f(tn+h2,yn+h2f(tn,yn))

modifizierte Euler-Formel d) yn+1 = yn+h6(4 f(tn,yn) +2 f(tn+1,yn+1) +h f0(tn,yn))

HINWEIS: Berechnen Sie jeweils zunächst die zugehörige Stabilitätsfunktion.

4. Für welche Schrittweiten wird das System

u0(t) =Au(t) mit A=

−10 9

10 −11

durch das klassische Runge-Kutta-Verfahren 4.Ordnung stabil integriert?

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