翻译成C代码 SUBROUTINE FORMJBBBPCCPAJFDAJXTYT IMPLICIT REAL8A-HO-Z DIMENSION BBP4CCP4AJF22XT4YT4 DIMENSION AJ22 AJ11=025D0-XT2-XT3+XT4 AJ12=025D0-YT3+YT4 AJ21=025D0XT

void FORMJB(double BBP[4], double CCP[4], double AJF[2][2], double* DAJ, double XT[4], double YT[4]){ double AJ[2][2]; AJ[0][0] = 0.25*(-XT[1]-XT[2]+XT[3]); AJ[0][1] = 0.25*(-YT[2]+YT[3]); AJ[1][0] = 0.25*(XT[1]-XT[2]-XT[3]); AJ[1][1] = 0.25*(-YT[2]-YT[3]); DAJ = AJ[0][0]AJ[1][1]-AJ[1][0]AJ[0][1]; AJF[0][0] = AJ[1][1]/(DAJ); AJF[0][1] = -AJ[0][1]/(DAJ); AJF[1][0] = -AJ[1][0]/(DAJ); AJF[1][1] = AJ[0][0]/(DAJ); BBP[0] = 0.25(AJF[0][0]+AJF[0][1]); BBP[1] = 0.25(-AJF[0][0]+AJF[0][1]); BBP[2] = 0.25(-AJF[0][0]-AJF[0][1]); BBP[3] = 0.25(AJF[0][0]-AJF[0][1]); CCP[0] = 0.25(AJF[1][0]+AJF[1][1]); CCP[1] = 0.25(-AJF[1][0]+AJF[1][1]); CCP[2] = 0.25(-AJF[1][0]-AJF[1][1]); CCP[3] = 0.25*(AJF[1][0]-AJF[1][1]);