20.07.2013 Aufrufe

Lösungen zu nicht besprochenen Aufgaben der 5. Übung

Lösungen zu nicht besprochenen Aufgaben der 5. Übung

Lösungen zu nicht besprochenen Aufgaben der 5. Übung

MEHR ANZEIGEN
WENIGER ANZEIGEN

Erfolgreiche ePaper selbst erstellen

Machen Sie aus Ihren PDF Publikationen ein blätterbares Flipbook mit unserer einzigartigen Google optimierten e-Paper Software.

ist aber, dass die Näherunfswerte wesentlich besser sind wie wir gleich sehen werden.<br />

Konkret für unser Beispiel ist f(x,y) = y, sodass sich bei uns die ki wie folgt berechnen:<br />

k1 = yn,<br />

k2 = yn + 1<br />

2 hk1,<br />

k3 = yn + 1<br />

2 hk2,<br />

k4 = yn +hk3.<br />

Wir führen nun das Runge-Kutta-Verfahren <strong>zu</strong>mindest mal für die Schrittweite h = 1/2<br />

durch.<br />

In diesem Fall ist x1 = 1<br />

2 . Bevor y1 berechnet werden kann, müssen die einzelnen ki<br />

berechnet werden:<br />

k1 = y0 = 1,<br />

k2 = y0 + 1<br />

2 hk1 = 1+ 1 1<br />

· ·1 = 1.25,<br />

2 2<br />

k3 = y0 + 1<br />

2 hk2 = 1+ 1 1<br />

· ·1.25 = 1.3125,<br />

2 2<br />

k4 = y0 +hk3 = 1+ 1<br />

·1.3125 = 1.6562<strong>5.</strong><br />

2<br />

Nun sind wir in <strong>der</strong> Lage, y1 aus<strong>zu</strong>rechnen:<br />

<br />

1<br />

y1 = y0 +h·<br />

6 k1 + 1<br />

3 k2 + 1<br />

3 k3 + 1<br />

6 k4<br />

<br />

= 1+ 1<br />

2 ·<br />

<br />

1 1 1 1<br />

·1+ ·1.25+ ·1.3125+<br />

6 3 3 6 ·1.65625<br />

<br />

= 1.648437<strong>5.</strong><br />

Dies ist also eine Näherung für y( 1<br />

2 ) – den wahren Funktionswert an <strong>der</strong> Stelle x1. Weiter<br />

istx2 = 1. Umy2 <strong>zu</strong> berechnen, sind erneut dieki <strong>zu</strong> ermitteln – nur dieses Mal ausgehend<br />

von (x1,y1):<br />

k1 = y1 = 1.6484375,<br />

k2 = y1 + 1<br />

2 hk1 = 1.6484375+ 1 1<br />

· ·1.6484375 = 2.06054688,<br />

2 2<br />

k3 = y1 + 1<br />

2 hk2 = 1.6484375+ 1 1<br />

· ·2.06054688 = 2.16357422,<br />

2 2<br />

k4 = y1 +hk3 = 1.6484375+ 1<br />

·2.16357422 = 2.73022461.<br />

2<br />

Für y2 ergibt sich:<br />

<br />

1<br />

y2 = y1 +h·<br />

6 k1 + 1<br />

3 k2 + 1<br />

3 k3 + 1<br />

6 k4<br />

<br />

= 1.6484375+ 1<br />

2 ·<br />

<br />

1 1 1 1<br />

·1.6484375+ ·2.06054688+ ·2.16357422+<br />

6 3 3 6 ·2.73022461<br />

<br />

= 2.71734619.<br />

8

Hurra! Ihre Datei wurde hochgeladen und ist bereit für die Veröffentlichung.

Erfolgreich gespeichert!

Leider ist etwas schief gelaufen!