Труды КНЦ вып.38 (ГЕЛИОФИЗИКА вып. 4/2016(38))

Обратная V-l матрица (І+“Г‘ имеет вид + <D j = — -— j - <Ь+ су ® су^ , где а ) 0 ( о — диада, образованная вектором т . 1+ |т |' Поэтому решение векторного уравнения (16) дается формулой / Л ѵ_1 1 и = + <эj W = ^ 2 ( і ѵ -[(О x f V] + (fV :o j)o j) Применяя эту формулу к Урав„е„„ю ( |6 ) при т = | Г 'В ' , п°лУч„м выражение для и через остальные слагаемые: ч2 W —- е ' Г B 1 X W ] + [ тГ I (W ;B 1) B 1 и = 1+ [ - Г 1 2 B 1 О ( т ) . Подставляя в эту формулу уже найденные Г 1 ) , В і(0) , В 1 0 , получим следующие формулы: W =W 0+ W 1(0), W 1(0) = т е'В1(0). и1(0) = W —те 'Г 1(—) [B 1(0) x W ] +( т Г 1(—) I ( w ;B 1(0)) B ?1(0) ^ і + [ т г 1(—) | b 1(0)' 2 O (т3) . (17) I.4. Используя полученные u 1(-0) и B l(0), В 1(0), с точностью О ( т 3) находим Г 1(0) по формуле (І0), а также к d t j по формуле ( І І ) и d t j по формуле (І3). Присваиваем д0 = 1. II. Обычная итерация К —— К + 1. Пусть уже известны X l K , и 1(К). г 1(к) в1 (К е 1(к) ^ d u Г ’ d 2 ^ 1<К| У d t j к d t j с точностью II.1. Находим XlK+1^ с точностью O ( т 4) из соотношения (І4): X1K+1) = X0 + т ( г ѵ + г 1(к)и1(к)) — — 2 ' ’ 12 CTX к d t j d y к d t j O ( т 4) . 2 150