C DVCOIN SOURCE CHAT 05/01/12 22:58:58 5004 c SUBROUTINE DVFLCOMACINV (SIGINV,VARINT,OTRVAR,XPP,XQQ,XQT,XTT) IMPLICIT INTEGER(I-N) real*8 siginv(3),xpp,xqq,xqt,xtt,auxq,wwd0,wwd1,wwd2,q, . wwfunc,wwfuncder,wwfuncderder,pcap,pcapfunc,p,ap, . varint,otrvar,etaconfunc real*8 exm,qaa,pcc,eww,xk1,xk2,etb,exv,xep,paa,cco, . exl,xnn,aaa,bbb,pc0,exr,cca,phi,alp,pcc2,paa2,qaa2, . exm2,eww2,cpr,cpm,cfr common /MRScone/ exm,qaa,pcc,eww,xk1,xk2,etb,exv,xep,paa,cco, . exl,xnn,aaa,bbb common /MRScapp/ pc0,exr,cca,phi,alp,pcc2,paa2,qaa2,exm2,eww2, . cpr,cpm,cfr wwd0=wwfunc(siginv(3),eww) wwd1=wwfuncder(siginv(3),eww) wwd2=wwfuncderder(siginv(3),eww) pcap=pcapfunc(otrvar) p=siginv(1) ap=alp*pcap q=siginv(2) auxq=(1.D0+q/qaa) xpp=abs(xnn)*etaconfunc(varint)*(2.D0*ap-2.D0*pcc)/ . (p+ap-2.D0*pcc)**2 xqq=(wwd0*(exm-1.D0)*exm*q*auxq**(exm-2.D0))/(qaa*qaa)+ . (2.D0*wwd0*exm*auxq**(exm-1.D0))/qaa xqt=wwd1*auxq**exm+(wwd1*exm*q*auxq**(exm-1.D0))/qaa xtt=q*auxq**exm*wwd2 return end