13 void polynm (
float xdat[],
float xcor[],
int n,
float polycf[],
int m,
int ilim) {
14 double a[25][25], b[25], temp;
21 for(i = 0; i < 25; i++) {
23 for(j = 0; j < 25; j++) a[i][j] = 0.0;
27 for(i = 0; i < n; i++) {
28 for(k = 0; k < m; k++) {
30 if(k+ilim != 0)temp = pow(xcor[i], (
float) (k+ilim));
33 for(j = 0; j <= k; j++) {
35 if(k+j+2*ilim != 0)temp = pow(xcor[i], (
float) (k+j+2*ilim));
41 for(k = 1; k < m; k++) {
42 for(j = 0; j < k; j++) a[k][j] = a[j][k];
48 for(i = 0; i < m; i++) polycf[i] = b[i];