RG=10D0R00 RF=10D0R90 RN=R45+05D010D0R00+10D0R90 FGF=RFRG+RF+RG FAA=RGEPN2-EPZFGF GAA=EPZ-RFEPN1FGF HAA=RFEPN1-RGEPN2FGF ANN=20D0EPN3EPN3RN EPT=RF+RG+10D0RFFAAF

*EP0,*EPN,*EPZ,*ANGLE,*R00,*R45,*R90 EP0=DSQRT(2.0D0EPT/3.0D0) EPN=EPNANGLE/PI EPZ=EPZANGLE/PI *R00=1.0D0/R00 *R90=1.0D0/R90 *R45=R45+0.5D0 *RN=R45(*R00+*R90) FGF=*R90**R00+*R90+R00 FAA=(R00EPN(2)-EPZ)/FGF GAA=(EPZ-R90EPN(1))/FGF HAA=(R90EPN(1)-R00EPN(2))/FGF ANN=2.0D0EPN(3)*EPN(3)/(R45(R00+R90)) EPT=(R90+R00+1.0D0)(R90FAAFAA+R00GAAGAA+HAAHAA+ANN)