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)

or

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)