Аппроксимация

    x1.clear(); // лист с точками аппроксимации по абциссе 
    y1.clear(); // лист с точками аппроксимации по ординате
 
    int degree=10; // количество участков на которое делится первоначальная функция
    int partSize=x.size()/degree; // размер участка
 
    for (int p=0;p<degree;p++)
    {
        double diskret=1./1e4; // дискретизация аппроксимирующей функции
        int additional=2; // количество точек первоначальной функции, которое добавляется с каждой стороны дополнительно к основному отрезку
        int start=(p==0) ? p*partSize : p*partSize-additional;
        int end=(p==degree-1) ? x.size() : (p+1)*partSize+additional;
 
        QList <double> x_part(x.begin()+start,x.begin()+end); // лист с дополнителными точками
        QList <double> y_part(y.begin()+start,y.begin()+end); // лист с дополнителными точками
 
        ///////////// поиск аппроксимирующего полинома
        int n=x_part.size();
        gsl_matrix *X=gsl_matrix_alloc(n,6);
        gsl_vector *yv=gsl_vector_alloc(n);
        gsl_vector *c=gsl_vector_alloc(6);
        gsl_matrix *cov=gsl_matrix_alloc(6,6);
 
        for (int i=0;i<n;i++)
        {
            for (int j=0;j<=5;j++)
            {
                gsl_matrix_set(X,i,j,pow(x_part[i],j));
            }
            gsl_vector_set(yv,i,y_part[i]);
        }
        {
            double chisq;
            gsl_multifit_linear_workspace *work=gsl_multifit_linear_alloc(n,6);
            gsl_multifit_linear(X,yv,c,cov,&chisq,work);
            gsl_multifit_linear_free(work);
        }
        ///////////// 
 
        for (double i=x_part.first();i<=x_part.last();i+=diskret) // построение аппроксимации по полиному
        {
            double y_val=0;
            for (int j=0;j<=5;j++)
            {
                y_val+=gsl_vector_get(c,j)*pow(i,j);
            }
            x1.append(i);
            y1.append(y_val);
        }
 
        gsl_matrix_free(X);
        gsl_vector_free(yv);
        gsl_vector_free(c);
        gsl_matrix_free(cov);
    }