Построение приближенного решения одного нелинейного дифференциального уравнения второго порядка в окрестности подвижной особой точки Текст научной статьи по специальности «Математика»

УДК 517.928.4

ПОСТРОЕНИЕ ПРИБЛИЖЕННОГО РЕШЕНИЯ ОДНОГО НЕЛИНЕЙНОГО ДИФФЕРЕНЦИАЛЬНОГО УРАВНЕНИЯ ВТОРОГО ПОРЯДКА В ОКРЕСТНОСТИ ПОДВИЖНОЙ ОСОБОЙ ТОЧКИ

В.Н. Орлов, Т.Ю. Леонтьева

Чувашский государственный педагогический университет им. И.Я.Яковлева,

Чебоксары, Российская Федерация

e-mail: orlowvn@rambler.ru; betty2784@yandex.ru

Дифференциальные уравнения представляют собой математические модели различных процессов и явлений. В отличие от линейных, нелинейные дифференциальные уравнения мало изучены. Сложности в изучении нелинейных дифференциальных уравнений связаны с наличием подвижных особых точек. Предложен метод приближенного решения нелинейного дифференциального уравнения с подвижными особыми точками, включающий в себя решение шести задач. Решены первые две задачи: теорема существования и единственности решения нелинейного дифференциального уравнения; построение приближенного решения и исследование влияния возмущения начальных условий на приближенное решение. Приведено доказательство теоремы существования и единственности решения рассматриваемого нелинейного дифференциального уравнения в окрестности подвижной особой точки и получена структура приближенного решения.

Ключевые слова: подвижная особая точка, нелинейное дифференциальное уравнение второго порядка, приближенное решение, метод мажорант, окрестность подвижной особой точки.

CONSTRUCTION OF APPROXIMATE SOLUTION OF ONE NONLINEAR DIFFERENTIAL SECOND-ORDER EQUATION IN THE NEIGHBORHOOD OF MOVABLE SINGULAR POINT

V.N. Orlov, T.Yu. Leont'eva

I. Yakovlev Chuvash State Pedagogical University, Cheboksary, Russian Federation

e-mail: orlowvn@rambler.ru; betty2784@yandex.ru

Differential equations are mathematical models of various processes and phenomena. In contrast to linear, the nonlinear differential equations have been poorly studied. Difficulties in the study of nonlinear differential equations are associated with movable singular points. The method to solve approximately a nonlinear differential equation with movable singularities comprising the decision of six problems, is proposed. First two problems — theorem of existence and uniqueness of solutions of nonlinear differential equations as well as construction of approximate solutions and study of the initial condition perturbation influence on the approximate solution, are considered. Proof of the theorem of the solution existence and uniqueness of this nonlinear differential equation in the neighbourhood ofmovable singular point, and the approximate solution structure, are given.

Keywords: movable singular point, nonlinear differential equation of second order, approximate solution, method of majorants, neighbourhood of movable singular point.

Введение. Нелинейные дифференциальные уравнения представляют большой интерес в связи с их приложением во многих областях

науки и техники [1]. Решения нелинейных уравнений связаны с большими трудностями, вызванными наличием подвижных особых точек у интегралов этих уравнений, которые и являются препятствием к использованию известных в настоящее время приближенных численных и аналитических методов решения [2, 3]. Отметим работы, направленные на разрешение в квадратурах нелинейных дифференциальных уравнений с подвижными особыми точками, но это удается сделать лишь в частных случаях [4].

Материалы и методы решения задачи и принятые допущения. Применяется метод решения нелинейных дифференциальных уравнений, предложенный в работах [1, 5], который состоит из шести этапов. В статьях [6-9] были решены первые две задачи указанного выше метода для области аналитичности.

В настоящей работе предложено решение следующих задач: доказана теорема существования и единственности решения рассматриваемого нелинейного дифференциального уравнения в окрестности подвижной особой точки и построено приближенное решение в указанной области. При доказательстве теоремы существования был применен метод мажорант, но не к правой части дифференциального уравнения, как это сделано в классической литературе [10], а к решению самого нелинейного дифференциального уравнения [1, 5, 11].

