double precision z1, b, discri

         discri = z1*z1-dble(fp)*dble(fpa)

         if (discri .lt. 0.d0) then

            discri = dsqrt(discri)

         discri = dble(fp)/z1

         discri = discri*dble(fpa)

         discri = z1-discri

         if (z1.ge.0.d0 .and. discri.ge.0.d0) then

            discri = dsqrt(z1)*dsqrt(discri)

         else if (z1.le.0.d0 .and. discri.le.0.d0) then

            discri = dsqrt(-z1)*dsqrt(-discri)

         discri = -discri

         t    = t + fp*(ta-t)/sngl(b+discri)

         anum = sngl(b-discri)