A novel semi-analytical solution method is proposed for 1D multilayer diffusion.
Method based on the Laplace transform and an orthogonal eigenfunction expansion.
Method avoids solving complex transcendental equations for the eigenvalues.
Method is accurate and efficient for a small or a large number of layers.
An application to macroscopic modelling of multilayer diffusion is explored.