В статье использовано понятие "нормальная форма" для нелинейных дифференциальных уравнений, введенное в работе [1], когда общая структура уравнения с полиномиальной правой частью у(к) (х) = а0 (х) уп (х) + а1 (х) уп-1 (х) + ... + ап-1 (х) у (х) + ап (х) с помощью некоторой замены приводится к виду у(к) (х) = уп (х) + + г (х), где к — порядок производной; п — степень функции. Эта терминология является продолжением для упрощенной структуры нелинейного дифференциального уравнения, упомянутой для уравнения Абеля в работе [12].

Результаты. Рассматривается нелинейное дифференциальное уравнение второго порядка

у" (х) = а0 (х) у5 (х) + а1 (х) у4 (х) + а2 (х) у3 (х) +

+ аз (х) у2 (х) + а4 (х) у (х) + а5 (х), (1) где а», % = 0,1,..., 5 — аналитические функции в рассматриваемой

С п " " Г \ и (х) а1 (х)

области. Заменой переменной у (х) = —. — -—— при выползу ао (х) 5ао (х)

нении условий

а0' (х) 5 /а0' (х)

а4 (x) +

a (x) а,2 (x) аэ (x) \ 4а°(x) 16 (x)

5а° (x) 2ai (x) а2 (x) аэ (x)

уравнение (1) приводится к нормальной форме [6] u" (ж) = u5 (ж) + + r (ж). Здесь

a,5 (ж) ^ao (ж) 3a0' (ж) ^ao (ж) 4/—--

r (ж) =--K5 4,4---on ^-+ a, (ж) V ao (ж) +

55a0(ж) 20a2 (ж)

a'/ (ж) a0 (ж) 2a0 (ж) a/ (ж) a0 (ж)

5a0(ж) 5a2(ж)

+ 2 (a0 (ж))2 a/ (ж) ^a0 (ж) (aQ (ж))2 a/ (ж) ^ap (ж) + a0 (ж) ^ (ж)

5a0 (ж) 16a0 (ж) 2a0 (ж) '

в каждой области, в которой a0 (ж) = 0. Рассмотрим задачу Коши

у" (ж) = у5 (ж) + r (ж); (2)

У (ж0) = У0,у' (ж0) = yi. (3)

Теорема 1. Пусть ж* — подвижная особая точка у (ж) задачи (2), (3) и функция r (ж) удовлетворяет следующим условиям:

1) r (ж) е C^ в области |ж* — ж| < р0, р0 = const > 0;

|r(n) (ж*) I

2) ЗМ0 : J-^^ < M0, M0 = const, n = 0,1, 2,...

n!

Тогда существует единственное решение задачи Коши (2), (3) в виде

ж

У (ж) = (ж* — ж)-1/2 £ Cn (ж* — ж)п/2 , C0 = 0, (4)

n=0

правильная часть которого сходится в области

|ж* — ж| < Р2, (5)

где р2 = min {р0,р1}, Р' = —, , M = max {М0,а}, а — пара-

4 У (+1)2 метр, зависящий от условий (3).

< В соответствии с условием теоремы имеем

r (ж) = £ Bn (ж* — ж)п . (6)

n=0

Учитывая структуру решения в окрестности подвижной особой точ-

„* _ „\Р ^^ П („* _ ^W2

n=0

ки, в общем случае у (ж) = (ж* — ж)р ^^ Cn (ж* — ж)п/2, C0 = 0, из

уравнения (2) получаем

£ Cn (ж* — ж)п/2+А = ( £ Cn (ж* — ж)п/2+р ) +£ Bn (ж* — ж)п

n=0 n=0 n=0

После выполнения соответствующих операций имеем

те

^ С (п/2 + р) (п/2 + р - 1) (X - х)п/2+р-2 =

п=0

те те

= ^ Бп* (х* - х)п/2+5р + ^ Бп (х* - х)п , (7)

п=0 п=0

п

п п

где Б** = V Б*-гСг, Б* = £ Бп-гБг, Бп = £ С-^, п = 0,1, 2,...

г=0 *=0 *=0

Равенство (7) обратится в тождество при выполнении условий

1) п/2 + р - 2 = п/2 + 5р; (8)

2) Сп (п/2 + р) (п/2 + р - 1) = Б**. (9)

Здесь

, Б**,п £ {0,1, 2, 3, 4, 6, 8,..., 2k} ; Б i Б** + Bn-5,n £{5, 7, 9,..., 2k + 1} . ( )

