Численные методы моделирования свойств нанокристаллов. Нагорнов Ю.С. - 15 стр.

UptoLike

Составители: 

15
÷
÷
ø
ö
ç
ç
è
æ
-
-=+m
÷
÷
ø
ö
ç
ç
è
æ
-
-
=l+
-
-
N
1NN
N
N
N1NN
1
2
12
2
211
h
yy
'y
h
6
K2K
'y
h
yy
h
6
KK2
(2.11)
где
'
y
1
=
'
f
(x
1
) и
'
y
N
=
'
f
(x
N
),
1
N1
=
m
=
l
.
Коэффициенты сплайна, заданного таким образом, можно вычислить
путем решения системы уравнений
ú
ú
ú
ú
ú
ú
ú
ú
ú
û
ù
ê
ê
ê
ê
ê
ê
ê
ê
ê
ë
é
=
ú
ú
ú
ú
ú
ú
ú
ú
ú
û
ù
ê
ê
ê
ê
ê
ê
ê
ê
ê
ë
é
ú
ú
ú
ú
ú
ú
ú
ú
ú
û
ù
ê
ê
ê
ê
ê
ê
ê
ê
ê
ë
é
m
lm
l
m
lm
l
-
-
-
-
--
-
N
1N
2N
3
2
1
N
1N
2N
3
2
1
N
1N1N
2N
3
22
1
d
d
d
...
d
d
d
K
K
K
...
K
K
K
20...000
2...000
02...000
.....................
000...20
000...2
000...02
(2.12)
для которой
;
h
yy
'y
h
6
d;'y
h
yy
h
6
d
N
1NN
N
N
N1
2
12
2
1
÷
÷
ø
ö
ç
ç
è
æ
-
-=
÷
÷
ø
ö
ç
ç
è
æ
-
-
=
-
1N...2j,
h
yy
h
yy
hh
6
d
j
1jj
1j
j1j
j1j
j
-=
÷
÷
ø
ö
ç
ç
è
æ
-
-
-
+
=
-
+
+
+
(2.13)
Систему уравнений (2.13) можно решить методом исключений Гаусса
(методом прогонки). С помощью прямого хода метода Гаусса получим
систему вида:
ú
ú
ú
ú
ú
ú
ú
ú
ú
û
ù
ê
ê
ê
ê
ê
ê
ê
ê
ê
ë
é
=
ú
ú
ú
ú
ú
ú
ú
ú
ú
û
ù
ê
ê
ê
ê
ê
ê
ê
ê
ê
ë
é
ú
ú
ú
ú
ú
ú
ú
ú
ú
û
ù
ê
ê
ê
ê
ê
ê
ê
ê
ê
ë
é
l
l
l
l
-
-
-
-
--
--
N
1N
2N
3
2
1
N
1N
2N
3
2
1
N
1N1N
2N2N
3
22
11
z
z
z
...
z
z
z
K
K
K
...
K
K
K
p00...000
p0...000
0p...000
.....................
000...p00
000...p0
000...0p
(2.14)
где диагональные элементы p
j
и правые части z
j
по рекуррентным
формулам
1.
.N...2j,
p
z
dz;dz
;N...2j,
p
2p;2p
j
1j
1j
jj11
j
1j
1j
j1
=m-==
=m
l
-==
-
-
-
-
(2.15)