Euler Method

Do 10 Steps. Solve Exactly.Compute the error. Show details.

1. y'+0.2y=0
y(0)= 5
h= 0.2

n Xn Yn Y(eksak) Error Y1=Y0+hf'(x,y)

0 0 5 5.0000 0 f'(x,y)=
1 0.2 4.8 4.8039 0.0039 f'(0,5)=
2 0.4 4.6080 4.6156 0.0076
3 0.6 4.4237 4.4346 0.0109 Y1=
4 0.8 4.2467 4.2607 0.0140 =
5 1 4.0769 4.0937 0.0168 =
6 1.2 3.9138 3.9331 0.0194
7 1.4 3.7572 3.7789 0.0217 Y (eksak)
8 1.6 3.6069 3.6307 0.0238 f(x,y)=
9 1.8 3.4627 3.4884 0.0257
10 2 3.3242 3.3516 0.0274

Do 10 Steps. Solve Exactly.Compute the error. Show details.

2. y'= 1/2phi. Sqrt(1-y^2) =0
y(0)= 0
h= 0.1

n Xn Yn Y(eksak) Error Y1=Y0+hf'(x,y)

0 0 0 0 0 f'(x,y)=
1 0.1 0.1571 0.1564 0.0006 f'(0,0)=
2 0.2 0.3122 0.3090 0.0032
3 0.3 0.4614 0.4540 0.0074 Y1=
4 0.4 0.6008 0.5878 0.0130 =
5 0.5 0.7264 0.7071 0.0193 =
6 0.6 0.8343 0.8090 0.0253
7 0.7 0.9209 0.8910 0.0299 Y (eksak)
8 0.8 0.9821 0.9511 0.0311 f(x,y)=
9 0.9 1.0117 0.9877 0.0240
10 1 1.0358 1.0000 0.0358

Do 10 Steps. Solve Exactly.Compute the error. Show details.

3. y'= (y-x)^2
y(0)= 0
h= 0.1

n Xn Yn Y(eksak) Error Y1=Y0+hf'(x,y)

0 0 0 0 0 f'(x,y)=
1 0.1 0.0000 0.0003 0.0003 f'(0,0)=
 2 0.2 0.0010 0.0026 0.0016
3 0.3 0.0050 0.0087 0.0037 Y1=
4 0.4 0.0137 0.0201 0.0064 =
5 0.5 0.0286 0.0379 0.0093 =
6 0.6 0.0508 0.0630 0.0121
7 0.7 0.0810 0.0956 0.0147 Y (eksak)
8 0.8 0.1193 0.1360 0.0167 f(x,y)=
9 0.9 0.1656 0.1837 0.0181
10 1 0.2196 0.2384 0.0188

Do 10 Steps. Solve Exactly.Compute the error. Show details.

4. y'= (y+x)^2
y(0)= 0
h= 0.1

n Xn Yn Y(eksak) Error Y1=Y0+hf'(x,y)

0 0 0 0.00000 0 f'(x,y)=
1 0.1 0 0.00033 0.0003 f'(0,0)=
2 0.2 0.0010 0.00271 0.0017
3 0.3 0.0050 0.00934 0.0043 Y1=
4 0.4 0.0143 0.02279 0.0084 =
5 0.5 0.0315 0.04630 0.0148 =
6 0.6 0.0598 0.08414 0.0244
7 0.7 0.1033 0.14229 0.0390 Y (eksak)
8 0.8 0.1678 0.22964 0.0618 f(x,y)=
9 0.9 0.2615 0.36016 0.0987
10 1 0.3964 0.55741 0.1610

Improved Euler Method

Do 10 Steps. Solve Exactly.Compute the error. Show details.

5. y'= y
y(0)= 1
h= 0.1