Из условия (8) следует р = -1/2, а соотношение (9) позволяет однозначно определить все коэффициенты Сп. При этом вариант С0 = 0 отбрасываем, так как в этом случае получаем противоречие с предположением относительно характера особой точки х*. Согласно соотно-

13 4

шению (9), С0 = ^3, С1 = С = Сз = С4 = 0, С5 = -4Б0, Се = а,

4 4 20 /27

С7 = -Бь С8 = 0, С9 = — Б2, С10 = -\ — Б2. Таким образом,

7 9 ь 8 ' 9 33 10 147 V 4 0 Р

получено однозначно формальное представление решения задачи (2),

(3) в окрестности точки х* в виде (4).

Докажем сходимость правильной части ряда (4) в области (5). Из

|г(п) (х*)|

условия теоремы следует М = эир -1-:-1, п = 0,1, 2,..., тогда для

п п!

коэффициентов ряда (6) справедлива оценка

|Бп |< М. (11) Из выражения (9) с учетом (10), (11) для коэффициентов Сп предполагаем оценку: |Сп| < 7-1-г2пМ (М + 1)[п/5] = ип, п =

(п + 2) (п - 6)

= 0, 1, 2, . . . Исходя из закономерности образования коэффициентов Сп, с помощью метода математической индукции докажем эту оценку для коэффициента Сп+1 при п + 1 = 5 (2к + 1) = 10к + 5:

|С1ок+5 ((10к + 5) /2 + р) ((10к + 5) /2 + р - 1 )| <

< |^1С1к+5 - < -^10к+5 + В10к - < |^10к+5 - + |В10к| <

<

<

10k+5

i=0

10k+5-iCi

+ |Biofc+5 | <

10k+5 / 10k+5-i

E I E D10k+5-i-j Dj i=0 \ j=0

—i—j D j I Сi

10k+5 / 10k+5-i /10k+5-i-j

