Ngày tải lên :
30/09/2013, 03:20
...
223
INTEGER IP1, IM1, JP1, JM1
N1=N -1
M1=M -1
DO I=2,N1
DO J=2,M1
IP1=I +1
IM1=I -1
JP1=J +1
JM1=J -1
LAPU(I,J)=(U(IP1,JP1)+U(IP1,JM1)+U(IM1,JP1) &
+U(IM1,JM1)+4.*(U(IP1,J)+U(IM1,J)+U(I,JP1) ... A(i1,i)
A(i1 ,1: N +1) =A(i1 ,1: N +1) -A(i ,1: N +1) *Tmp
endif
enddo
211
()
⎟
⎟
⎟
⎟
⎟
⎠
⎞
⎜
⎜
⎜
⎜
⎜
⎝
⎛
=
mmmm
m
m
jk
RRR
RRR
RRR
R
21
222 21
112 11
(9.3 .13 ’)
Một cách tổng quát, hệ số ...
j=
0 ,1, 2
, ,NT
.
Từ (9.4.25) ta có hệ phương trình:
0, 1 1, 1 2, 1 0, 1, 2,
1, 1 2, 1 3, 1 1, 2, 3,
2, 1 3, 1 4, 1 2, 3, 4,
2 ,1 1 ,1 ,1 2,
(2 2) (2 2)
(2 2 ) (2 2 )
(2 2 ) (2 2 )
(2 2 )
+++
+++
+++
−+...