1 typedef float rvec[3]; calc_similar_ind(int nind,rvec x[])2float calc_similar_ind (int nind, rvec x[]) 3 { 4 int i, j, d; 5 float m, tm, xd, rd; 6 for(j=0; j<nind; j++) { 7 i = x[j][0]; 8 tm += m; 9 for(d=0 ; d<3; d++) { 10 xd = x[i][d] - x[i][d]; 11 rd += m * xd; 12 } 13 } 14 return rd/tm; 15 } 16