Page 100 - 203_
P. 100

2 k  2     b   a
                                         f   dxx    y 2k    4y 2 k  1   y 2 k  2  ,
                                                   6n
                                       2 x k
        де  y   f   , kx    , 2 , 1 , 0    2 ,  . n
             k      k
           Додаючи почленно ці наближені рівності, одержимо наближену формуму
         b
                   b   a
           f   dxx    y 0   y 2n    2    y 2   y 4    y 2 n  2   4  y 1   y 3    y 2n  1   .
         a          6n
           Ця формула називається формулою парабол або формулою Сімпсона.
           Геометрично зміст формули Сімпсона очевидний: площа криволінійної трапеції під
        графіком функції   xf   на відрізку  ba,   наближено замінюється сумою площ фігур, яка
        лежать під параболами (прямими).
           Можна показати, що якщо функція  xf    має на  ba,   неперервну похідну, то абсолютна
        величина похибки формули Сімпсона не більша за
                                                            5
                                                      b    a
                                                   M          ,
                                                      2880n  4
        де  M  - найбільше значення  f   IV    x  на відрізку  ba,  . Вище відзначалось, що похибка
        формул прямокутників і трапецій оцінюється числом
                                               ab   3           
                                            k         k  max  f    x
                                                                   
                                              12 n 2      ba,      .
                      4
                                             2
           Оскільки  n  росте швидше, ніж  n , то похибка формули Сімпсона з ростом n
        зменшується значно швидше. Цим пояснюється, що формула Сімпсона дозволяє отримати
        більшу точність в порівнянні з формулами прямокутників і трапецій.

                                                                    1  dx
           Приклад 2. Обчислити за формулою Сілпсона інтеграл             при n =4.
                                                                    0  1  x
           Розіб’ємо відрізок  1,0  на чотири рівних частини точками
         x    , 0 x    1  , x    1  , x    3  , x    . 1  Обчислемо наближено значення функції   xf    1   в
          0     1    4   2    2  3    4   4                                                      1   x
        цих точках:
             y    , 1  000 , y    , 0  8000 , y    , 0  6667 , y    , 0  5714 , y    , 0  5000 .
             0         1           2            3           4
           За формулою Сімпсона одержимо
         b         b   a
           f   dxx    y 0   y 4    2 y 2     4 y 1   y 3 
         a          6n
         1  0
              000,1    , 0  5000  2  , 0   6667   4  8000,0    , 0  5714  0  , 69325 .
          12
                                                                                           1
           Оцінимо похибку одержаного результата. Для підінтегральної функції  xf            маємо
                                                                                          1   x
                     24
         f   IV   x   .
                   1 x   6
           Звідси випливає, що на відрізку  1,0   f   IV    x  24. Отже, можна взяти M =24, і похибка

                                                24
        результата не перевищує величини                 . 0  0004 . Порівнюючи наближене значення з
                                              2880  4   4
        тосним, робимо висновок, що абсолютна похибка результата, отриманого за формулою
        Сімпсона, менша за 0,00011. Це відповідає даній вище оцінці похибки і, крім таго, свідчить,
        що формула Сімпсона значно точніша і порівнянні з формулами прямокутників і трапецій.
   95   96   97   98   99   100   101   102   103   104   105