Ngày tải lên :
30/09/2013, 03:20
...
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) &
+U(I,JM1))-20.*U(I,J))/(6.*H**2) ...
NX = L/h
. Do đó
i=
0 ,1, 2
, ,NX
,
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 ... I
A(i ,1: N +1) = A(i ,1: N +1) / Tmp
do i1 =1, N !
Khử các phần tử ngoài đường chéo
if (i1 /= i) then
Tmp = 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
...