a
i,j
+1
=
K
q
(
l
i
−
1
,j
)
2
+ (
m
i
−
1
,j
)
2
+ (
n
i
−
1
,j
)
2
,
K
=
(
cos(
α
i
−
1
,j
−
α
i,j
)
для случая преломления
,
cos(2
α
i
−
1
,j
)
для случая отражения
,
(14)
b
i,j
+1
= cos (
α
i,j
)
vuut Ã
∂
ˆ
F
i
∂x
!
2
+
Ã
∂
ˆ
F
i
∂y
!
2
+
Ã
∂
ˆ
F
i
∂z
!
2
.
(
15
)
Уравнение
(13)
примет вид
l
i,j
+1
Ã
a
i,j
+1
∂
ˆ
F
i
∂x
−
b
i,j
+1
l
i
−
1
,j
!
+
m
i,j
+1
Ã
a
i,j
+1
∂
ˆ
F
i
∂y
−
b
i,j
+1
m
i
−
1
,j
!
+
+
n
i,j
+1
Ã
a
i,j
+1
∂
ˆ
F
i
∂z
−
b
i,j
+1
n
i
−
1
,j
!
= 0
.
(16)
Из условия
(10),
проведя разложение определителя по второй стро
-
ке
,
получим
l
i,j
+1
Ã
m
i
−
1
,j
∂
ˆ
F
i
∂z
−
n
i
−
1
,j
∂
ˆ
F
i
∂y
!
+
m
i,j
+1
Ã
l
i
−
1
,j
∂
ˆ
F
i
∂z
−
n
i
−
1
,j
∂
ˆ
F
i
∂x
!
+
+
n
i,j
+1
Ã
l
i
−
1
,j
∂
ˆ
F
i
∂y
−
m
i
−
1
,j
+1
∂
ˆ
F
i
∂x
!
= 0
.
(17)
Уравнения
(16), (17)
образуют систему
,
содержащую три неизвест
-
ных
:
l
i,j
+1
,
m
i,j
+1
,
n
i,j
+1
.
Поскольку для определения направления вектора
~a
i,j
+1
(
l
i,j
+1
, m
i,j
+1
, n
i,j
+1
)
величина модуля не имеет значения
,
можем
одну из компонент вектора выбрать произвольно
.
Пусть выбрано не
-
нулевое значение компоненты
n
i,j
+1
=
n
0
i,j
+1
.
Тогда
,
решая систему
уравнений
(16), (17)
относительно
l
0
i,j
+1
, m
0
i,j
+1
,
получим
l
0
i,j
+1
= ˆ
n
0
i,j
+1
à µ
m
i
−
1
,j
∂
ˆ
F
i
∂x
−
l
i
−
1
,j
∂
ˆ
F
i
∂y
¶µ
a
i,j
+1
∂
ˆ
F
i
∂y
−
b
i,j
+1
m
i
−
1
,j
¶
−
−
µ
n
i
−
1
,j
∂
ˆ
F
i
∂x
−
l
i
−
1
,j
∂
ˆ
F
i
∂z
¶µ
b
i,j
+1
n
i
−
1
,j
−
a
i,j
+1
∂
ˆ
F
i
∂z
¶ !
,
(18)
ISSN 1812-3368.
Вестник МГТУ им
.
Н
.
Э
.
Баумана
.
Сер
. “
Естественные науки
”. 2004.
№
4
29