25 #include "floatmath.h"
30 extern void imcore_polynm (
float xdat[],
float xcor[],
int n,
float polycf[],
32 double a[25][25], b[25], temp;
39 for(i = 0; i < 25; i++) {
41 for(j = 0; j < 25; j++) a[i][j] = 0.0;
45 for(i = 0; i < n; i++) {
46 for(k = 0; k < m; k++) {
48 if(k+ilim != 0)temp = pow(xcor[i], (
float) (k+ilim));
51 for(j = 0; j <= k; j++) {
53 if(k+j+2*ilim != 0)temp = pow(xcor[i], (
float) (k+j+2*ilim));
59 for(k = 1; k < m; k++) {
60 for(j = 0; j < k; j++) a[k][j] = a[j][k];
64 imcore_solve(a, b, m);
66 for(i = 0; i < m; i++) polycf[i] = b[i];