MATLAB code for John Newman's BAND numerical method.
Modified from original FORTRAN code in Appendix C of "Electrochemical Systems" 3rd ed. by John Newman and Karen Thomas-Alyea.
Original publication for BAND method: Newman, John. Numerical Solution of Coupled, Ordinary Differential Equations. I&EC Fundamentals. 7:3, 1968. pp. 514-517
Test problem is a ternary diffusion problem using Stefan-Maxwell equations and two Dirichlet boundary conditions, from Chapter 2 in "Multicomponent Mass Transfer" by Ross Taylor & R. Krishna (example 2.1.1).
This program includes 5 files:
autoband_test: run this file. defines operating conditions, number of variables, initial guess, intializes variables
autoband: calculates derivatives of each governing equation by each variable, at mesh points j-2, j-1,j,j+1,j+2
band: solves for and returns change variables
eqn: includes governing equations (must sum to zero)
matinv: inverts banded matrix in band algorithm