Решение линейных задач математической физики на основе методов взвешенных невязок. Анкилов А.В - 93 стр.

UptoLike

93
Найдем коэффициенты системы дифференциальных уравнений
A
t
H
d
d
CH
B
для отыскания функций H
k
(t) с начальными условиями
AH0() D1
i 1
n

B
i 1
a
b
x
Kx() V2 0 x()
x
Kx() V1 0 x()
d
d
x() V 0 x() gx()
Wix()
d
i1
n

j1
n

A
i 1 j 1
a
b
x
Vjx()Wix()
d
C
i 1 j 1
a
b
x
Kx() V2 j x()
x
Kx()
d
d
V1 j x()
x() Vjx()
Wix()
d
i1
n

D1
i 1
a
b
x
fx() V 0 x()()Wi x()
d
Приведем систему к виду
t
H
d
d
A1 H
B
1
с начальными условиями
H 0() D
2
A1 A
1
C

B1 A
1
B

D2 A
1
D1
Найдем решение получившейся системы дифференциальных уравнений:
HD
2

DtH()A1 H
B
1
Y rkfixed H 0 T 100 D()
Следовательно, при t=T получим следующие коэффициенты
Y
100 k
2.8878
0.0481
2.0958 10
3
6.9289 10
5
1.3289 10
6