Mine

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