“…Hence, for the system to have a non-zero solution, a i,i must be equal to zero. Therefore, we can regard a i,i as a function of (k, Re, c) and apply an iterative method (we used a modified Newton-Raffson method with regularization multiplier described, for example, in Kogan et al (2000) to find roots of a i,i (c) Z0 for given (k, Re). The code built based on the proposed scheme was tested against results published by Orzag (1971), Blennerhassett (1980) and Yiantsios & Higgins (1988) and demonstrated good accuracy and convergence characteristics.…”