Interpolação Polinomial em C

Um algoritimo fácil para interpolação polinomial linear em C
Cálculo Numérico
Os valores são lidos nos vetores x[i] e fx[i]. Nesse caso, usei com 10 combinações lineares dessas(10 pontos)

     somax = 0;//inicializando os valores dos somatorios.
     somax2 = 0;
     somaf = 0;
     somaxf = 0;
     n=10;
     for(i=1;i<=10;i++){//Somatorios de x, x², f(x) e x*f(x)
                        somax = somax + x[i];
                        somax2 = somax2 + x[i]* x[i];
                        somaf = somaf + fx[i];
                        somaxf = somaxf + x[i]* fx[i];
                        }              
                      
                       //calculo dos valores de a e b da regressão linear.
                       a = (somaf * somax2 - somax * somaxf)/(n * somax2 - somax *somax);
                       b = (n * somaxf - somax * somaf)/(n * somax2 - somax * somax);
                      
                       printf ("\nA curva que melhor ajusta os 10 dados fornecidos eh \ng(x)= %f + %f*x\n\n", a, b);