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

UptoLike

Рубрика: 

151
M1
0.000000
0.000000
0.000000
0.000000
0.000000
0.000000
0.000000
0.000000
0.000000
0.000000
0.000000
0.000000
0.000000
0.000000
0.000000
0.000000
0.000000
0.000000
0.000000
0.000000
0.000000
0.000000
0.000000
0.000000
0.000000
=
C1 A
1
C:=
C1
5.471344
21.446746
43.838749
42.803134
16.255657
5.117974
10.942688
2.158284
21.559281
15.205775
4.404850
45.206238
130.400363
157.454762
70.670055
11.523921
88.877037
243.679989
294.236717
132.061581
18.964933
141.920564
395.415783
484.225799
217.333939
=
B1 A
1
B:=
B1
T
00000()=
D2 A
1
D1:=
D2
T
3.194 9.53 10
12
× 2.363 10
11
× 2.021 10
11
× 3.535 10
12
×
()
=
N2 A
1
N1:=
N2
T
00000()=
Приведем к нормальной системе дифференциальных уравнений первого
порядка
t
H
d
d
AA H BB+
с начальными условиями
H 0() D2
i1n..:=
D2
ni+ 1
N2
i 1
:=
D2
T
3.194 9.53 10
12
× 2.363 10
11
× 2.021 10
11
× 3.535 10
12
× 00000
( )
=
i1n..:=
j1n..:=
AA
i 1 j 1,
0:=
AA
ni+ 1 nj+ 1,
M1
i 1 j 1,
:=
AA
ni+ 1 j 1,
C1
i 1 j 1,
:=
AA
i 1 nj+ 1,
if ij1, 0,():=