Document Type


Subject Area(s)

Chemical Engineering


Newman's BAND(J) subroutine, which has been used widely to solve models of various electrochemical systems, is extended to solve a system of coupled, ordinary differential equations with interior boundary conditions. A set of coupled, linear ordinary differential equations is used to demonstrate the solution procedure. The results show that the extended technique has the same accuracy as that of using pentadiagonal BAND(J), but the execution speed is about five times faster than that of pentadiagonal BAND(J). Using sparse matrix solver Y12MAF to solve the same set of equations takes even longer time than pentadiagonal BAND(J).