Boundary conditions ( i = 0 , n , j = 1 , 2 , ⋯ , m )
u ( 0 , j ) = 0 , u ( n , j ) = 0
Initial Conditions ( i = 1 , 2 , ⋯ , n − 1 , j = 0 )
− 2 ω r u i − 1,1 + 2 ( 1 + 2 ω r ) u i ,1 − 2 ω r u i + 1,1 = ( 1 − 2 ω ) r u i − 1,0 + ( 2 + 4 ω r − 2 r ) u i ,0 + ( 1 − 2 ω ) r u i + 1,0 − 2 ω k r g i − 1 + 2 ( 1 − 2 ω r ) k g i − 2 ω k r g i + 1 (CD) (32)
or
− ω r u i − 1,1 + ( 1 − 2 ω r ) u i ,1 − ω r u i + 1,1 = ( 1 − ω ) r u i − 1,0 + ( 1 + 2 ω r − 2 r ) u i ,0 + ( 1 − ω ) r u i + 1,0 + ω k r ( − g i − 1 + 2 g i − g i + 1 ) + 1 2 ω c 2 k 2 r ( f ″ i − 1 − 2 f ″ i + f ″ i ) + k g i − 1 2 c 2 k 2 f ″ i (BD) (33)
u i , 1 = f i + k g i + 1 2 k 2 c 2 f ″ i (FD) (34)
Rest Values ( i = 1 , 2 , ⋯ , n − 1 , j = 1 , 2 , ⋯ , m )
− ω r u i − 1, j + 1 + ( 1 + 2 ω r ) u i , j + 1 − ω r u i + 1, j + 1 = ( 1 − 2 ω ) r u i − 1, j + ( 2 + 4 ω r − 2 r ) u i , j + ( 1 − 2 ω ) r u i + 1, j + ω r u i − 1, j − 1 + ( − 1 − 2 ω r ) u i , j − 1 + ω r u i + 1, j − 1 (35)