6
b1= 0.04833419895509438083025798424226646350104-0.3205959705202491275579925812180370610420i
b2= 0.9516658010449056191697420157577335364990-1.696774337833590152863879033397701636829i
c= 0.1730887968652108268236509028841256629989-0.1694095699539015207491928002966606686529i
a
= 0.5359744564304914728722121662633122018534-0.9665922748484189311614005696233999104339i
Коэффициенты для непрерывного метода 3-го порядка
p11=3.619032671469712553690516593029702603963-0.2004225284051934625395768416824377885643i
p12=-8.414446230029770923182842133180291069678-0.6613406243058764343155995103891857483354i
p13=4.843747757515152750322583524392854929215+0.5411671821908207692971837708535864758580i
p21=-2.619032671469712553690516593029702603963-9.183665759393742328007731132919548994127i
p22=8.414446230029770923182842133180291069678+29.50534469034316892770669484128541508572i
p23=-4.843747757515152750322583524392854929215-22.01845326878301675256284274176356772842i
2.3. CROS2_Optimal
Данный метод взят из [4]. Это двухстадийный метод с комплексными коэффициентами 4-го порядка. Он
имеет оптимальное L – затухание как для функции устойчивости, так и для функции внутренней
устойчивости. Метод задается следующими формулами:
2
2
1
1
1
Re
k
k
y
y
n
n
)
(
)
(
1
1
n
n
y
y
f
k
y
f
E
1
21
1
21
1
21
2
2
))
Re(
(
))
Re(
(
)
(
k
k
y
f
k
y
f
k
y
f
E
n
y
n
n
y
Уточненные параметры даны в таблице 3. Там так же приводится непрерывный метод, основанный на
формуле (7).
Таблица 3.
Параметры CROS2_Optimal
α1=0.4573733434972975729886712378597810560027-0.2351004879985426732069962332538443381790i
α2=0.04262665650270242701132876214021894399728-0.3946329531721133798109289632446957436482i
γ21=0.9147466869945951459773424757195621120054+0.6546908396281089262421198339615985775839i
δ21=0.7451903446015528669012176569684448709433+3.042955553310975129001697189423158388626i
π21=-0.2317220370463119616442110173716549151421+0.07081125027258135212790263920701540770149i
ω1=0.2784065608064581775941260630978839609622-0.9196225438536243933605777235086716993250i
ω2=0.7215934391935418224058739369021160390378+0.2014773772275646399542954916451697899548i
Коэффициенты для непрерывного метода 3-го порядка
p11=1.662352335632077153904266284657282585433+3.058378731739278701194655437881970784617i
p12=-1.145501425793116439658894990826412127076-6.527071235918501156496073309360094210372i
p13=-0.2384443490325025366512452307329864973956+2.549069960325598061940840147969451726429i
p21=-0.6623523356320771539042662846572825854333-3.084120463545428936697275176270423611394i
p22=1.145501425793116439658894990826412127076+5.333814343596499665083881944741250473414i
p23=0.2384443490325025366512452307329864973956-2.048216502823506088432311276825657072065i