“…Then, starting from initial guesses, we obtain approximate solutions for uðx, z, tÞ, wðx, z, tÞ, Tðx, z, tÞ, θðx, z, tÞ, and Cðx, z, tÞ by iteratively solving the matrix equations ( 73), ( 81), ( 91), (99), and (107), in turn, for s = 0, 1, 2, ⋯. We only confirm convergence to the expected solutions by measuring the norm residual errors of the difference between successive iterations as defined in (64).…”