One of the promising frontiers of bioengineering is the controlled release of a therapeutic drug from a vehicle across the skin (transdermal drug delivery). In order to study the complete process, a two-phase mathematical model describing the dynamics of a substance between two coupled media of different properties and dimensions is presented. A system of partial differential equations describes the diffusion and the binding/unbinding processes in both layers. Additional flux continuity at the interface and clearance conditions into systemic circulation are imposed.

We consider a class of integral equations of Volterra type with constant coefficients containing a logarithmic difference kernel. This equation can be transformed into an equivalent singular equation of Cauchy type which allows us to give the explicit formula for the solution. The numerical method proposed in this paper consists of applying the Lagrange interpolation to the inner Cauchy type singular integral in the latter formula after subtracting the singularity. For the error of this method weighted norm estimates as well as estimates on discrete subsets of knots are given.

We consider a class of integral equations of Volterra type with constant coefficients containing a logarithmic difference kernel. This class coincides for a=0 with the Symm's equation. We can transform the general integral equation into an equivalent singular equation of Cauchy type which allows us to give an explicit formula for the solution g. The numerical method proposed in this paper consists in substituting the Lagrange polynomial interpolating the known function f in the expression of the solution g.