Page 10 - DOC-20220324-WA0012.-488-500
P. 10

474    Chapter 11 Systems of Ordinary Differential Equations

                         (3) The Runge-Kutta method of order 4 is
                                                         h
                                           X(t + h) = X +  (K 1 + 2K 2 + 2K 3 + K 4 )
                                                         6
                         where
                                               ⎧
                                                 K 1 = F(t, X)
                                               ⎪
                                               ⎪
                                               ⎪             1      1
                                               ⎨
                                                 K 2 = F t + h, X + hK 1
                                                             2      2

                                                             1      1
                                               ⎪ K 3 = F t + h, X + hK 2
                                               ⎪             2      2
                                               ⎪
                                               ⎩
                                                 K 4 = F(t + h, X + hK 3 )
                         Here, X = X(t), and all quantities are vectors with n components except variables t and h.
                         (4) We can absorb the t variable into the vector by letting x 0 = t and then writing the
                         autonomous form for the system of ordinary differential equations in vector notation as

                                                      X = F(X)

                                                      X(a) = S, given
                         where vectors are defined to have n + 1 components. Then
                                            ⎧                       T
                                                 X = [x 0 , x 1 , x 2 ,..., x n ]
                                            ⎪
                                            ⎪
                                            ⎪
                                                 X = [x , x , x ,..., x ]
                                            ⎨                        T
                                                       0  1  2     n
                                                 F = [1, f 1 , f 2 ,..., f n ] T
                                            ⎪
                                            ⎪
                                            ⎪
                                            ⎩                               T
                                               X(a) = [a, x 1 (a), x 2 (a),..., x n (a)]
                         (5) The Runge-Kutta method of order 4 for the system of ordinary differential equations
                         in autonomous form is
                                                         h
                                           X(t + h) = X +  (K 1 + 2K 2 + 2K 3 + K 4 )
                                                         6
                         where
                                                  ⎧
                                                     K 1 = F(X)
                                                  ⎪
                                                  ⎪
                                                  ⎪              1
                                                  ⎨
                                                     K 2 = F X + hK 1
                                                                 2

                                                                 1
                                                  ⎪ K 3 = F X + hK 2
                                                  ⎪              2
                                                  ⎪
                                                  ⎩
                                                     K 4 = F(X + hK 3 )
                         Here, X = X(t), and all quantities F and K i are vectors with n + 1 components except the
                         variables t and h.
           Problems 11.1
                         a 1. Consider

                                                 x = y               x(0) =−1

                                                           with
                                                 y = x               y(0) = 0

                            Write down the equations, without derivatives, to be used in the Taylor series method
                            of order 5.
   5   6   7   8   9   10   11   12   13