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

UptoLike

118
02
0
0.01
0.02
UT x T()Ux()
x
Замените старое значение меры точности
11
наибольшим значением
UT x() Ux()
на отрезке [a,b] (для этого необходимо кликнуть мышью по
графику, где в левом верхнем углу появится наибольшее значение,
скопировать это значение и заменить на него уже имеющееся ниже
значение
11
)
11 0.019
Получим матрицу предыдущего (для 4
n) пробного решения
AP submatrix A 0 n 2 0 n 2
()
MP submatrix M 0 n 2 0
n 2
()
CP submatrix C 0 n 2 0 n 2
()
BP submatrix B 0 n 2 0 0
()
D1P submatrix D1 0 n 2 0
0
()
N1P submatrix N1 0 n 2 0
0
()
M1P AP
1
MP

C1P AP
1
C
P

B1P AP
1
BP

D2P AP
1
D1
P

N2P AP
1
N1

i1n1
j1n1

AAP
i 1 j 1
0

AAP
ni 2 nj 2
M1P
i 1 j 1

AAP
ni 2 j 1
C1P
i 1 j 1

AAP
i 1 nj 2
if ij1 0()
i1n1
D2P
ni 2
N2P
i 1

BBP
i 1
0

BBP
ni 2
B1P
i 1

HP D2
P

DtHP()AAP HP
B
B
P

YP rkfixed HP 0 T 100 D()
Следовательно, предыдущее пробное решение U(x) для
n 5
при t=T имеет
вид