Télécharger nrqfun.eso

Retour à la liste

Numérotation des lignes :

  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.  
  27.  

© Cast3M 2003 - Tous droits réservés.
Mentions légales