C DVCAIN SOURCE CHAT 05/01/12 22:58:54 5004 c SUBROUTINE DVFLCAMACINV (SIGINV,VARINT,OTRVAR,XPP,XQQ,XQT,XTT) IMPLICIT INTEGER(I-N) real*8 siginv(3),xpp,xqq,xqt,xtt,q,qqa,fr,pcap,gtet, . gtetd,gtetdd,pcapfunc,etaconfunc, . wwfunc,wwfuncder,wwfuncderder,varint,otrvar real*8 pc0,exr,cca,phi,alp,pcc,paa,qaa,exm,eww, . cpr,cpm,cfr common /MRScapp/ pc0,exr,cca,phi,alp,pcc,paa,qaa,exm,eww, . cpr,cpm,cfr pcap =pcapfunc(varint) gtet =wwfunc(siginv(3),eww) gtetd =wwfuncder(siginv(3),eww) gtetdd=wwfuncderder(siginv(3),eww) fr =(cfr*pcap-pcc)*etaconfunc(otrvar) q =siginv(2) qqa =(1.D0+q/qaa) xpp=2.D0/(cpr*pcap)**2 xqq=(2.D0*gtet*gtet/fr/fr)*(qqa**(2.D0*exm)+ . exm*(2.D0*exm-1.D0)*q*q*qqa**(2.D0*exm-2.D0)/qaa/qaa+ . 4.D0*exm*q*qqa**(2.D0*exm-1.D0)/qaa) xqt=(4.D0*gtet*gtetd*q/fr/fr)* . (qqa**(2.D0*exm)+exm*q*qqa**(2.D0*exm-1.D0)/qaa) xtt=(2.D0*q*q*qqa**(2.D0*exm))/fr/fr*(gtetd**2+gtet*gtetdd) return end