double precision z1, b, discri

      z1 = dble(fp) + dble(fpa) - 3.d0*dble(fa-f)/dble(ta-t)

      b  = z1 + dble(fp)

      if (abs(z1).le.1.) then

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

         discri = dble(fp)/z1

         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)

         den  = sngl(z1+b+dble(fpa))