Esope source files

   1 : C NRQFUN    SOURCE    CHAT      05/01/13    02:01:54     5004           
   2 : c
   3 :              FUNCTION NR_QFUNC(QINI,QA,M,CT,PRECIS)
   4 :       IMPLICIT INTEGER(I-N)
   5 :       real*8 qini,qa,m,ct,precis,qqa,valor,dq,NR_qfunc
   6 :       integer it
   7 :       valor=qini
   8 :       it=0
   9 : 10    continue
  10 :       it=it+1
  11 :       qqa=(1.D0+valor/qa)
  12 :       dq=(valor*qqa**m-ct)/(qqa**m+m*valor/qa*qqa**(m-1.D0))
  13 :       valor=valor-dq
  14 :       if (abs(dq).lt.precis*(qini+1.D0)) then
  15 :          NR_qfunc=valor
  16 :          return
  17 :       endif
  18 :       if (it.gt.20) then
  19 :         write(*,*)' NR_qfunc error',it,abs(dq),precis
  20 :         stop
  21 :       endif
  22 :       goto 10
  23 :       end
  24 :  
  25 :  
  26 :