翻译成C语言 SUBROUTINE PREXYBIXYDXXYYXAYANK IMPLICIT REAL8A-HO-Z DIMENSION XXNKYYNK DIMENSION E12E22E32E42 IXYD=0 E11=XX1-XA E12=YY1-YA E21=XX2-XA E22=YY

void PREXYB(int *IXYD, double XX[], double YY[], double XA, double YA, int NK) { double E1[2], E2[2], E3[2], E4[2]; *IXYD = 0; E1[0] = XX[0] - XA; E1[1] = YY[0] - YA; E2[0] = XX[1] - XA; E2[1] = YY[1] - YA; E3[0] = XX[2] - XA; E3[1] = YY[2] - YA; E4[0] = XX[3] - XA; E4[1] = YY[3] - YA; double Z12 = E1[0] * E2[1] - E2[0] * E1[1]; if (Z12 >= -1.0e-10) { double Z23 = E2[0] * E3[1] - E3[0] * E2[1]; if (Z23 >= -1.0e-10) { double Z31 = E3[0] * E1[1] - E1[0] * E3[1]; if (Z31 >= -1.0e-10) { *IXYD = 1; goto exit; } } } double Z13 = E1[0] * E3[1] - E3[0] * E1[1]; if (Z13 >= -1.0e-10) { double Z34 = E3[0] * E4[1] - E4[0] * E3[1]; if (Z34 >= -1.0e-10) { double Z41 = E4[0] * E1[1] - E1[0] * E4[1]; if (Z41 >= -1.0e-10) { *IXYD = 2; } } } exit: return;