Page 19 - 4363
P. 19

rkfixed(y,  a,  Повертає матрицю з р + 1 стовпцями і n +
                               b, n, D)      1  рядками  (р  -  кількість  рівнянь  або
                                             порядок  рівняння,  n  -  число  кроків  на
                                             інтервалі  [a,  b])  -  таблиця  розв’язків
                                             системи:  перший  стовпець  -  це  значення
                                             аргументу х, а наступні стовпці - значення
                                             ординат  розв’язку.  y  -  вектор  початкових
                                             умов  розмірності  n.  D  (x,  y)  -  функція-
                                             вектор  з  n  елементів,  що  містять  перші
                                             похідні невідомих функцій.
                                  Можна  розв’язати  задачу  більш  точно  (більш  швидко),
                            якщо зменшити крок h там, де похідна змінюється швидко, і
                            збільшити  крок  там,  де  вона  веде  себе  більш  спокійно.  Для
                            цього  передбачена  функція  Rkadapt  (adaption  -  адаптація).
                            Аргументи  і  матриця,  повернена  функцією  Rkadapt,  так,  як
                            при rkfixed (див.рис.6.1). Розв’язок системи ЗДР показано на
                            рис. 6.3 (приклад 2).
                                                       Крайові задачі
                                  Крайова  задача  формулюється  таким  чином:  нехай  на
                            відрізку  [a,  b]  потрібно  знайти  розв’язок  диференціального
                            рівняння (для простоти викладеного будемо вести на прикладі
                            ЗДР другого порядку):
                                                    d  2 y
                                                          f  (x ,  , y  ) ' y          (2.6)
                                                    dx 2
                                  при граничних умовах у (а) = А, у (b) = В.
                                  У  цьому  випадку  Mathcad  пропонує  використовувати
                            функцію sbval, щоб знайти недостатні початкові умови в точці
                            а.
                                  У  цьому  випадку  Mathcad  пропонує  використовувати
                            функцію sbval, щоб знайти відсутні початкові умови в точці а.







                                                           17
   14   15   16   17   18   19   20   21   22   23   24