转换成C代码 SUBROUTINE NOLSHEAEEZE31E42XXYYZZ IMPLICIT REAL8A-HO-Z DIMENSION EZ3XX4YY4ZZ4 DIMENSION E313E423 E311=XX3-XX1 E312=YY3-YY1 E313=Z

#include <math.h> void NOLSHEAE(double EZ[], double E31[], double E42[], double XX[], double YY[], double ZZ[]) { double DEX; E31[0] = XX[2] - XX[0]; E31[1] = YY[2] - YY[0]; E31[2] = ZZ[2] - ZZ[0]; E42[0] = XX[3] - XX[1]; E42[1] = YY[3] - YY[1]; E42[2] = ZZ[3] - ZZ[1]; EZ[0] = E31[1] * E42[2] - E42[1] * E31[2]; EZ[1] = E31[2] * E42[0] - E42[2] * E31[0]; EZ[2] = E31[0] * E42[1] - E42[0] * E31[1]; DEX = sqrt(EZ[0] * EZ[0] + EZ[1] * EZ[1] + EZ[2] * EZ[2] + 1.0e-80); EZ[0] = EZ[0] / DEX; EZ[1] = EZ[1] / DEX; EZ[2] = EZ[2] / DEX;