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

t > t 0 . Обозначим t 12 = t 0 + т / 2 , для функций вида f ( t ) и F ( x , t ) введем обозначения f a= f ( t a ) и F a= F ( x ( t a ) , t a ) , a = 0 , 1 , 1 . I. Стартовая итерация 1.1. Находим поля E 0= e ( x 0, t ) , t 0= B ^ x 0, t ) , затем по ним и формулам (8), находим Ь° , , 0 С° , Е® = ( е °;Ь° ^6 ° , ѵ® = (v°;b° ^6°, 1.2. По первой формуле в (9) находим x 12(0) = X ( t 12, t 0, x 0, v0) ,[ b° - < ] VL = v” - V l- vE в постоянных полях E 0, B 0 : 12 ( 0 ) 0 / 0 0 \ т т ir.0 V L ' ч ^ >= x +[Vp +v,| - + — e E ; + — si 1 E J2 8 1 0 0 т V 2 , 1- cos Находим поля E l2(0) = E (x l2(0), t 2) , t = b ( x 120), t ) : по ним и формулам (8), находим b V2(0) l/2(0) а затем ® W 0), Ei У20) = (e V2(0); b ^2(0)) ^ 0) j 2 ( 0 ) _ l ^ 0 . hP(0)\ ^2(0) ^ / 2(0) = v0 - ^ 2 ^ _ v | 2(0) = 1 E = (v0;b l2(0)) = v - v. ■vc I.3. По первой формуле в (9) находим x 1(0)= x ( t х, t 0, x 0, v0) в постоянных полях E 12(0), B12(0): д№) = х 0 + (vV2W+vV2W y + ^ e,E l/2W + 1/2W + -sin (®12(К)т) [ b - e' — X v l/2(^j ^ V c c " 1/2( к ) co' ( 1 - cos (®у2(К)т) ) : (10) где К = 0 . Находим поля E 1(0) = E (x 1(0), t 1) , B 1(0) = B (x 1(0), t 1) . I.4. Находим при К = 0 поля e V2(k+1) = + ++ 0 + e 1к)) B 12К+1) = —+B 0 + ° 1к)) далее по ним и формулам в (8), находим b І 2 ( К + 1 ) v 1 2 ( к + 1) ^ 1 / 2 ( к + 1 ) / 7 У 2 ( к + 1) Е E l і / 2 ( к + і ) = ( v 0 . b 1 2 ( K + 1 ) ) b 1 2 ( K + 1 ) v \ і / 2( к + і ) _ , ,0 і/2К + і ) _ , 1/2(К+1) = v - v - v E . Затем по формуле (10) при к = 1 находим x 1(1j= x ( t 1, t 0, x 0, v0) в постоянных полях E 12(1), B ^ 2(1). Далее находим поля E 1(1) = E (x 1(1), t 1) , B 11 = B (x 11, t1) , а также параметр точности 8 1 = x 1(1)- x 1(0) . Если достигнута заданная точность 8 : Е > L c v 147