ccccccccccccccccccccccccc     Program 2.7     cccccccccccccccccccccccccc
c
cccccccccccccccccccccccccccccccccccccccccccccccccccccccccccccccccccccccc
c                                                                      c
c Please Note:                                                         c
c                                                                      c
c (1) This computer program is part of the book, "An Introduction to   c
c     Computational Physics," written by Tao Pang and published and    c
c     copyrighted by Cambridge University Press in 1997.               c
c                                                                      c
c (2) No warranties, express or implied, are made for this program.    c
c                                                                      c
cccccccccccccccccccccccccccccccccccccccccccccccccccccccccccccccccccccccc
c
      PROGRAM NEWTON
C
C This program uses the Newton method to find the root of
C f(x)=exp(x)*ln(x)-x*x=0.
C
      DL = 1.0E-06
      A  = 1.0
      B  = 2.0
      DX = B - A
      X0 = (A+B)/2.0
      ISTEP = 0
      DO    100  WHILE (ABS(DX).GT.DL)
        X1 = X0 - F(X0)/DF(X0)
        DX = X1 - X0
        X0 = X1
        ISTEP = ISTEP + 1
  100 END DO
      WRITE (6,999) ISTEP,X0,DX
      STOP
  999 FORMAT (I4,2F16.8)
      END
C
      FUNCTION F(X)
        F = EXP(X)*ALOG(X) - X*X
      RETURN
      END
C
      FUNCTION DF(X)
        DF = EXP(X)*(ALOG(X) + 1./X) - 2.*X
      RETURN
      END