Abstract
A number of parallel algorithms based on higher-order Padé approximations are developed for the numerical solution of a mathematical model of percutaneous drug absorption. These methods are L o-stable and require only the application of tridiagonal solvers with complex arithmetic.