CC BY
100
Известия высших учебных заведений. Поволжский регион
УДК 519.622
E. A. Новиков
АЛГОРИТМ ПЕРЕМЕННОЙ структуры С ПРИМЕНЕНИЕМ ТРЕХСТАДИЙНЫХ МЕТОДОВ ТИПА РУНГЕ - КУТТЫ И РОЗЕНБРОКА1
Аннотация.
Актуальность и цели. Проблема решения задачи Коши для жестких систем большой размерности возникает при моделировании физических и химических процессов, при аппроксимации уравнений в частных производных системой обыкновенных дифференциальных уравнений и во многих других важных приложениях. Учет большого числа факторов при построении математических моделей приводит к расширению класса задач, описываемых жесткими системами большой размерности. Сложность практических задач приводит к возрастающим требованиям к вычислительным алгоритмам.
Материалы и методы. В случае большой размерности жесткой системы дифференциальных уравнений основные затраты приходятся на декомпозицию матрицы Якоби. В некоторых алгоритмах применяется замораживание матрицы Якоби, т.е. одна матрица используется на нескольких шагах интегрирования. Проблема замораживания матрицы достаточно просто решается в методах, в которых стадии вычисляются с применением итерационного процесса. Для безытерационных численных формул это существенная проблема.
В данной работе сокращение затрат достигается за счет комбинирования явных и Z-устойчивых методов по критерию устойчивости в процессе расчетов.
Результаты. Создан алгоритм интегрирования переменной структуры на основе явной схемы типа Рунге - Кутты и Z-устойчивого метода типа Розен-брока третьего порядка. На каждом шаге эффективная численная формула выбирается по критерию устойчивости. Оценка максимального собственного числа матрицы Якоби, необходимая для переключения между методами, для явных численных схем определяется степенным методом через ранее вычисленные стадии, а для метода типа Розенброка - через норму матрицы Якоби. Построены неравенства для контроля точности и устойчивости. Приведены результаты расчетов.
Выводы. Алгоритм интегрирования предназначен для решения жестких задач большой размерности. Результаты расчетов подтверждают эффективность построенного алгоритма.
Ключевые слова: жесткая система, схемы типа Рунге - Кутты и Розенбро-ка, контроль точности и устойчивости, автоматический выбор метода.
E. A. Novikov
VARIABLE STRUCTURE ALGORITHM APPLYING 3-STAGE METHODS OF RUNGE-KUTTA AND ROSENBROCK TYPES
Abstract.
Background. The Cauchy problem for large-scale stiff systems arises in simulation of physical and chemical processes, in approximation of partial differential equations by a system of ordinary differential equations and in plenty of other im-
1Работа выполнена при поддержке Российского фонда фундаментальных исследований (проект № 14-01-00047).
50
University proceedings. Volga region
№ 3 (35), 2015
Физико-математические науки. Математика
portant applications. Taking into consideration a large number of factors in model development leads to expansion of a class of problems, determined by stiff systems of high dimension. The complexity of practical problems leads to additional requirements to computational algorithms.
Materials and methods. Dealing with high dimensionality of a stiff system of ordinary differential equations, the main computational expenses concern the Jacobi matrix decomposition. In some algorithms one may use freezing of the Jacobi matrix, i.e. applying the same matrix over several integration steps. The problem of freezing is solved rather easy in those methods, the stages of which are computed using some iterative processes. For non-iterative numerical formulas the freezing is quite a difficult problem. In this study, the costs reduction was achieved by combining explicit and L-stable methods using the criterion of stability in calculations.
Results. The author has created an algorithm of variable structure integration, based on the explicit scheme of the Runge-Kutta type and the L-stable method of the Rosenbrock type. Both schemes have the third order of accuracy. An efficient numerical formula was chosen according to the criterion of stability at each step of integration. Estimation of maximum eigen value, which is necessary to switch between the methods, for explicit numerical schemes was determined by power iterations using already computed stages, and using the Jacobi matrix norm for the Ros-enbrock type method. The researcher also formulated inequalities for accuracy and stability control. The article adduces the results of calculations.
Conclusions. The integration algorithm is aimed at solving stiff problems of high dimension. Numerical results confirm the efficiency of the constructed algorithm.
Key words: stiff system, schemes of Runge-Kutta and Rosenbrock types, accuracy and stability control, automatic selection of a method.
Введение
Проблема решения задачи Коши для жестких систем обыкновенных дифференциальных уравнений большой размерности возникает во многих важных прикладных задачах [1-3]. Для численного решения жестких задач обычно применяются Z-устойчивые методы [1]. При реализации таких численных схем на каждом шаге решается линейная система алгебраических уравнений с применением LU-разложения некоторой матрицы, размерность которой совпадает с размерностью вектора решения. В силу плохой обусловленности матрицы Якоби жесткой системы дифференциальных уравнений, решение алгебраической системы осуществляется с выбором главного элемента по строке или столбцу, а иногда и по всей матрице. При большой размерности исходной задачи общие вычислительные затраты фактически полностью определяются временем декомпозиции данной матрицы.
В последнее время при решении жестких задач широкое распространение получили методы типа Розенброка [4], которые относятся к одношаговым безытерационным численным формулам. В отличие от неявных или полуявных методов типа Рунге - Кутты, в данных численных схемах матрица Якоби введена непосредственно в вычислительную формулу. В результате вместо решения нелинейных систем алгебраических уравнений на каждом шаге несколько раз решается линейная система. Введение матрицы Якоби в численную схему является и недостатком, потому что возникают принципиальные проблемы с замораживанием матрицы Якоби [3]. Известно, что максимальный порядок точности численных формул с замораживанием равен двум [5]. Это означает, что применение методов типа Розенброка ограни-
Physical and mathematical sciences. Mathematics
51
Известия высших учебных заведений. Поволжский регион
чено либо задачами небольшой размерности, либо расчетами с небольшой точностью [6].
Некоторым аналогом замораживания матрицы Якоби является применение в расчетах алгоритмов интегрирования на основе явных и Z-устойчивых методов с автоматическим выбором численной схемы [7]. В этом случае эффективность алгоритма может быть повышена за счет расчета переходного участка, соответствующего максимальному собственному числу матрицы Якоби, явным методом. В качестве критерия выбора эффективной численной формулы естественно применять неравенство для контроля устойчивости [8]. Применение таких комбинированных алгоритмов полностью не снимает проблему замораживания матрицы Якоби, потому что явным методом можно просчитать, вообще говоря, только погранслойное решение, соответствующее максимальному собственному числу матрицы Якоби.
Здесь на основе явного метода типа Рунге - Кутты и Z-устойчивой схемы типа Розенброка третьего порядка точности построен явно-неявный алгоритм переменной конфигурации. Выбор эффективной численной формулы на каждом шаге осуществляется по критерию устойчивости. Приведены результаты расчетов, подтверждающие эффективность алгоритма интегрирования переменной структуры.
1. Метод Рунге - Кутта - Фельберга третьего порядка
Для численного решения задачи Коши
у, = f ( У К У (о ) = Уo, (о < t < tk, (1)
будем применять явный трехстадийный метод следующего вида: yn+1 = yn + Pp1k1 + Р p2k2 + pp3k3,
k1 = hf (n, Уп ), k2 = hf (tn + К Уп + k1 ),
кз = hf
1 , 1 , 1 ,
fn + 2 h, yn + 4 k1 + 4 k2
(2)
При значениях параметров
P31 = P32 = 1/6, P33 = 2/3, (3)
формула (2) совпадает с методом Фельберга третьего порядка [9], а при коэффициентах
Р21 = Р22 = °.5, Р23 = 0,
формула (2) имеет второй порядок точности. Введем обозначения
(4)
£n,3
S (P3i - P2i)
i=1
= 1ll2k3
где Ц-Ц - некоторая норма в RN .
Тогда для контроля точности схемы (2) с коэффициентами (3) можно применять неравенство £n 3 < е, где £ - требуемая точность интегрирования.
52
University proceedings. Volga region
№ 3 (35), 2015
Физико-математические науки. Математика
2. Контроль устойчивости
Известно, что метод (2), (3) при решении нежестких задач является эффективным. Однако его применение для расчета жестких задач на участке установления приводит к большому числу повторных вычислений решения (возвратов) за счет возникающей неустойчивости. Этих возвратов частично можно избежать за счет дополнительного контроля устойчивости численной схемы [8]. Построим неравенство для контроля устойчивости. Для этого применим метод (2), (3) для решения тестовой задачи Дальквиста [10]:
у =Ху, Re(X) < 0 , у(0) = уо, t > 0 . (5)
В результате получим уп+1 = 2з(x)yn, где x = h - X, а многочлен устойчивости 2з(x) имеет вид Q4(x) = 1 + x + x /2 + x /6. На рис. 1 приведены линии уровня |^з(x)| = q при значении q, равном 1; 0,7 и 0,3. Многочлен
устойчивости Q2(x) метода (2), (4) второго порядка точности имеет вид 2
Q2( x) = 1 + x + 0, 5x , а его область устойчивости приведена на рис. 2.
Рис. 1. Область устойчивости метода третьего порядка
Для построения неравенства для контроля устойчивости метода (2), (3) запишем стадии £1, £2 и £3 применительно к задаче у = Ay, где A есть матрица с постоянными коэффициентами. Получим
£1 = Xyn , £2 =( X + X2) уп , £3 =( X + 0,5 X2 + 0,25 X3) Уп, где X = hA.
Нетрудно видеть, что выполнены соотношения
£2 - £1 = X2 уп , 2 - (2£3 - £2 - £1) = X3 уп .
Physical and mathematical sciences. Mathematics
53
Известия высших учебных заведений. Поволжский регион
Рис. 2. Область устойчивости метода второго порядка
Тогда оценку максимального собственного числа vn = h A,n max| матрицы Якоби задачи (1) можно вычислить степенным методом по приближенной формуле
vn = 2 ■ max
1<i < N
2kl3 - k2 - k{
k2
(6)
В результате для контроля устойчивости схемы (2) можно применять неравенство vn < D, где константу D можно выбрать равной 2,5, т.е. примерно равной длине интервала устойчивости схемы (2), (3) (см. рис. 1).
Оценка (6) максимального собственного числа матрицы Якоби задачи (1) является грубой, потому что вовсе не обязательно максимальное собственное число сильно отделено от остальных, в степенном методе применяется мало итераций и дополнительные искажения вносит нелинейность задачи (1). Поэтому здесь контроль устойчивости используется как ограничитель на размер шага интегрирования. В результате прогнозируемый шаг hn+1 вычисляется следующим образом. Новый шаг hac по точности определяется по формуле hac = q^hn, где hn есть последний успешный шаг интегрирования, а
q1, учитывая соотношение en 3 = O (h3), задается уравнением q3 ■ en 3 =£ .
Шаг hst по устойчивости задается формулой hst = q2hn , где q2, учитывая соотношение vn = O (hn), определяется из равенства q2vn = 2,5 . В результате прогнозируемый шаг hn+1 вычисляется по следующей формуле:
hn+1 = max
hn ,min (hac, hst)
(7)
54
University proceedings. Volga region
№ 3 (35), 2015
Физико-математические науки. Математика
Если шаг по устойчивости должен быть уменьшен, то в качестве следующего шага будет применяться последний успешный шаг hn . Данная формула позволяет стабилизировать поведение шага на участке установления решения, где определяющую роль играет устойчивость. Из результатов расчетов алгоритмом с дополнительным контролем устойчивости следует повышение эффективности примерно в полтора-два раза, однако точности вычислений существенно выше задаваемой. Это естественно, потому что старые ошибки подавляются за счет контроля устойчивости, а новые ошибки на участке установления невелики.
3. Метод типа Розенброка
Для решения задачи (1) рассмотрим численную формулу вида Уп+1 = Уп + Ml + P2k2 + P3k3 ,
Dnk1 = hf (tn,Уп ) Dnk2 = hf (n +в21h,Уп + p21k1)
Dnk3 = hf (n + [( +Р32 ]h, yn + ^31k1 +Р32 k2 ) , (8)
где h - шаг интегрирования; Dn = E - ahfn, E - единичная матрица;
fn = df (tn,yn)/ Эу - матрица Якоби системы (1); a , pt и - числовые ко-
эффициенты. Разложим стадии kt, 1 < i < 3 , по степеням h и подставим в первую формулу (8). Полагая уп = у (tn) и сравнивая ряды для точного и
приближенного решений до членов с h включительно, получим условия третьего порядка точности схемы (8), которые после несложных преобразований записываются так:
P1 + P2 + P3 = l P21P2 +PP3 = ■“ - a, (9)
P2bp2 +Р2P3 = 3 Р21Р32P3 = Т-a + a2, (10)
3 6
где Р = Р31 +Р32.
Локальная ошибка 5п 3 схемы (8) при условии (10) имеет вид
n,3 = I 1 3 a 3 2 1 + — a
{ 24 22
" 1 1 (1
-a -p21 I — a+a
_ 12 3 21 16
h4 f ,3f +
11 I 1 2 i„
---a-I---a + a I-B
8 3 16 1
h 4 fff 2 +
h4 fff2 +
-4-P21 p2 -P3P3
h4f f3 + O (h5 ) .
4. Исследование устойчивости
Применяя численную формулу (8) для решения тестовой задачи (5), имеем yn+1 = Qr (x)yn, где x = h -X, а функция устойчивости Qr (x) имеет вид
Physical and mathematical sciences. Mathematics
55
Известия высших учебных заведений. Поволжский регион
Qr (x) =-
1 + (1-За)х + За2 -2а + Р21 Pi + (31 +Рзз)Р3
(1 - ах)
{а3 -а2р1 + а(21 -а)р2 - а2 -а(31 +^32) + Р21Р32 рз}х
(1 - ах)
Необходимым условием L -устойчивости является требование того, чтобы степень многочлена в числителе была на единицу меньше степени полинома в знаменателе. Нетрудно видеть, что это требование будет выполне-3 2
но, если а - 3а + 1,5а -1/6 = 0. При записи условия L -устойчивости учтены соотношения (10). Известно [11], что схема (8) будет A -устойчивой, если параметр а удовлетворяет неравенству 1/3 < а < 1,0685790. Поэтому среди корней уравнения выбираем а = 0,435866521508459, при котором численная схема (8) будет L -устойчивой.
При расчетах по схеме (8) в двух промежуточных точках tn + вц^ и tn +(31 +Р31) вычисляются приближенные решения yna= yn +Р2А и yn р = yn +^31^1 +^32^2 . Применяя их для решения тестового уравнения (5), имеем yn а = Qa(х)yn и yn р = Qp (х)yn , где функции устойчивости внутренних или промежуточных численных схем имеют вид
Q а (х)
1 + (21 - а )х 1 - ах
• Уп ,
Q р(х ) =
1 + (31 +в32-2а )х + а2-аР + Р21Р32 (1 - ах )2
"Уп .
Для уменьшения влияния ошибок в промежуточных вычислениях потребуем L -устойчивости промежуточных численных формул, что приводят
к двум соотношениям а - ав + в21в32 = 0 и Р21 - а = 0 . Теперь с учетом равенства Р21 = а условия порядка (10) и устойчивости записываются в виде
P1 + P2 + P3 = 1, ар2 +вР3 = 0,5(1 - 2а), а2 P2 +р2 P3 =1/3,
Р32P3 = а-1 (6а2 -6а +1)/6, Р32 =Р-а , Р = Р31 +Р32 .
Исследуя совместность этой нелинейной алгебраической системы, получим коэффициенты L -устойчивого метода (8) третьего порядка:
а = 0.435866521508459, р1 =р21 =р31 = а , в = а (6а2 - 6а +1) (6а2 - 3а + 2), Р32 = в - а ,
P3 = а 1 (в-а) 1 (6а2 -6а +1)/6,
56
University proceedings. Volga region
№ 3 (35), 2015
Физико-математические науки. Математика
Р2 = 0.5а 1 (1 - 2а - 2(3рз), р = 1- р2 - ръ.
5. Контроль точности вычислений
Контроль точности вычислений построим по аналогии [12]. Для оценки ошибки используем вложенный метод второго порядка
Уп+1,2 = Уп + b1k1 + b2k2 , Dnk1 = hf (Уп ) , Dnk2 = hf (Уn + ak1),
где приближение yn вычислено по формуле (8), b = 0,5а 1 (4а -1) и b2 = 0,5а-1 (1 — 2а). В результате оценка ошибки вычисляется по формуле
еп (Jn ) Dn " ( уп+1 + Уп+1,2 ). Теперь неравенство для контроля точности имеет вид
D
1 Jn
(yn+1 + уп+1,1 ) ) , 1 - jn < 2,
где К| - некоторая норма в RN; е - требуемая точность интегрирования, а параметр jn выбирается наименьшим значением, при котором выполняется неравенство.
6. Алгоритм переменной структуры
Оценку максимального собственного числа vn о = h • |^п max | матрицы
Якоби системы (1), необходимую для перехода на явную формулу, оценим через ее норму по формуле
vn,0 =
df (Уп )
дУ
= max {
1-i - N
N
z
j=1
dfi (Уп )
дУ,
> .
Норма d в неравенствах для контроля точности вычисляется по фор-
муле
IId 11= max fcl /
1<i - N l
+ r
где r - положительный параметр. В случае выполнения неравенства
уП
< r
по i -й компоненте решения контролируется абсолютная ошибка r е . В противном случае контролируется относительная ошибка е .
На основе построенных методов легко сформулировать алгоритм, в котором на каждом шаге выбирается наиболее эффективная численная схема. Расчеты всегда начинаются явным методом, потому что в нем не используется матрица Якоби. Нарушение неравенства vn < 2,5 вызывает переход на L -устойчивую схему, где оценка vn вычисляется по формуле (6). Передача
Physical and mathematical sciences. Mathematics
57
Известия высших учебных заведений. Поволжский регион
управления явным методам происходит в случае выполнения неравенства vn о < 2,5, где оценка vn 0 вычисляется через норму матрицы Якоби.
7. Результаты расчетов
Вычисления проводились с двойной точностью, потому что для жестких задач характерен большой разброс входящих в описание задачи констант. В расчетах параметр r выбирался таким образом, чтобы по всем компонентам решения фактическая точность была не хуже задаваемой. Расчеты проводились с точностью 8 = 10-4 , при которой наиболее эффективны методы третьего порядка. Ниже через if и id обозначены соответственно суммарное
число вычислений правой части и количество декомпозиций матрицы Якоби задачи (1), которые позволяют объективно оценить эффективность алгоритма интегрирования. В качестве тестового примера выбрана простейшая математическая модель с периодическим решением для описания реакции Белоусова -Жаботинского (орегонатор) [13]:
у1 = 77,27(У2 -yiУ2 + Л -8,375 • 10-6у2),
У2 = 7727 (2 -У1У2 + УзУ У3 = 0,161(У1 -УзУ
tе [0,300], У1 (0) = Уз(0) = 4, У2(0) = 1,1, h = 10-3.
Коэффициент жесткости равен примерно 10 . Для построенного алгоритма переменной структуры вычислительные затраты следующие: if = 3 983
и id = 400, а при расчетах только по численной схеме Розенброка следующие: if = 3 179 и id = 706. Из сравнения затрат следует сокращение числа
декомпозиций матрицы Якоби за счет комбинирования методов почти в два раза. Данную задачу можно посчитать явным методом Рунге - Кутта. При расчетах алгоритмом без контроля устойчивости if = 11 011 774, при расчетах с дополнительным контролем устойчивости if = 8 920 580 .
Заключение
В построенном алгоритме с помощью признака можно задавать различные режимы расчета:
1) явным методом с контролем или без контроля устойчивости;
2) Z-устойчивым методом с аналитической или численной матрицей Якоби;
3) с автоматическим выбором численной схемы.
Все это позволяет применять данный алгоритм для решения как жестких, так и нежестких задач. При расчетах с автоматическим выбором численной схемы вопрос о том, является задача жесткой или нет, перекладывается на алгоритм интегрирования. Из результатов численных расчетов следует существенное сокращение числа декомпозиций матрицы Якоби за счет ком-
58
University proceedings. Volga region
№ 3 (35), 2015
Физико-математические науки. Математика
бинирования методов, что при решении систем дифференциальных уравнений большой размерности сокращает время вычислений почти в два раза.
Список литературы
1. Hairer, E. Solving Ordinary Differential Equations I: Nonstiff and Differential-Algebraic Problems / E. Hairer, S. P. Norsett, G. Wanner. - Berlin : Springer-Verlag, -1993. - 528 p.
2. Hairer, E. Solving Ordinary Differential Equations II: Stiff and Differential-Algebraic Problems / E. Hairer, G. Wanner. - Berlin : Springer-Verlag, 1996. - 614 p.
3. Новиков, E. А. Компьютерное моделирование жестких гибридных систем / Е. А. Новиков, Ю. В. Шорников. - Новосибирск : Изд-во НГТУ, 2012. - 451 с.
4. Rosenbrock, H. H. Some general implicit processes for the numerical solution of differential equations / H. H. Rosenbrock // Computer. - 1963. - № 5. - P. 329-330.
5. Новиков, E. A. Замораживание матрицы Якоби в методе типа Розенброка второго порядка точности / В. А. Новиков, Е. А. Новиков, Л. А. Юматова // Журнал вычислительной математики и математической физики. - 1987. - Т. 27, № 3. -С. 385-390.
6. Novikov, E. A. Numerical Integration of Stiff Systems with Low Accuracy / A. E. Novikov, E. A. Novikov // Mathematical Models and Computer Simulations. -2010. - Vol. 2, № 4. - P. 443-452.
7. Novikov, E. A. Construction of algorithm for the integrating stiff differential equations on no uniform schemes / E. A. Novikov // Soviet Math. Dokl. - 1984. - Vol. 30, № 2. - P. 358-361.
8. Новиков, Е. А. Явные методы для жестких систем / Е. А. Новиков. - Новосибирск : Наука, 1997. - 198 с.
9. Fehlberg, E. Low order classical Runge-Kutta formulas with stepsize control and their application to some heat transfer problems / E. Fehlberg. - Washington : NASA T.R.R., 1969. - 316 p.
10. Dahlquist, G. A special stability problem for linear multistep methods / G. Dahlquist // BIT. - 1963. - m 3. - P. 23-43.
11. Демидов, Г. В. Исследование точности неявных одношаговых методов / Г. В. Демидов, Л. А. Юматова. - Новосибирск : Изд-во ВЦ СО АН СССР, 1976. -(Препринт № 25).
12. Демидов, Г. В. Оценка ошибки одношаговых методов интегрирования обыкновенных дифференциальных уравнений / Г. В. Демидов, Е. А. Новиков // Численные методы механики сплошной среды. - 1985. - Т. 16, № 1. - С. 27-39.
13. Enright, W. H. Comparing numerical methods for the solutions of systems of ODE’s / W. H. Enright, T. E. Hull // BIT. - 1975. - Vol. 15. - P. 10-48.
References
1. Hairer E., Norsett S. P., Wanner G. Solving Ordinary Differential Equations I: Nonstiff and Differential-Algebraic Problems. Berlin: Springer-Verlag, 1993, 528 p.
2. Hairer E., Wanner G. Solving Ordinary Differential Equations II: Stiff and Differential-Algebraic Problems. Berlin: Springer-Verlag, 1996, 614 p.
3. Novikov E. A., Shornikov Yu. V. Komp’yuternoe modelirovanie zhestkikh gibridnykh sistem [Computer simulation of stiff hybrid systems]. Novosibirsk: Izd-vo NGTU, 2012, 451 p.
4. Rosenbrock H. H. Computer. 1963, no. 5, pp. 329-330.
5. Novikov E. A., Novikov V. A., Yumatova L. A. Zhurnal vychislitel’noy matematikt i matematicheskoy fikhiki [Journal of calculus mathematics and mathematical physics]. 1987, vol. 27, no. 3, pp. 385-390.
Physical and mathematical sciences. Mathematics
59
Известия высших учебных заведений. Поволжский регион
6. Novikov E. A., Novikov A. E. Mathematical Models and Computer Simulations. 2010, vol. 2, no. 4, pp. 443-452.
7. Novikov E. A. Soviet reports on mathematics. 1984, vol. 30, no. 2, pp. 358-361.
8. Novikov E. A. Yavnye metody dlya zhestkikh sistem [Explicit methods for stiff systems]. Novosibirsk: Nauka, 1997, 198 p.
9. Fehlberg E. Low order classical Runge-Kutta formulas with stepsize control and their application to some heat transfer problems. Washington: NASA T.R.R., 1969, 316 p.
10. Dahlquist G. A special stability problem for linear multistep methods. BIT. 1963, vol. 3, pp. 23-43.
11. Demidov G. V., Yumatova L. A. Issledovanie tochnosti neyavnykh odnoshagovykh metodov [Accuracy research for non-explicit one-step methods]. Novosibirsk: Izd-vo VTs SO AN SSSR, 1976, (Preprint № 25).
12. Demidov G. V., Novikov E. A. Chislennye metody mekhaniki sploshnoy sredy [Numerical methods of continuum mechanics]. 1985, vol. 16, no. 1, pp. 27-39.
13. Enright W. H., Hull E. Comparing numerical methods for the solutions of systems of ODE’s. BIT. 1975, vol. 15, pp. 10-48.
Новиков Евгений Александрович
доктор физико-математических наук, профессор, главный научный сотрудник, Институт вычислительного моделирования Сибирского отделения РАН (Россия, г. Красноярск, Академгородок, 50, стр. 44)
Novikov Evgeniy Aleksandrovich Doctor of physical and mathematical sciences, professor, chief researcher, Institute of computational modeling of the Siberian branch of RAS (building 44, 50 Akademgo rodok street, Krasnoyarsk, Russia)
E-mail: [email protected]
УДК 519.622 Новиков, Е. А.
Алгоритм переменной структуры с применением трехстадийных методов типа Рунге - Кутты и Розенброка / E. A. Новиков // Известия высших учебных заведений. Поволжский регион. Физико-математические науки. -2015. - № 3 (35). - С. 50-60.
60
University proceedings. Volga region
