Abstract: The creep of a skin layer under a distributed surface pressure was solved by ananalysical method using Hankel transform and Laplace transform.The surface stressboundary conditions lead io a Volterra integral equation of the first kind, which was then solved by a numerical method.The IMSL subroutines DINLAP and DGORUL were employed to numerically obtain the Hankel-Laplace inversion. The calculateddisplacements at two distinctive moments were compared respectively with those obtained by an elastic solution for either incompressible or compressible solid. Thetransient creep responses of the skin layer were also presented.

Key words: numerical method, computation of shock waves, upwind compact scheme