The evolution of unsteady boundary layers on oscillating airfoils is investigated by solving the governing equations by the Characteristic Box scheme. The difficulties associated with computing the first profile on a given time line, and the velocity profiles with partial flow reversal are solved. A sample calculation is performed for an external velocity distribution typical of those found near the leading edge of thin airfoils. The viability of the calculation procedure is demonstrated.