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).
Journal of the Electrochemical Society, 1991, pages 1688-1691.
© The Electrochemical Society, Inc. 1991. All rights reserved. Except as provided under U.S. copyright law, this work may not be reproduced, resold, distributed, or modified without the express permission of The Electrochemical Society (ECS). The archival version of this work was published in the Journal of the Electrochemical Society.