Sources Esope | Cast3M

vycoin
C VYCOIN    SOURCE    CHAT      05/01/13    04:10:36     5004
 
      SUBROUTINE VYICOMACINV (SIGINV,VARINT,XP,XQ,XT)
      IMPLICIT INTEGER(I-N)
      real*8    siginv(3),varint,xp,xq,xt,aux1,aux2,
     .          wwfunc,etaconfunc,wwfuncder
        real*8    exm,qaa,pcc,eww,xk1,xk2,etb,exv,xep,paa,cco,
     .          exl,xnn,aaa,bbb
      common /MRScone/ exm,qaa,pcc,eww,xk1,xk2,etb,exv,xep,paa,cco,
     .                 exl,xnn,aaa,bbb
      aux1=wwfunc(siginv(3),eww)
      aux2=(1.D0+siginv(2)/qaa)
      xp=-etaconfunc(varint)
      xq=aux1*aux2**exm+aux1*exm*siginv(2)/qaa*aux2**(exm-1.D0)
      xt=siginv(2)*aux2**exm*wwfuncder(siginv(3),eww)
      return
      end