“…Therefore, the following set of differential equations for f j are obtained as: 11) where j = 0, 1, 2 and f −1 = f −2 = 0. By using the method of eigenfunction expansion, we have solved the coupled system of partial differential equations (3.9) and (3.11) to obtain the unknown functions f 0 (t, r), K 0 (t), f 1 (t, r), K 1 (t), f 2 (t, r) and K 2 (t) in this specific order [7,9]. By following the analysis of Sankarasubramanian & Gill [6] and Rana & Murthy [7][8][9], the function f 0 (t, r) and exchange coefficient K 0 (t) are obtained as follows:…”