I ^^ ( C10k+5-i—j—mCm

i=0 V j=0 \ m=0

ECj-iCi I Ci

ч l = 0

+ |Вюк+5 | < + |Вюк+5 | .

Тогда при p = — 1/2 после преобразований получим

|Сюк+51 <

22

(10k + 7) (10k - 1) 10k (10k-i (10k-i-j ЕЕ I Е I Е

i= 1 у j=1 у m= 1

10k-i-j-m

210k-i-j-mM (M + 1) 5 (10k — i — j — m + 2) (10k — i — j — m — 6)*

2mM (M + 1) " (m + 2) (m — 6)*

(j

ч!=1

2j-lM (M +1)^ 2lM (M +1)5 (j — l + 2) (j — l — 6)* (l + 2) (l — 6)*

x 2W(^1)* + m (i + 2) (i — 6)*

<

22

(10k + 7) (10k — 1)

10k 10k-i

ЕЕ (Е( 2

i= j=

10k-i-j j >10k-i-jm2 (m + 1) 5 2jM2 (M +1) 5 x

(10^—^ — j Е

m=1

(m + 2) (m — 6) (10k — i — j — m + 2) (10k — i — j — m — 6)

Е-

(l + 2) (l — 6)* (j — l + 2) (j — l — 6)*

2iM (M + 1) 5 (i + 2) (i — 6)*

<

210kM5 (M + 1)2k 2 I 1 +

(10k + 7) (10k — 1)

210k M4 (M + 1)

+ M) <

2k I <

2k

<

210k+2M (M + 1) (10k + 7) (10k — 1)

2k+1

22

210k+5M (M + 1) (10k + 7) (10k — 1)

где

(10k — i — j — m — 6)* =

_/ 1, (10k — г — j — m) = 1, 2, 3, 4, 5,6;

(10k — г — j — m — 6), (10k — г — j — m) = 7, 8, 9,.

1, (j — /) = 1, 2, 3,4, 5, 6, (j — l — 6), (j — 1) = 7, 8, 9,...;

1, m = 1, 2, 3, 4, 5, 6, (m — 6), m = 7, 8, 9,...;

1,/ = 1, 2, 3,4, 5, 6,

(j — / — 6)* = (m — 6)* = (/ — 6)* =

(/ — 6),/ = 7, 8, 9,

X

X

X

X

X

X

X

1

X

1

X

2

2

1

(i_ 6)* Г M = 1 2, 3, 4, 5,6, (i 6) \ (i - 6) ,i = 7, 8,9,...

Аналогичное доказательство имеем и при n +1 = 10k + 1, n + + 1 = 10k + 2, n + 1 = 10k + 3, n + 1 = 10k + 4. Рассмотрим ряд

ж

£ Vn |x* - x|(n-1)/2 , (12)

n=1

ж

Cn |x* - x|(n-1)/2.

который является мажорирующим для ряда ^ сп

п=1

Учитывая закономерность структуры коэффициентов СП, ряд (12) представляем в виде

те

* |(п-1)/2 _

У^Уп |x* - x|(

n=1

ж ж

= Е |x* - x|(5k-5)/2 + £ V5k+1 |x* - x|(5k-4)/2 +

k=1 k=1 жж

+ E V5k+2 |x* - x|(5k-3)/2 + ^ V5fe+3 |x* - x|(5k-2)/2 +

k=1 k=1

ж

+ £ V5k+4 |X* - x|(5k-1)/2

k=1

На основании признака Даламбера устанавливаем сходимость ряда (12) в области р1 = — . Полагая р2 = min {ро,р1}, получаем

4 5 (+1)2

сходимость ряда (4) в области (5), что доказывает теорему. ►

Полученные в теореме 1 оценки позволяют построить приближенное решение задачи (2), (3):

N

Vn (x) = (x* - x)-1/2£ Cn (x* - x)n/2 , Co = 0. (13)

п=0

Теорема 2. Пусть выполняются пп. 1 и 2 теоремы 1, тогда для приближенного решения (4) задачи (2), (3) в области |х* — х| < р2 справедлива оценка погрешности

Аум (х) = |у (х) — ум (х)|< А, (14)

где

25пМ (М + 1)п |х* — х|(5п-1)/2 А < ^То *

1 — 25 (М + 1) |х* — х|5/2

х ----ГГ +

2 Ix* — x|1/2

(5n + 2) (5n - 6) (5n + 3) (5n - 5)

22 |x* _ xl 23 Ix* _ xl3/2 24 Ix* _ xI2

¿j \Jb <JU \ ¿J \<Jb <JU \ ¿J \Jb <JU \

+ 73-^^-+ 77-' W '-^ +

(5n + 4) (5n - 4) (5n + 5) (5n - 3) (5n + 6) (5n - 2) Г

в случае N + 1 = 5n. Для вариантов N + 1 = 5n +1, N + 1 = 5n + 2, N + 1 = 5n + 3, N + 1 = 5n + 4 соответственно:

< 25n+1M (M + 1)n |x* - x|5n/2

< 1 - 25(M + 1) |x* - x|5/2

1 2 |x* - x|1/2 x I ^ . ^-- +

(5n + 3) (5n - 5) (5n + 4) (5n - 4)

22 |x* - x| 23 |x* - x|3/2 24 (M +1) |x* - x|2'

+ (5n + 5) (5n - 3) + (5n + 6) (5n - 2) + (5n + 7) (5n - 1)

< 25n+2M (M + 1)n |x* - x|(5n+1)/2

< 1 - 25 (M +1) |x* - x|5/2

1 2 |x* - x|1/2 x „ ....... +

(5n + 4) (5n - 4) (5n + 5) (5n - 3)

22 |x* - x| 23 (+1) |x* - x|3/2 24 (+1) |x* - x|2'

+ (5n + 6) (5n - 2) + (5n + 7) (5n - 1) + (5n + 8) 5n

< 25n+3M (M + 1)n |x* - x|(5n+2)/2

< 1 - 25 (M +1) |x* - x|5/2

( 1 2 |x* - x|1/2

X I _ _|__!_!__L

у (5n + 5) (5n - 3) (5n + 6) (5n - 2) 22 (+1) |x* - x| 23 (+1) |x* - x|3/2 + 24 (+1) |x* - x|

(5n + 7) (5n - 1) (5n + 8) 5n (5n + 9) (5n + 1)

A 2 +4M (M +1)n |x* - x|(5n+3)/2

А < -"-—-'—/2;- X

1 - 25 (M + 1) |x* - x|5/2

1 + 2 (+1) |x* - x|1/2 +

(5n + 6) (5n - 2) (5n + 7) (5n - 1)

2

22 (+1) |x* - x| 2э (+1) |x* - x|3/2 + 24 (+1) |x* - x|

(5n + 8) 5n (5n + 9) (5n + 1) (5n + 10) (5n + 2) / '

при этом p2 = min ^ p°,-1 = M = max {M°, a}.

4 ^ (M + 1)

M По определению

AyN (x) = |У (x) - yw (x)| =

N

Cn (x* - x)(n-1)/2 - £ Cn (x* - x)(n-1)/2

X

< £ |Cn||x* - x|

n=° n=°

oo

Сп (х* — х)(п-1)/2

п=М +1

Учитывая закономерность образования коэффициентов Сп из теоремы 1, имеем

x*- x|(n-1)/2

n=N+1

£ |Cn||x* - x|(n-1)/2 =

=N+1

X X

= £ |C5k ||x* - x|(5k-1)/2 + £ |C5k+1||x* - x|5k/2 +

fe=N+1 fe=N+1

XX

+ £ |C5k+2 ||x* - x|(5k+1)/2 + £ |C5k+3||x* - x|(5k+2)/2 +

fe=N +1 fc=N+1

oo

+ £ |C5k+4||x* - x|(5k+3)/2 ,

| x - x|

fe=N+1

при N +1 = 5k находим

-1 ^ 25kM (M + 1)k |x* - x|(5k-1)/2 k=N+1 1 - 25 (M +1) |x* - x|

A= E | C5k ||x* - x| 2 ^ oR„, , ,5/2

1 2 |x* - x|1/2 22 |x* - x|

(5k + 2) (5k - 6) + (5k + 3) (5k - 5) + (5k + 4) (5k - 4) +

23 |x* _ x |э/2 24 |x* _ x|

¿-t \<Jb <JU \ ¿J \Jb <JU \

+ ' , '--T +

(5k + 5) (5k - 3) (5k + 6) (5k - 2)

Аналогично получаем оценки для коэффициентов N +1 = 5k + 1;

N + 1 = 5k + 2; N + 1 = 5k + 3; N + 1 = 5k + 4. ►

Пример. Найдем приближенное решение задачи (2), (3) в случае г (х) =0 при начальных данных у (0, 5) = 1, у' (0, 5) = 1/\/3

I 73

и а = 0,001. Эта задача имеет точное решение у = \ --¡=-.

V 1 + 73 — 2х

Найдем радиус аналитичности р2 ^ 0,18946457, точное значение подвижной особой точки х* = (1 + \[Щ /2. Выберем значение аргумента х = 1, 3. С использованием (13) при N =12 вычислим приближенные значения функции:

х у 12 у Ду 12 Ду Д1У

1,3 3,6216787 3,6216776 0,00012 0,0000011 0,000002

Здесь у12 — приближенное решение (13); у — значение точного решения; Ау12 — оценка погрешности приближенного решения (14) по теореме 2; Ау — абсолютная погрешность приближенного решения у12; А1у — апостериорная оценка погрешности, которая определяется путем решения обратной задачи теории погрешности. При е = 0,000002 на основании апостериорной оценки убеждаемся в том, что в структуре приближенного решения (13) должно быть значение N = 18. В таком случае нет необходимости рассчитывать сумму 18 слагаемых, так как для номеров п = 1—5, 7—11, 13—17 коэффициенты Сп = 0. Добавка в структуре приближенного решения для N = 18 не превышает требуемой точности. Поэтому в структуре приближенного решения можем ограничиться значением N =12, при котором приближенное решение, содержащее три отличных от нуля слагаемых, будет иметь погрешность е = 0,000002.

Обсуждение полученных результатов и сопоставление их с ранее известными. Полученные результаты позволяют строить приближенное решение в окрестности подвижной особой точки с любой наперед заданной точностью и в настоящее время не имеют аналогов для рассматриваемого уравнения. Для минимизации структуры приближенного аналитического решения используется апостериорная оценка погрешности.

Заключение. Доказана теорема существования и единственности решения рассматриваемого нелинейного дифференциального уравнения второго порядка, технология доказательства которого позволяет построить приближенное решение в окрестности подвижной особой точки. Получено представление решения нелинейного дифференциального уравнения в ряд по смешанным степеням.

ЛИТЕРАТУРА

1. Орлов В.Н.Метод приближенного решения первого, второго дифференциальных уравнений Пенлеве и Абеля. М.: МПГУ, 2013. 174 с.

2. Коллатц Л.Численные методы решения дифференциальных уравнений. М.: Изд-во иностр. литер., 1953. 460 с.

3. Озерецковский В.Б. Ряды Тейлора как метод решения дифференциальных, интегральных и интегро-дифференциальных уравнений математической физики. М., 1994. 71 с.

4. Чичурин А.В. Использование системы Mathematica при поиске конструктивных методов интегрирования уравнения Абеля // Вучоныя зашсю БрДУ iмя А.С. Пушкша. Брест, 2007. Т. 3. Ч. 2. С. 24-38.

5. Орлов В.Н.Метод приближенного решения скалярного и матричного дифференциальных уравнений Риккати. Чебоксары: Перфектум, 2012. 112 с.

6. Орлов В.Н., Леонтьева Т.Ю. Построение приближенного решения нелинейного дифференциального уравнения в области аналитичности // Междунар. научно-практич. конф. "Фундаментальные и прикладные проблемы механики деформируемого твердого тела, математического моделирования и информационных технологий". 12-15 августа 2013 г. Чебоксары. С. 47-52.

7. Орлов В.Н. Леонтьева Т.Ю. Влияние возмущения начальных данных на приближенное решение одного нелинейного дифференциального уравнения второго порядка в области аналитичности / В.Н. Орлов, Т.Ю. Леонтьева // Вестник Чуваш. гос. пед. ун-та им. И.Я. Яковлева. Сер. Механика предельного состояния. 2013. № 3 (17). С. 103-109.

8. Орлов В.Н., Леонтьева Т.Ю. Построение приближенного решения одного нелинейного дифференциального уравнения второго порядка в области голоморфности // Вестник Чуваш. гос. пед. ун-та им. И.Я.Яковлева. Сер. Естественные и технические науки. 2013. № 4 (80). С. 156-162.

9. Орлов В.Н., Леонтьева Т.Ю. Исследование влияния возмущения начальных данных на приближенное решение одного нелинейного дифференциального уравнения второго порядка в области голоморфности // Вестник РГСУ 2013. № 1 (28). С. 108-111.

10. Голубев В.В. Лекции по аналитической теории дифференциальных уравнений. М.; Л.: Гостехиздат, 1950. 436 с.

11. Орлов В.Н., Гузь М.П. Исследование влияния возмущения подвижной особой точки на приближенное решение задачи Коши одного нелинейного дифференциального уравнения // Междунар. научно-практич. конф. "Фундаментальные и прикладные проблемы механики деформируемого твердого тела, математического моделирования и информационных технологий". 12-15 августа 2013 г. Чебоксары. С. 36-46.

12. Камке Э. Справочник по обыкновенным дифференциальным уравнениям. М.: Наука, 1971. 576 с.

REFERENCES

[1] Orlov V.N. Metod priblizhennogo resheniya pervogo, vtorogo differentsial'nykh uravneniy Penleve i Abelya [Method of approximate solution of the first, second differential Painleve's and Abe'sl equations]. Moscow, MPGU Publ., 2013. 174 p.

[2] Kollatts L. Chislennye metody resheniya differentsial'nykh uravneniy [Numerical methods of differential equations solution]. Moscow, Inostr. Liter. Publ., 1953. 460 p.

[3] Ozeretskovskiy V.B. Ryady Teylora kak metod resheniya differentsial'nykh, integral'nykh i integro-differentsial'nykh uravneniy matematicheskoy fiziki [Taylor's series as a method to solve differential, integral and integro-differential equations]. Moscow, 1994. 71 p.

[4] Chichurin A.V. Using of Mathematica system to find a constructive method for Abel's equation integration. Vuchonyya zapiski BrDU imya A.S. Pushkina [Scientific Articles, Brest State University named after A.S. Pushkin]. Brest, 2007, vol. 3, part. 2, pp. 24-38 (in Russ.).

[5] Orlov V.N. Metod priblizhennogo resheniya skalyarnogo i matrichnogo differentsial'nykh uravneniy Rikkati [Approximate method to solve scalar and matrix Riccati's differential equations]. Cheboksary, Perfektum Publ., 2012. 112 p.

[6] Orlov V.N., Leont'eva T.Yu. Construction of approximate solution of nonlinear differential equation in analyticity region. Mezhdunar. Nauch.-Praktich. Konf. "Fundamental'nye i prikladnye problemy mekhaniki deformiruemogo tverdogo tela, matematicheskogo modelirovaniya i informatsionnykh tekhnologiy" [International Scientific-Training Conf. "Fundamental and Applied Problems of Deformable Solids Mechanics, Mathematical Simulation, and Information Technologies]. Cheboksary. August 12-15, 2013, pp. 47-52 (in Russ.).

[7] Orlov V.N. Influence of initial data perturbation on approximate solution of one nonlinear second-order equation in analyticity region. Vestn. Chuvash. Gos. Ped. Univ. im. I.Ya. Yakovleva. Ser. Mekh. predel'n. sost. [Herald of the I. Yakovlev Chuvash State Pedagogical University, Limit-State Mech.], 2013, no. 3 (17). pp. 103109 (in Russ.).

[8] Orlov V.N., Leont'eva T.Yu. Construction of approximate solution of nonlinear differential equation in holomorphy region. Vestn. Chuvash. Gos. Ped. Univ. im. I.Ya. Yakovleva. Ser. Estestv. i tekhn. nauki [Herald of the I. Yakovlev Chuvash State Pedagogical University, Nat. and Eng. Sci.], 2013, no. 4 (80), pp. 156-162 (in Russ.).

[9] Orlov V.N., Leont'eva T.Yu. Analysis of the influence of initial data perturbation on approximate solution of one nonlinear second-order equation in the holomorphy region]. Vestnik RGSU [Bull. of the Russian State Social University], 2013, no. 1 (28), pp. 108-111 (in Russ.).

[10] Golubev V.V. Lektsii po analiticheskoy teorii differentsial'nykh uravneniy [Lectures on analytical theory of differential equations]. Moscow, Leningrad, Gostekhizdat Publ., 1950. 436 p.

[11] Orlov V.N., Guz' M.P. Analysis of the movable singularity perturbation influence on Cauchy problem approximate solution of one nonlinear differential equation. Mezhdunar. Nauch.-Praktich. Konf. "Fundamental'nye i prikladnye problemy mekhaniki deformiruemogo tverdogo tela, matematicheskogo modelirovaniya i informatsionnykh tekhnologiy" [Inter. Scientific-Training Conf. "Fundamental and Applied Problems of Deformable Solids Mechanics, Mathematical Simulation, and Information Technologies"]. Cheboksary. August 12-15, 2013, pp. 36-46 (in Russ.).

[12] Kamke E. Spravochnik po obyknovennym differentsial'nym uravneniyam [Reference book on ordinary differential equations]. Moscow, Nauka Publ., 1971. 576 p.

Статья поступила в редакцию 22.05.2014

Орлов Виктор Николаевич — д-р физ.-мат. наук, доцент, заведующий кафедрой алгебры и геометрии Чувашского государственного педагогического университета им. И.Я. Яковлева. Автор более 140 научных работ, в том числе монографии и патентов на изобретение, в области аналитической теории дифференциальных уравнений, вычислительной математики, математического моделирования. Чувашский государственный педагогический университет им. И.Я. Яковлева, Российская Федерация, 428000, Чебоксары, ул. К. Маркса, д. 38.

Orlov V.N. — Dr. Sci. (Phys.-Math.), head of the Algebra and Geometry department of the I. Yakovlev Chuvash State Pedagogical University. Author of more than 140 publications, including monographs and patents in the field of analytical theory of differential equations, computational mathematics, mathematical simulation.

I. Yakovlev Chuvash State Pedagogical University, ul. K. Marksa 38, Cheboksary, 428000 Russian Federation.

Леонтьева Татьяна Юрьевна — аспирант кафедры алгебры и геометрии Чувашского государственного педагогического университета им. И.Я. Яковлева. Автор шести научных работ в области аналитической теории дифференциальных уравнений, вычислительной математики и математического моделирования. Чувашский государственный педагогический университет им. И.Я. Яковлева, Российская Федерация, 428000, Чебоксары, ул. К. Маркса, д. 38.

Leont'eva T.Yu. — post-graduate of Algebra and Geometry department of the I. Yakovlev Chuvash State Pedagogical University. Author of six publications in the field of analytical theory of differential equations, computational mathematics, mathematical simulation. I. Yakovlev Chuvash State Pedagogical University, ul. K. Marksa 38, Cheboksary, 428000 Russian Federation.