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

UptoLike

115
График точного решения при t=T
02
0
1
2
UT x T()
x
Скопируйте график полученной интегральной кривой в файл отчета.
Получение приближенного решения
Введите порядок пробного решения
U
n
=V(0,x)+
n
k
tkHxkV
1
),(),( .
n 5
1. Введите систему пробных функций:
k 1
n

V0 k x() xa()
k
xb()
Нормируем их. Для этого вычислим нормировочные коэффициенты
i 1
n

VV
i 1
a
b
x
V0 i x()()
2
d

Получили нормированные пробные функции
k 0
n

Vkx()if k 0
V0 k x()
VV
k 1
b0 a2
b
b2 a0
a
a0 b0 ba()
b2 a0
b0 a2
()x
a0 b0 ba()

Введите функции V1(k,x) и V2(k,x), равные первой и второй
производной от функции V(k,x)
k 0
n

V1 k x()if k 0
xa()
k
xb()k xa()
k 1
VV
k 1
b2 a0 b0 a2
a0 b0 ba()

V2 k x()if k 0 if k 1
2k xa()
k 1
xb()k k 1() xa()
k 2
VV
k 1
2
VV
k 1
0

Введите систему поверочных функций:
Wkx
()Vkx
()