翻译成C代码 SUBROUTINE EQUISTRESGTSG0ANGLER00R45R90 IMPLICIT REAL8A-HO-Z DIMENSION SG03 RG=10D0R00 RF=10D0R90 RN=R45+05D010D0R00+10D0R90 SGT=RG+1D0SG01SG01+RF+1D0

void EQUISTRE(double *SGT, double *SG0, double ANGLE, double R00, double R45, double R90){ double RG = 1.0 / R00; double RF = 1.0 / R90; double RN = (R45 + 0.5) * (1.0 / R00 + 1.0 / R90); *SGT = (RG + 1.0) * SG0[0] * SG0[0] + (RF + 1.0) * SG0[1] * SG0[1] - 2.0 * SG0[0] * SG0[1] + 2.0 * RN * SG0[2] * SG0[2]; *SGT = sqrt(1.5 / (RF + RG + 1.0) * (*SGT)); }