function fexp(x)
      external Factor
      real sum, x, a
      sum = 0
      do while(abs(sum-exp(x)) .gt. 10**(-5))
         do a =1, 50
         sum = sum+x**a/factor(x)
         enddo
      enddo
      end