Scribd
Upload
24 views1 page

Runge Kutta

Uploaded by

2018 01097
Copyright
© © All Rights Reserved
We take content rights seriously. If you suspect this is your content, claim it here.
Available Formats
Download as PDF, TXT or read online on Scribd
24 views1 page

Runge Kutta

Uploaded by

2018 01097
Copyright
© © All Rights Reserved
We take content rights seriously. If you suspect this is your content, claim it here.
Available Formats
Download as PDF, TXT or read online on Scribd
You are on page 1/ 1

Runge Kutta Method

In[328]:= f[x_, y_] := 2 x + y


x[0] = 0; y[0] = 1;
h = 0.2; x[n] = x[0] + n h;
Forn = 0, n ≤ 2, n++,
k1 = h f[x[n], y[n]];
1 1
k2 = h fx[n] + h, y[n] + k1;
2 2
1 1
k3 = h fx[n] + h, y[n] + k2;
2 2
k4 = h f[x[n] + h, y[n] + k3];
k1 + 2 k2 + 2 k3 + k4
k[n] = ;
6
y[n + 1] = y[n] + k;
Print["y(", x[n + 1], ")= ", y[n + 1]] // Quiet

y(0.2)= 1.2642

y(0.4)= 1.5284

y(0.6)= 1.7926

In[340]:= f[x_, y_] := 2 x + y


x[0] = 0; y[0] = 1;
h = 0.2; x[n] = x[0] + n h
Forn = 0, n ≤ 2, n++,
k1 = h f[x[n], y[n]];
1 1
k2 = h fx[n] + h, y[n] + k1;
2 2
1 1
k3 = h fx[n] + h, y[n] + k2;
2 2
k4 = h f[x[n] + h, y[n] + k3];
k1 + 2 k2 + 2 k3 + k4
k[n] = ;
6
y[n + 1] = y[n] + k;
Print["y(", x[n + 1], ")= ", y[n + 1]] // Quiet

Out[342]= 0.6

y(0.2)= 1.2642

y(0.4)= 1.5284

y(0.6)= 1.7926

You might also like