A mathematical absorption model (e.g. transit compartment model) is useful to describe complex absorption process. However, in such a model, an assumption has to be made to introduce multiple doses that a prior dose has been absorbed nearly completely when the next dose is administered. This is because the drug input cannot be determined from drug depot compartment through integration of the differential equation system and has to be analytically calculated. We propose a method of dose superimposition to introduce multiple doses; thereby eliminating the assumption. The code for implementing the dose superimposition in WinNonlin and NONMEM was provided. For implementation in NONMEM, we discussed a special case (SC) and a general case (GC). In a SC, dose superimposition was implemented solely using NM-TRAN abbreviated code and the maximum number of the doses that can be administered for any subject must be pre-defined. In a GC, a user-supplied function (FUNCA) in FORTRAN code was defined to perform dose superimposition to remove the restriction that the maximum number of doses must be pre-defined.