Chương 3: Nội suy | Bài giảng môn Phương pháp tính và matlab CTTT | Đại học Bách khoa hà nội

Chương 3: Nội suy. Tài liệu trắc nghiệm môn Phương pháp tính và matlab CTTT giúp bạn tham khảo, ôn tập và đạt kết quả cao. Mời bạn đọc đón xem!

Trường:

Đại học Bách Khoa Hà Nội 2.8 K tài liệu

Thông tin:
27 trang 3 tháng trước

Bình luận

Vui lòng đăng nhập hoặc đăng ký để gửi bình luận.

Chương 3: Nội suy | Bài giảng môn Phương pháp tính và matlab CTTT | Đại học Bách khoa hà nội

Chương 3: Nội suy. Tài liệu trắc nghiệm môn Phương pháp tính và matlab CTTT giúp bạn tham khảo, ôn tập và đạt kết quả cao. Mời bạn đọc đón xem!

57 29 lượt tải Tải xuống
C
C
C
h
h
h
ư
ư
ư
ơ
ơ
ơ
n
n
n
g
g
g
I
I
I
I
I
I
I
I
I
:
:
:
N
N
N
O
O
O
Ä
Ä
Ä
I
I
I
S
S
S
U
U
U
Y
Y
Y
1
1
)
)
)
N
N
N
o
o
o
ä
ä
ä
i
i
i
s
s
s
u
u
u
y
y
y
ñ
ñ
ñ
a
a
a
t
t
t
h
h
h
ö
ö
ö
ù
ù
ù
c
c
c
2
2
)
)
)
N
N
N
o
o
o
ä
ä
ä
i
i
i
s
s
s
u
u
u
y
y
y
S
S
S
p
p
p
l
l
l
i
i
i
n
n
n
e
e
e
b
b
b
a
a
a
ä
ä
ä
c
c
c
3
3
3
3
3
3
)
)
)
P
P
P
h
h
h
ö
ö
ö
ô
ô
ô
n
n
n
g
g
g
p
p
p
h
h
h
a
a
a
ù
ù
ù
p
p
p
b
b
b
ì
ì
ì
n
n
n
h
h
h
p
p
p
h
h
h
ö
ö
ö
ô
ô
ô
n
n
n
g
g
g
t
t
t
o
o
o
á
á
á
i
i
i
t
t
t
h
h
h
i
i
i
e
e
e
å
å
å
u
u
u
Ngô Thu Lương – Phương Pháp Tính
1
1
1
1
.
.
.
1
1
1
)
)
)
N
N
N
o
o
o
ä
ä
ä
i
i
i
s
s
s
u
u
u
y
y
y
ñ
ñ
ñ
a
a
a
t
t
t
h
h
h
ö
ö
ö
ù
ù
ù
c
c
c
t
t
t
h
h
h
e
e
e
o
o
o
L
L
L
a
a
a
g
g
g
r
r
r
a
a
a
n
n
n
g
g
g
e
e
e
a
a
a
)
)
)
N
N
N
o
o
o
ä
ä
ä
i
i
i
d
d
d
u
u
u
n
n
n
g
g
g
:
Bieát caùc giaù trò
( )
i i
y f x
=
cuûa haøm
( )
y f x
=
taïi caùc ñieåm
i
x
theo baûng
Ngô Thu Lương – Phương Pháp Tính
2
m haøm laïi haøm
( )
f x
Lôøi giaûi : V
V
V
o
o
o
â
â
â
s
s
s
o
o
o
á
á
á
h
h
h
a
a
a
ø
ø
ø
m
m
m
Tìm
( ) ( )
f x P x
=
c laø ñ
ñ
ñ
a
a
a
t
t
t
h
h
h
ö
ö
ö
ù
ù
ù
c
c
c
b
b
b
a
a
a
ä
ä
ä
c
c
c
n
n
n
thoûa
i
i
y
x
P
=
)
(
øi giaûi ld
d
d
u
u
u
y
y
y
n
n
n
h
h
h
a
a
a
á
á
á
t
t
t
Caùc böôùc tìm ña thöùc
)
(
x
P
B
B
ö
ö
ô
ô
ù
ù
c
c
1
1
:
: Thieát laäp ñ
ñ
ñ
a
a
a
t
t
t
h
h
h
ö
ö
ö
ù
ù
ù
c
c
c
c
c
c
ô
ô
ô
s
s
s
ô
ô
ô
û
û
û
L
L
L
a
a
a
g
g
g
r
r
r
a
a
a
n
n
n
g
g
g
e
e
e
=
=
n
ikk
ki
k
i
xx
xx
xL
,0
)(
)(
)(
V
V
í
í
d
d
u
u
ï :
=
)
(
0
x
L
)
)..(
)(
)...(
(
x
x
x
x
x
x
x
x
Ngô Thu Lương – Phương Pháp Tính
3
))..()()...((
)
)..(
)(
)...(
(
0
0
1
0
1
0
11
n
i
i
nii
xxxxxxxx
x
x
x
x
x
x
x
x
=
B
B
ö
ö
ô
ô
ù
ù
c
c
2
2
:
:
C
C
C
o
o
o
â
â
â
n
n
n
g
g
g
t
t
t
h
h
h
ö
ö
ö
ù
ù
ù
c
c
c
t
t
t
í
í
í
n
n
n
h
h
h
P
P
P
(
(
(
x
x
x
)
)
)
=
=
n
i
ii
xLyxP
0
)()(
=
)
(
...
)
(
)
(
1100
x
L
y
x
L
y
x
L
y
nn
+
+
+
b
b
b
)
)
)
S
S
S
a
a
a
i
i
i
s
s
s
o
o
o
á
á
á
:
( ) ( )
f x P x
Ngô Thu Lương – Phương Pháp Tính
4
( ) ( )
f x P x
))....()((
)!1(
10
)
1
(
n
n
xxxxxx
n
M
+
+
c
c
c
)
)
)
N
N
N
h
h
h
a
a
a
ä
ä
ä
n
n
n
x
x
x
e
e
e
ù
ù
ù
t
t
t
:
*) Soá moác noäi suy caøng lôùn thì sai soá caøng nhoû , tuy
nhieân baäc cuûa ña thöùc seõ lôùn, tính toaùn seõ daøi .
*)Sai sph thuoäc vaøo
)
1
(
+
n
M
, tïc t khoâng bieát
hm
( )
f x
chöa bieát
*)Ña thöùc noäi suy
)
(
x
P
laø duy nhaát
Ngô Thu Lương – Phương Pháp Tính
5
V
V
V
í
í
í
d
d
d
u
u
u
ï
:
Tìm ña tùc noäi suy P(x) ø baûng soá lieäu
1
,
0
,
1
2
1
0
=
=
=
x
x
x
3,1,
3
1
210
=== yyy
Tính gaàn ñuùng giaù trò cuûa bng taïi
7
.
0
=
x
Ngô Thu Lương – Phương Pháp Tính
6
G
G
G
i
i
i
a
a
a
û
û
û
i
i
i
: Ta tìm caùc ña thöùc Lagrange
2)11)(01(
)1)(0(
)(
2
0
xxxx
xL
=
=
1
1
)10)](1(0[
)1)](1([
)(
2
1
=
=
xxx
xL
)
0
)](
1
(
[
)
(
2
x
x
x
x
x
L
+
=
=
Ngô Thu Lương – Phương Pháp Tính
7
2)01)](1(1[
)
0
)](
1
(
[
)
(
2
x
x
x
x
x
L
+
=
=
3
342
)(3)(1)(
3
1
)(
2
210
++
=++=
xx
xLxLxLxP
2.26
3
3)7.0.(4)7.0.(2
)7.0(
2
=
++
=P
d
d
d
)
)
)
T
T
T
y
y
y
û
û
û
s
s
s
a
a
a
i
i
i
p
p
p
h
h
h
a
a
a
â
â
â
n
n
n
Tyû sai phaân baäc 0 cuûa
f
taïi
0
x
:
)
(
]
[
0
0
x
f
x
f
=
Tyû sai phaân baäc 1 cuûa
f
taïi
1
,
0
x
x
:
0
1
01
10
]
[
]
[
],[
x
x
x
f
x
f
xxf
=
Ngô Thu Lương – Phương Pháp Tính
8
0
1
x
x
Tyû sai phaân baäc 2 cuûa
f
taïi
2
,
,
10
x
x
x
02
1021
210
]
,
[
]
,
[
],,[
xx
x
x
f
x
x
f
xxxf
=
Töông töï cho tyû sai phaân baäc cao hôn
e
e
e
)
)
)
B
B
B
a
a
a
n
n
n
g
g
g
t
t
t
y
y
y
s
s
s
a
a
a
i
i
i
p
p
p
h
h
h
a
a
a
â
â
â
n
n
n
Ngô Thu Lương – Phương Pháp Tính
9
f
f
f
)
)
)
N
N
N
o
o
o
ä
ä
ä
i
i
i
s
s
s
u
u
u
y
y
y
N
N
N
e
e
e
w
w
w
t
t
t
o
o
o
n
n
n
t
t
t
i
i
i
e
e
e
á
á
á
n
n
n
t
t
t
h
h
h
e
e
e
o
o
o
b
b
b
a
a
a
û
û
û
n
n
n
g
g
g
t
t
t
y
y
y
û
û
û
s
s
s
a
a
a
i
i
i
p
p
p
h
h
h
a
a
a
â
â
â
n
n
n
Ña thöùc
)
(
x
P
coù theå tìm döôùi daïng
..
)
)(
(
)
(
)
(
1
0
2
0
1
0
+
+
+
=
x
x
x
x
a
x
x
a
a
x
P
)
)..(
)(
(
..
1
1
0
+
n
n
x
x
x
x
x
x
a
Ngô Thu Lương – Phương Pháp Tính
10
Ngô Thu Lương – Phương Pháp Tính
11
1
3
42
3
2
)0)(1(
3
2
)1(
3
2
3
1
)( ++=++++= xxxxxxP
g
g
)
)
N
N
o
o
ä
ä
i
i
s
s
u
u
y
y
N
N
e
e
w
w
t
t
o
o
n
n
l
l
u
u
ø
ø
i
i
..
)
)(
(
)
(
)
(
1210
+
+
+
=
nnn
x
x
x
x
a
x
x
a
a
x
P
)
)..(
)(
(
...
1
1
x
x
x
x
x
x
a
n
n
n
+
]
[
0
n
x
f
a
=
Ngô Thu Lương – Phương Pháp Tính
12
0
n
=
1
a
]
[
1
,
n
n
x
x
f
.
]
,
,
[
212
=
nnn
x
x
x
f
a
]
,
,
[
1...,1 knknnnk
x
x
x
x
f
a
+
=
1 3 2 1 0
[ , ,.. , , , ]
n n n
a f x x x x x x
=
Ngô Thu Lương – Phương Pháp Tính
13
)0)(1(
3
2
)1(23)( ++= xxxxP
2
2 4
1
3 3
x x
+ +
=
2
2
2
)
)
)
N
N
N
o
o
o
ä
ä
ä
i
i
i
s
s
s
u
u
u
y
y
y
S
S
S
p
p
p
l
l
l
i
i
i
n
n
n
e
e
e
b
b
b
a
a
a
ä
ä
ä
c
c
c
3
3
3
a
a
a
)
)
)
N
N
N
o
o
o
ä
ä
ä
i
i
i
d
d
d
u
u
u
n
n
n
g
g
g
:
:
:
Cho baûng soá lieäu
Tìm
m
m
o
o
ä
ä
t
t
h
h
a
a
ø
ø
m
m
)
(
x
S
thoûa caùc ñieàu kieän :
Ngô Thu Lương – Phương Pháp Tính
14
Tìm
m
m
o
o
ä
ä
t
t
h
h
a
a
ø
ø
m
m
)
(
x
S
thoûa caùc ñieàu kieän :
)
(
x
S
: Đi qua các đim đã cho trong bng
)
(
x
S
laø
ñ
ñ
a
a
t
t
h
h
ö
ö
ù
ù
c
c
b
b
a
a
ä
ä
c
c
3
3
t
t
r
r
e
e
â
â
n
n
m
m
o
o
ã
ã
i
i
ñ
ñ
o
o
a
a
ï
ï
n
n
n
n
h
h
o
o
û
û
]
,
[
1+jj
x
x
( caùc ña thöùc naøy coù caùc heä soá khaùc nhau)
Goïi
( )
j
S x
laø ña thöùc treân moãi ñoaïn nhoû
1
[ , ]
j j
x x
+
( )
j
S x
thoûa
caùc ñieàu kieän
:
a
a
a
)
)
)
jjj
y
x
S
=
)
(
1
1
)
(
++
=
jjj
y
x
S
b
b
b
)
)
)
)()(
1
/
1
1
/
+
+
+
=
j
j
j
j
xSxS
Ngô Thu Lương – Phương Pháp Tính
15
1
1
1
+
+
+
j
j
j
j
c
c
c
)
)
)
)()(
1
//
1
1
//
+
+
+
=
j
j
jj
xSxS
d
d
d
)
)
)
// //
0 0 1
( ) ( )
n n
S x S x
=
ñ
ñ
ñ
i
i
i
e
e
e
à
à
à
u
u
u
k
k
k
i
i
i
e
e
e
ä
ä
ä
n
n
n
b
b
b
i
i
i
e
e
e
â
â
â
n
n
n
t
t
t
ö
ö
ö
ï
ï
ï
n
n
n
h
h
h
i
i
i
e
e
e
â
â
â
n
n
n
jjj
x
x
h
=
+1
jj
y
a
=
j
jj
j
h
c
c
d
3
)
(
1
=
+
3
)
2
(
)
(
11 jjj
j
jj
j
c
c
h
h
a
a
b
+
=
++
Ngô Thu Lương – Phương Pháp Tính
16
3
j
h
Ñeå tìm
j
c
ta giaûi töø heä
b
Ax
=
+
+
+
=
100000
)(2.00
0...00
0.)(20
0..)(2
000001
1122
2211
1100
nnnn
hhhh
hhhh
hhhh
A
c
Ngô Thu Lương – Phương Pháp Tính
17
=
n
n
c
c
c
c
x
1
1
0
.
=
0
)
21
(
2
3
)
1
(
1
3
.
.
)
01
(
0
3
)
12
(
1
3
0
n
a
n
a
n
h
n
a
n
a
n
h
aa
h
aa
h
B
V
V
í
í
d
d
u
u
ï
ï
:
:
Noäi suy Spline baäc 3 cuûa baûng
3
2
1
0
3
2
1
0
=
=
=
=
x
x
x
x
0
4
1
0
3
2
1
0
=
=
=
=
y
y
y
y
1
0
1
1
0
0
=
=
=
=
y
a
y
a
4
=
=
y
a
0
3
3
=
=
y
a
Caùc heä s
i
c
tính theo heä phöông trình
0
0
0
0
1
0
c
0
0
c
Ngô Thu Lương – Phương Pháp Tính
18
=
0
21
6
0
1000
1410
0141
0
0
0
1
3
2
1
0
c
c
c
c
=
0
6
3
0
3
2
1
0
c
c
c
c
0 1 2
0 3 0
b b b
= = =
2
3
1
2
1
0
=
=
=
d
d
d
Ta coù haøm :
=
)
(
x
S
+
++
32)2(2)2(64
21)1(3)1(3)1(31
10)0(1
32
32
3
xxx
xxxx
xx
Ngô Thu Lương – Phương Pháp Tính
19
Spline vi ñieàu kieän bieân raøng buoäc
/ /
0 0 0 1
( ) '( ), ( ) '(
)
)
n n n
S x f x S x fd
x
= =
trong ñoù
)
(
'
,
)
(
'
0
n
x
f
x
f
laø caùc ñaïi löôïng cho tröôùc
0 0
0 0 1 1
2 0 0 0 0
2( ) . . 0
h h
h h h h
+
+
Ngô Thu Lương – Phương Pháp Tính
20
1 1 2 2
2 2 1 1
1 1
0 2( ) . 0
0 0 . . . 0
0 0 . 2( )
0 0 0 0 2
n n n n
n n
h h h h
A
h h h h
h h
+
=
+
| 1/27

Preview text:

Chương III : NOÄI SUY
1) Noäi suy ña thöùc
2) Noäi suy Spline baäc 3
3) Phöông phaùp bình phöông toái thieåu
Ngô Thu Lương – Phương Pháp Tính 1
1.1) Noäi suy ña thöùc theo Lagrange
a) Noäi dung : Bieát caùc giaù trò y = f (x ) cuûa haøm i i y = f ( )
x taïi caùc ñieåm x i theo baûng
Tìm haøm laïi haøm f (x)
Lôøi giaûi : Voâ soá haøm
Tìm f (x) = P(x) chæ laø ña thöùc baäc n
thoûa P( ix) = iy
Lôøi giaûi laø duy nhaát
Ngô Thu Lương – Phương Pháp Tính 2
Caùc böôùc tìm ña thöùc P(x)
Böôùc 1 : Thieát laäp ña thöùc cô sôû Lagrange n ( x − ) L ( x ) = ∏ xk i ( x 0, i k = k i xk )
Ví duï : L ( ) 0 x = (x x ). ) . . . . ( . x ( − − − − 1 x x )( ) x i−1 − x ).)..(.x ( i x ) n =
(x x )...(x x )(x
x )..(x x ) 0 1 0 i 1 0 i 0 n
Ngô Thu Lương – Phương Pháp Tính 3
Böôùc 2 : Coâng thöùc tính P(x) n
P( x) = ∑ yiL i (x) = i =0
y L ( x) + y L ( x) + ... + y L ( x) 0 0 1 1 n n b) Sai soá :
f ( x) − P ( x) ≤ (n + ) 1 M
( x x )( x x )....( x x ) 0 1 n (n + ) 1 !
Ngô Thu Lương – Phương Pháp Tính 4 c) Nhaän xeùt :
*) Soá moác noäi suy caøng lôùn thì sai soá caøng nhoû , tuy
nhieân baäc cuûa ña thöùc seõ lôùn, tính toaùn seõ daøi .
*)Sai soá phuï thuoäc vaøo (n+ ) 1 M , thöïc teá khoâng bieát
vì haøm f (x) chöa bieát
*)Ña thöùc noäi suy P(x) laø duy nhaát
Ngô Thu Lương – Phương Pháp Tính 5 Ví duï :
Tìm ña thöùc noäi suy P(x) töø baûng soá lieäu
x = −1 , x = 0 , x = 1 0 1 2 1 y = , y = 1 , y = 3 0 1 2 3
Tính gaàn ñuùng giaù trò cuûa bảng taïi x = 0.7
Ngô Thu Lương – Phương Pháp Tính 6
Giaûi : Ta tìm caùc ña thöùc Lagrange 2 − − − (x ) 0 (x ) 1 x x L (x) 0 = = ( 1 − − ) 0 ( 1 − − ) 1 2 2 x − − x x − [ ( ) 1 ]( ) 1 1 L (x) 1 = = 0 [ − (− ) 1 ] 0 ( − ) 1 −1 2 [x − (− ) 1 ](x − ) 0 x + [x ( ) 1 ](x ) 0 x x L (x) 2 = = 1 [ − (− ) 1 ] 1 ( − ) 0 2 1 2 2 x + 4 + 3 ( P ) x = x L ( ) x +1L ( ) x + 3L ( ) x = 0 1 2 3 3 2.(0.7)2 + 4.(0.7) + 3 P(0.7) = = 2.26 3
Ngô Thu Lương – Phương Pháp Tính 7 d) Tyû sai phaân
Tyû sai phaân baäc 0 cuûa f taïi 0 x :
f [x ] = f (x ) 0 0
Tyû sai phaân baäc 1 cuûa f taïi x , : 0 1 x
f [x ] − f [x ] 1 0 f [ 0 x , x1] = x1 − 0 x
Tyû sai phaân baäc 2 cuûa f taïi x , x , 0 1 x 2
f [x1, x 2] − f [ 0 x , x1] f [ 0
x , x1, x2] = x 2 − 0 x
Töông töï cho tyû sai phaân baäc cao hôn
Ngô Thu Lương – Phương Pháp Tính 8
e) Baûng tyû sai phaân
Ngô Thu Lương – Phương Pháp Tính 9
f) Noäi suy Newton tieán theo baûng tyû sai phaân
Ña thöùc P(x) coù theå tìm döôùi daïng
P(x) = a + a (x x ) + a (x x )(x x ) + .. 0 1 0 2 0 1 .. + a (
n x x )(x x )..(x x ) 0 1 n 1 −
Ngô Thu Lương – Phương Pháp Tính 10 1 2 2 2 2 4 ( P ) x = + (x + ) 1 + (x + ) 1 (x − ) 0 = x + x +1 3 3 3 3 3
Ngô Thu Lương – Phương Pháp Tính 11
g) Noäi suy Newton luøi ( P )
x = a + a (x x )
n + a (x x )( n x x ) − n +.. 0 1 2 1
... + a (x x )(x x )..( x − − x ) n n n 1 1
a0 = f [ x ] 0 n . 1
a = f [x , x ] n n 1 − a
= f [ x , x , x ] 2 n n 1 − n−2
a = f [x , x x , x ] k n n− , 1 ... nk 1 + nk
a = f [ x , x
1 , ..x3 , x2 , 1 x , x − 0 ] n n n
Ngô Thu Lương – Phương Pháp Tính 12 2
P( x) = 3 + 2( x − ) 1 + ( x − ) 1 ( x − 0) 3 2 2 4 = x + x + 1 3 3
Ngô Thu Lương – Phương Pháp Tính 13
2) Noäi suy Spline baäc 3
a) Noäi dung : Cho baûng soá lieäu Tì T m ì moäot t h a
h øam S(x) th t oûoa ca c ùc c ñie i àeu kie i äen : :
S(x) : Đi qua các điểm đã cho trong bảng
S(x) laø ña thöùc baäc 3 treân moãi ñoaïn nhoû [ x , ] j x j +1
( caùc ña thöùc naøy coù caùc heä soá khaùc nhau)
Ngô Thu Lương – Phương Pháp Tính 14 Goïi S ( )
j x laø ña thöùc treân moãi ñoaïn nhoû [ x , x 1] j j + S ( )
j x thoûa caùc ñieàu kieän :
a) S j (x j ) = y j S ( j x j+1) = y j 1 + b) / S (x ) / j j+1 = S (x ) +1 + j j+ j j +1 1 c) // S (x ) // j j = S (x ) +1 + j j +1 1 d) // / / S = 0 ( 0 x ) S 1(x ) nn
ñieàu kieän bieân töï nhieân
Ngô Thu Lương – Phương Pháp Tính 15 h = +1 − j x j x j a j = y j (c j c 1 − j ) d j = + h 3 j
(a + − a ) h (c + + 2c ) j 1 j j j 1 j b = − j h 3 j
Ngô Thu Lương – Phương Pháp Tính 16
Ñeå tìm c ta giaûi töø heä j Ax = b  1 0 0 0 0 0    h 2( 0 0 h + h ) h . . 0 1 1   0 h 2( 1 h1 + h ) h . 0  A =  2 2   0 0 . . . 0    0 0 . n h − 2( 2 n h −2 + n h − ) 1 n h −  1  0 0 0 0 0 1   c  0  0    c  3 3   1   ( aa ) − ( aa ) 2 1 1 0  h h    1 0  x =   .  .  B =      .  cn −1   3 3  ( aa ) − ( aa )  n n −1 n −1 n − 2 c    h hn n −1 n − 2    0 
Ngô Thu Lương – Phương Pháp Tính 17
Ví duï : Noäi suy Spline baäc 3 cuûa baûng x = 0 x = 1 x = 2 x = 3 0 1 2 3 y = 0 y = 1 y = 4 y = 0 0 1 2 3 a = y = 0 a = y = 1 0 0 1 1 a = y = 4 2 2 a = y = 0 3 3
Caùc heä soá ic tính theo heä phöông trình 1  0 0  0 c   0  0 cc   0  0 0 c           1 4 1  0  1 c   6  =  1 c   3  =  0 1 4 1 c   c   2 − 6 2 − 2  1           0 0 0 1  3 c   0   3 c   0  = = = 0 b 0 1 b 3 2 b 0 d = 1 d = −3 d = 2 0 1 2
Ngô Thu Lương – Phương Pháp Tính 18
Ta coù haøm : S (x) =  ( 1 x − ) 0 3 0≤ x ≤1   1+ ( 3 x − ) 1 + ( 3 x − ) 1 2 − ( 3 x − ) 1 3 1≤ x ≤ 2  4 − ( 6 x − ) 2 2 + ( 2 x − ) 2 3 2≤ x ≤  3
Ngô Thu Lương – Phương Pháp Tính 19
Spline với ñieàu kieän bieân raøng buoäc / / d ) S = = 0 ( 0 x ) f '( 0
x ) , S 1(x ) f '( x ) nn n
trong ñoù f '(x ), f '( ) laø caùc ñaïi löôïng cho tröôùc 0 n x 2   0 h 0 h 0 0 0 0    + 0 h 2( 0 h h1) 1 h . . 0    0 + 1 h 2( 1 h 2 h ) 2 h . 0 A =    0 0 . . . 0    0 0 . +  n h 2 2( − n h 2 − n h 1) − n h 1 −   0 0 0 0   n h 1 2 − n h 1 − 
Ngô Thu Lương – Phương Pháp Tính 20