A framework model of single-proton conduction through gramicidin was previously designed to incorporate potentials of mean force and diffusion coefficients computed by the molecular dynamics simulations of Pom s and Roux (1). The resulting diffusion model was solved analytically using the lumped state approximation (LSA), allowing a detailed comparison to be made with conductance data from gramicidin A and two Trp→ Phe analogs (2). The comparison included a sensitivity analysis which required over 1 million current evaluations. A numerical method for constructing framework models is now introduced which involves finding the steady states of random walks using a trapezoid rule closely related to the rule for numerical integration. The method is described and then applied directly to the LSA. Convergence of the results to the analytical solution is seen as the number of random walk sites increase. The numerical method is then used to construct a more elaborate framework model which avoids the LSA. This is also in very good agreement with the analytical solution under the experimental conditions, confirming the accuracy of the LSA. The numerical method remains fast enough to allow an extensive comparison with conductance data.