Log in Register Contribute

Custom Search Ask a Question

MATH NOTES CALCULATORS WEBASSIGN ANSWERS MATH GAMES AND LOGIC PUZZLES

SOLVED QUESTIONS

Home /  Calculators /  Differential Equations Calculators /  Math Problem Solver (all calculators)

IMPROVED EULER (HEUN'S) METHOD CALCULATOR

The calculator will nd the approximate solution of the rst-order differential equation using the
Improved Euler (Heun's) method, with steps shown.

Fiken
ÅPNE
Faktura, timeføring, altinn

SHOW INSTRUCTIONS

Enter a function: y ′ = f (x, y) or y ′ = f (t, y) =

-y+t+1

Enter the Step Size h= : 0.1

Enter the initial condition: y ( 0 )= 1

Find y ( 1 )
If the calculator did not compute something or you have identi ed an error, please write it in
comments below.

Write all suggestions in comments below.

Show steps

CALCULATE

SOLUTION
1
Your input: nd y (1) for y ′ = t− y + 1 , when y (0) = 1 h =, using the Improved
10
Euler's method.

h
~
The Improved Euler method states that yn+1 = y
n
+ (f (tn , y ) + f (tn+1 , y
n n+1
)) , where
2
~
y
n+1
= y
n
+ h ⋅ f (tn , y )
n
and tn+1 = tn + h .

1
We have that h = , t0 = 0 y , 0
,
= 1 f (t, y) = t − y + 1
10

Step 1.

1 1
t1 = t0 + h = 0 + =
10 10

~
y = y + h ⋅ f (t0 , y ) = 1 + h ⋅ f (0, 1) =
1 0 0

1
= (1) + ⋅ (0) = 1
10

h
~
y = y + (f (t0 , y ) + f (t1 , y )) =
1 0 0 1
2
 h 1
= 1 + (f (0, 1) + f ( , 1)) =
2 10

10
1
= 1 + (0 + ( )) = 1.005
2 10

Step 2.

1 1 1
t2 = t1 + h = + =
10 10 5

1
~
y = y + h ⋅ f (t1 , y ) = 1.005 + h ⋅ f ( , 1.005) =
2 1 1
10

1
= (1.005) + ⋅ (0.0950000000000002) = 1.0145
10

h
~
y2 = y1 + (f (t1 , y1 ) + f (t2 , y 2 )) =
2

h 1 1
= 1.005 + (f ( , 1.005) + f ( , 1.0145)) =
2 10 5

10
= 1.005 + (0.0950000000000002 + (0.1855)) = 1.019025
2

Step 3.

1 1 3
t3 = t2 + h = + =
5 10 10

1
~
y = y + h ⋅ f (t2 , y ) = 1.019025 + h ⋅ f ( , 1.019025) =
3 2 2
5

1
= (1.019025) + ⋅ (0.180975) = 1.0371225
10

h
~
y = y + (f (t2 , y ) + f (t3 , y )) =
3 2 2 3
2

h 1 3
= 1.019025 + (f ( , 1.019025) + f ( , 1.0371225)) =
2 5 10
 1

10
= 1.019025 + (0.180975 + (0.2628775)) = 1.041217625
2

Step 4.

3 1 2
t4 = t3 + h = + =
10 10 5

3
~
y = y + h ⋅ f (t3 , y ) = 1.041217625 + h ⋅ f ( , 1.041217625) =
4 3 3
10

1
= (1.041217625) + ⋅ (0.258782375) = 1.0670958625
10

h
~
y = y + (f (t3 , y ) + f (t4 , y )) =
4 3 3 4
2

h 3 2
= 1.041217625 + (f ( , 1.041217625) + f ( , 1.0670958625)) =
2 10 5

10
= 1.041217625 + (0.258782375 + (0.3329041375)) = 1.070801950625
2

Step 5.

2 1 1
t5 = t4 + h = + =
5 10 2

2
~
y = y + h ⋅ f (t4 , y ) = 1.070801950625 + h ⋅ f ( , 1.070801950625) =
5 4 4
5

1
= (1.070801950625) + ⋅ (0.329198049375) = 1.1037217555625
10

h
~
y = y + (f (t4 , y ) + f (t5 , y )) =
5 4 4 5
2

h 2 1
= 1.070801950625 + (f ( , 1.070801950625) + f ( , 1.1037217555625))
2 5 2
=

10
= 1.070801950625 + (0.329198049375 + (0.3962782444375))
2

= 1.10707576531562
Step 6.

1 1 3
t6 = t5 + h = + =
2 10 5

1
~
y = y + h ⋅ f (t5 , y ) = 1.10707576531562 + h ⋅ f ( , 1.10707576531562) =
6 5 5
2

1
= (1.10707576531562) + ⋅ (0.392924234684375) = 1.14636818878406
10

h
~
y = y + (f (t5 , y ) + f (t6 , y )) =
6 5 5 6
2

= 1.10707576531562

h 1 3
+ (f ( , 1.10707576531562) + f ( , 1.14636818878406)) =
2 2 5

10
= 1.10707576531562 + (0.392924234684375 + (0.453631811215938))
2

= 1.14940356761064

Step 7.

3 1 7
t7 = t6 + h = + =
5 10 10

3
~
y = y + h ⋅ f (t6 , y ) = 1.14940356761064 + h ⋅ f ( , 1.14940356761064) =
7 6 6
5

1
= (1.14940356761064) + ⋅ (0.45059643238936) = 1.19446321084958
10

h
~
y = y + (f (t6 , y ) + f (t7 , y )) =
7 6 6 7
2

= 1.14940356761064

h 3 7
+ (f ( , 1.14940356761064) + f ( , 1.19446321084958)) =
2 5 10

10
= 1.14940356761064 + (0.45059643238936 + (0.505536789150424))
2

= 1.19721022868763
Step 8.

7 1 4
t8 = t7 + h = + =
10 10 5

7
~
y = y + h ⋅ f (t7 , y ) = 1.19721022868763 + h ⋅ f ( , 1.19721022868763) =
8 7 7
10

1
= (1.19721022868763) + ⋅ (0.50278977131237) = 1.24748920581887
10

h
~
y = y + (f (t7 , y ) + f (t8 , y )) =
8 7 7 8
2

= 1.19721022868763

h 7 4
+ (f ( , 1.19721022868763) + f ( , 1.24748920581887)) =
2 10 5

10
= 1.19721022868763 + (0.50278977131237 + (0.552510794181133))
2

= 1.2499752569623

Step 9.

4 1 9
t9 = t8 + h = + =
5 10 10

4
~
y = y + h ⋅ f (t8 , y ) = 1.2499752569623 + h ⋅ f ( , 1.2499752569623) =
9 8 8
5

1
= (1.2499752569623) + ⋅ (0.550024743037695) = 1.30497773126607
10

h
~
y = y + (f (t8 , y ) + f (t9 , y )) =
9 8 8 9
2

= 1.2499752569623

h 4 9
+ (f ( , 1.2499752569623) + f ( , 1.30497773126607)) =
2 5 10

10
= 1.2499752569623 + (0.550024743037695 + (0.595022268733925))
2

= 1.30722760755089
Step 10.

9 1
t10 = t9 + h = + = 1
10 10

9
~
y = y + h ⋅ f (t9 , y ) = 1.30722760755089 + h ⋅ f ( , 1.30722760755089) =
10 9 9
10

1
= (1.30722760755089) + ⋅ (0.592772392449114) = 1.3665048467958
10

h
~
y = y + (f (t9 , y ) + f (t10 , y )) =
10 9 9 10
2

= 1.30722760755089

h 9
+ (f ( , 1.30722760755089) + f (1, 1.3665048467958)) =
2 10

10
= 1.30722760755089 + (0.592772392449114 + (0.633495153204203))
2

= 1.36854098483355

Answer: y (1) = 1.36854098483355

If you like the website, please share it anonymously with your friend or teacher by
entering his/her email: SHARE

About Contact Terms of use Privacy Policy Links Request a Note

Copyright (c) 2020. All rights reserved.