Среднеквадратичное отклонение

#include <iostream>
#include <conio.h>
#include <math.h>
using namespace std;
long double fdiff(long double x, long double y, int k)
{
    return k*sqrt(y)+sqrt(x);
}
void masX(long double* x, long double a, long double h, int n)
{
    for(int i=0; i<=n; i++)
        x[i]=a+i*h;
}
void masY(long double* x, long double* y, long double h, int n, int k)
{
    y[0]=1;
    for(int i=0; i<n; i++)
        y[i+1]=y[i]+h*fdiff(x[i], y[i], k);
}
void approx(long double* x, long double* y, int k, long double& c0, long double& c1, long double& c2)
{
    int n=3, m=4;
    long double a[3][4]={{0,0,0,0},{0,0,0,0},{0,0,0,0}},
                a0[3][4]={{0,0,0,0},{0,0,0,0},{0,0,0,0}};
    a0[0][0]=a[0][0]=k+1;
    for(int i=0; i<=k; i++)
        a0[0][1]=a0[1][0]=a[0][1]=a[1][0]=a[1][0]+x[i];
    for(int i=0; i<=k; i++)
        a0[0][2]=a0[1][1]=a0[2][0]=a[0][2]=a[1][1]=a[2][0]=a[2][0]+pow(x[i],2);
    for(int i=0; i<=k; i++)
        a0[1][2]=a0[2][1]=a[1][2]=a[2][1]=a[2][1]+pow(x[i],3);
    for(int i=0; i<=k; i++)
        a0[2][2]=a[2][2]=a[2][2]+pow(x[i],4);
    for(int i=0; i<=k; i++)
        a0[0][3]=a[0][3]=a[0][3]+y[i];
    for(int i=0; i<=k; i++)
        a0[1][3]=a[1][3]=a[1][3]+x[i]*y[i];
    for(int i=0; i<=k; i++)
        a0[2][3]=a[2][3]=a[2][3]+pow(x[i],2)*y[i];
    int i,j,k;
    float buf,x[3];
    for(i=0;i<n-1;i++)
    for(j=i+1;j<n;j++)
    {
        buf=a[i][i]/a[j][i];
        for(k=0;k<=n;k++) a[j][k]=a[j][k]*buf-a[i][k];
    }
    x[n-1]=a[n-1][n]/a[n-1][n-1];
    for(i=n-2;i>=0;i--)
    {
        buf=0;
        for(j=i+1;j<n;j++)
            buf+=a[i][j]*x[j];
        x[i]=(a[i][n]-buf)/a[i][i];
    }
    c0=x[0];
    c1=x[1];
    c2=x[2];
    cout<<endl<<"Решение (полином второй степени)"<<endl;
    cout<<"P2(x) = "<<x[2]<<" * x^2 + "<<x[1]<<" * x + "<<x[0]<<endl;
}
long double sko(long double* y, long double* yt, int n, long double* x)
{
    long double summ=0;
    return summ;
}
void main()
{
    long double h=0.1, a=0, b=1, c0, c1, c2, q=0;
    int n=(b-a)/h, k=1;
    setlocale(LC_ALL, "rus");
    long double x[11], y[11], yt[11];
    masX(x, a, h, n);
    masY(x, y, h, n, k);
    for(int i=0; i<=n; i++)
        cout<<endl<<"x= "<<x[i]<<"\t"<<"y= "<<y[i]<<endl;
    approx(x, y, n, c0, c1, c2);
    q=sko(y, yt, n,x);
    cout<<"q= "<<q<<endl;
    k=-1;
    _getch();
}