Geoscience Reference
In-Depth Information
From ( 2.59 ) and ( 2.60 ) we get:
Z
x k þ 1
k h k þ h k X
X
n 1
s¼2 e
n
h s 2
k
s
0
ð s Þ
k
½P i y ð n i Þ P i y ð n i Þ dx
e k þ 1 ¼ e k þ e
!
i¼1
x k
n
Z
Z
x k þ 1
x k þ 1
½f if k dx þ k
if W½y W k ½y g dx
þ
x k
x k
n
n
We know that
Z
x k þ 1
Z
x k þ 1
if k dx ¼ if k h k
n !
W k ½ y dx ¼ W k ½y
n !
;
x k
x k
n
n
Let us denote
;
;
ð j Þ
k
y ð n i Þ
E k ¼ max
j
e
h max ¼ max
k
h k ;
p i ¼ max
½a ; b
P j ;
M n i ¼ max
½a ; b
;
;
;
P ik
if k
y ð n i Þ
L i ¼ max
½a ; b
N n i ¼ max
½a ; b
F ¼ max
½a ; b
j ;
G 0 ¼ max
½a ; b
j X
;
r
M j max
G
W k ½y
T ¼ max
G
s ¼ b a
j
max
G
K j ð x
; nÞ
j
W½y
j;
j¼0
"
#
X
n
1
n
l ¼
ð
p i M n i þ L i N n i
Þ þ F þ G 0 þ jjð T þ s Þ
;
!
i¼1
g ¼ 1 þ X
n 1
h s 1
ax
s
;
!
s¼2
and M n i
N n i . Then we get the following recurrent error estimation
E k þ 1 ð 1 þ gh max Þ E k þ lh n max
ð 2
:
61 Þ
Hence we get
k
k
e 0 þ lh n 1
g 1
E k ð 1 þ gh max Þ
½ ð 1 þ gh max Þ
1
ð 2
:
62 Þ
ax
where
ʵ 0 is the maximum error in the initial data. Obviously, if
ʵ 0 = 0, then from
( 2.62 ) it follows that if h max
0 then E k
0, i.e.
~
y ð x k Þ! y ð x k Þ .
Search WWH ::




Custom Search