Page 118 - 140

e D a  ( K i , K ) j   ( K i K , j )
(N K )  . (2.162)
L   ( K i K , j )  ( K i K , j )

Інтегруємо матричне рівняння (2.156) по часу з викорис-
танням методу трапецій


K  K  5 , 0   ( K  K  o ) , (2.163)


o

де K - невідомі значення концентрацій продуктів у кінці

кроку обчислень; K - відомі початкові значення концентра-
o
ції.
Формула (2.156) з врахуванням (2.163) набуває вигляду

2 2
[ N ( K  )  M ]K   [ N ( K o )  M ]K (2.164)
o

   

або узагальненої форми

A н ( K  )K  B н ( K o )K . (2.165)
o


Для розв'язування нелінійного рівняння (2.165) викорис-
товуємо метод Ньютона-Рафсона. Знаходимо нев'язку за ре-
зультатами i -тої апроксимації

i i i
R  A н ( K  )K  B н ( K o )K . (2.166)
o


Накладаємо таку умову

i 1 i i 1
R  R  R  0 . (2.167)

Це дає таке рекурентне співвідношення:

- 112 -