Error of Error of
n Xn Yn* Predictor Yn Corrector Y(eksak) Improved Euler
Euler
0 0 1 1 1.0000 0 0 Yn*=y0+h(yo)
1 0.1 1.1000 1.1050 1.1052 0.0002 0.0052 Y1*=
2 0.2 1.2155 1.2210 1.2214 0.0004 0.0059 =
3 0.3 1.3431 1.3492 1.3499 0.0006 0.0067 Yn=Y0+h.1/2 ((yo)+(Yn*
4 0.4 1.4842 1.4909 1.4918 0.0009 0.0077 Y1=
5 0.5 1.6400 1.6474 1.6487 0.0013 0.0087 =
6 0.6 1.8122 1.8204 1.8221 0.0017 0.0099 =
7 0.7 2.0025 2.0116 2.0138 0.0022 0.0113
 8 0.8 2.2127 2.2228 2.2255 0.0028 0.0128 Y (eksak)
9 0.9 2.4451 2.4562 2.4596 0.0034 0.0145 f(x,y)=
10 1 2.7018 2.7141 2.7183 0.0042 0.0165

Do 10 Steps. Solve Exactly.Compute the error. Show details.

6. y'= 2(1+y^2)
y(0)= 0
h= 0.05

Error of Error of
n Xn Yn* Predictor Yn Corrector Y(eksak) Improved Euler
Euler
0 0 0 0 0 0 0 Yn*=y0+h(2(1+yo^2)
1 0.05 0.1000 0.1005 0.1003 0.0002 0.0003 Y1*=
2 0.1 0.2015 0.2030 0.2027 0.0003 0.0012 =
3 0.15 0.3072 0.3098 0.3093 0.0005 0.0022 Yn=Y0+h.1/2 ((2(1+yo^
4 0.2 0.4194 0.4234 0.4228 0.0006 0.0034 Y1=
5 0.25 0.5413 0.5470 0.5463 0.0007 0.0050 =
6 0.3 0.6769 0.6849 0.6841 0.0008 0.0072 =
7 0.35 0.8318 0.8429 0.8423 0.0007 0.0105
8 0.4 1.0140 1.0299 1.0296 0.0002 0.0156 Y (eksak)
9 0.45 1.2360 1.2593 1.2602 0.0009 0.0242 f(x,y)=
10 0.5 1.5179 1.5538 1.5574 0.0036 0.0395

Do 10 Steps. Solve Exactly.Compute the error. Show details.

7. y'-xy^2=0
y(0)= 1
h= 0.1

Error of Error of
n Xn Yn* Predictor Yn Corrector Y(eksak) Improved Euler
Euler
0 0 1 1 1.0000 0 0 Yn*=y0+h(x0*y0^2)
1 0.1 1.0000 1.0050 1.0050 0.0000 0.0050 Y1*=
2 0.2 1.0151 1.0204 1.0204 0.0001 0.0053 =
3 0.3 1.0412 1.0470 1.0471 0.0001 0.0059 Yn=Y0+h.1/2 ((x0*y0^2
4 0.4 1.0799 1.0868 1.0870 0.0002 0.0070 Y1=
5 0.5 1.1340 1.1426 1.1429 0.0003 0.0088 =
6 0.6 1.2078 1.2190 1.2195 0.0005 0.0117 =
7 0.7 1.3081 1.3234 1.3245 0.0011 0.0164
8 0.8 1.4460 1.4684 1.4706 0.0022 0.0245 Y (eksak)
9 0.9 1.6409 1.6758 1.6807 0.0049 0.0398 f(x,y)=
10 1 1.9285 1.9881 2.0000 0.0119 0.0715

Do 10 Steps. Solve Exactly.Compute the error. Show details.

8. y'=y-y^2
y(0)= 0.2
 h= 0.1

Error of Error of
n Xn Yn* Predictor Yn Corrector Y(eksak) Improved Euler
Euler
0 0 0.2 0.2 0.2000 0 0 Yn*=y0+h(y0-y0^2)
1 0.1 0.2160 0.2165 0.2165 0.0000 0.0005 Y1*=
2 0.2 0.2334 0.2339 0.2339 0.0000 0.0005 =
3 0.3 0.2518 0.2523 0.2523 0.0000 0.0005 Yn=Y0+h.1/2 ((x0*y0^2
4 0.4 0.2711 0.2716 0.2716 0.0001 0.0005 Y1=
5 0.5 0.2914 0.2918 0.2919 0.0001 0.0005 =
6 0.6 0.3125 0.3129 0.3130 0.0001 0.0005 =
7 0.7 0.3344 0.3348 0.3349 0.0001 0.0005
8 0.8 0.3570 0.3574 0.3575 0.0001 0.0005 Y (eksak)
9 0.9 0.3803 0.3806 0.3808 0.0001 0.0004 f(x,y)=
10 1 0.4042 0.4045 0.4046 0.0001 0.0004
 Grafik Euler Method y'+0.2y=0
5

4.8

4.6

4.4
Y1=Y0+hf'(x,y)
-0.2y 4.2

-0.2(5) = -1 4

Y
3.8
Y0 + h f'(x,y) 3.6
5 + 0.2 -1
3.4
4.8
3.2

3
0 0.1 0.2 0.3 0.4 0.5 0.6 0.7 0.8 0.9 1 1.1 1.2 1.3 1.4 1.5 1.6 1.7
5/e^(x/5)
X

Grafik Euler Method y'= 1/2phi. Sqrt(1

1.4

1.2

Y1=Y0+hf'(x,y)
1
1/2phi rad.sqrt(1-y^2)
'1/2phi rad.sqrt(1-0^2) = 1.570796 0.8
Y

Y0 + h f'(x,y) 0.6
0 + 0.1 1.570796
0.15707963267949 0.4

0.2

sin((phi rad*x)/2) 0
0 0.1 0.2 0.3 0.4 0.5 0.6 0.7 0.8

Grafik Euler Method y'= (y-x)^2

0.3

Y1=Y0+hf'(x,y) 0.25
(y-x)^2
(0-0)^2= 0 0.2

0.15
Y
 0.25

0.2

Y0 + h f'(x,y)
0.15
0 + 0.1 0

Y
0
0.1

x-tanhx 0.05

0
0 0.1 0.2 0.3 0.4 0.5 0.6 0.7 0.8

Grafik Euler Method y'= (y+x)^2

0.6
Y1=Y0+hf'(x,y)
(y+x)^2 0.55

(0+0)^2= 0 0.5

0.45
Y0 + h f'(x,y) 0.4
0 + 0.1 0 0.35
0
0.3
Y

0.25

0.2
tanx-x
0.15

0.1

0.05

0
0 0.1 0.2 0.3 0.4 0.5 0.6 0.7 0.8

Grafi
2.8
2.7
2.6
2.5
2.4
2.3
Yn*=y0+h(yo) 2.2
2.1
1+ 0.1 1
2
1.1 1.9
Y

Yn=Y0+h.1/2 ((yo)+(Yn*)) 1.8

Y0 + h 0.5 ( y0 + y1* ) 1.7
1.6
1 + 0.1 0.5 ( 1+ 1.1 ) 1.5
1.1050 1.4
1.3
1.2
1.1
1
0 0.1 0.2 0.3
 1.7
1.6
1.5
1.4
1.3
1.2
1.1
e^x 1
0 0.1 0.2 0.3

Grafik Im
1.6
1.5
Yn*=y0+h(2(1+yo^2) 1.4
0+ 0.05 2 1.3
0.1 1.2
Yn=Y0+h.1/2 ((2(1+yo^2))+(2(1+yn*^2))) 1.1
1
Y0 + h 0.5 ( (2(1+yo^2) + (2(1+yn*^2)) )
0.9
0 + 0.05 0.5 ( 2+ 2.02 ) 0.8

Y
0.1005 0.7
0.6
0.5
0.4
tan(2x)
0.3
0.2
0.1
0
0 0.1 0

Grafik Imp
2
Yn*=y0+h(x0*y0^2)
1.9
1+ 0.10 0
1 1.8

Yn=Y0+h.1/2 ((x0*y0^2))+((x1*(y1*)^2))) 1.7

Y0 + h 0.5 ( (x0*y0^2) + ((x1*(y1*)^2)) )
1.6
1 + 0.1 0.5 ( 0+ 0.1 )
1.0050 1.5
Y

1.4

1.3
2/(2-x^2)
1.2

1.1

1
0 0.1 0.2 0.3 0.
 1.1

1
0 0.1 0.2 0.3 0.

Grafik Im
Yn*=y0+h(y0-y0^2)
0.2 + 0.10 0.16 0.44
0.216 0.42
Yn=Y0+h.1/2 ((x0*y0^2))+((x1*(y1*)^2))) 0.4
Y0 + h 0.5 ( (y0-y0^2) + (((y1*)-(y1*)^2)) )
0.38
0.2 + 0.1 0.5 ( 0.16 + 0.169344 )
0.36
0.2165
0.34

0.32

Y
e^x/(4+e^x) 0.3
0.28
0.26
0.24
0.22
0.2
0 0.1 0.2 0.3
er Method y'+0.2y=0

Yn
Y(eksak)

1 1.1 1.2 1.3 1.4 1.5 1.6 1.7 1.8 1.9 2 2.1

od y'= 1/2phi. Sqrt(1-y^2) =0

Yn
Y(eksak)

0.5 0.6 0.7 0.8 0.9 1

r Method y'= (y-x)^2

Yn
Y(eksak)
 Yn
Y(eksak)

0.5 0.6 0.7 0.8 0.9 1

Method y'= (y+x)^2

Yn
Y(eksak)

0.5 0.6 0.7 0.8 0.9 1

Grafik Improved Euler Method y'= y

2.8
2.7
2.6
2.5
2.4
2.3
2.2
2.1
2 Yn Euler
1.9 Yn Improved Euler
1.8
1.7
Y(eksak)
1.6
1.5
1.4
1.3
1.2
1.1
1
0 0.1 0.2 0.3 0.4 0.5 0.6 0.7 0.8 0.9 1

X
 1.7
1.6
1.5
1.4
1.3
1.2
1.1
1
0 0.1 0.2 0.3 0.4 0.5 0.6 0.7 0.8 0.9 1

Grafik Improved Euler Method y'= 2(1+y^2)

1.6
1.5
1.4
1.3
1.2
1.1
1
0.9 Yn Euler
0.8 Yn Improved Euler
0.7
Y(eksak)
0.6
0.5
0.4
0.3
0.2
0.1
0
0 0.1 0.2 0.3 0.4 0.5

Grafik Improved Euler Method y'-xy^2= 0

2

1.9

1.8

1.7

1.6
Yn Euler
1.5 Yn Improved Euler
1.4 Y(eksak)
1.3

1.2

1.1

1
0 0.1 0.2 0.3 0.4 0.5 0.6 0.7 0.8 0.9 1

X
1.1

1
0 0.1 0.2 0.3 0.4 0.5 0.6 0.7 0.8 0.9 1

Grafik Improved Euler Method y'-y^2= 0

0.44
0.42
0.4
0.38

0.36
0.34 Yn Euler
0.32 Yn Improved Euler
0.3 Y(eksak)
0.28
0.26
0.24
0.22
0.2
0 0.1 0.2 0.3 0.4 0.5 0.6 0.7 0.8 0.9 1