Моделирование функции заданной с погрешностью и сглаживание выбросов

#include <iostream>
#include <ctime>
#include <math.h>
#include <cmath>
#include <cstdlib>
using namespace std;
 
double f(double x)
{
    return 2*x;
}       
int pow(int a, int b) {
    int result = 1;
    for (int i = 0; i < b; i++)
        result *= a; 
    return result;
    }
int main()
{
    double m;
    double h;
    double sum;
    double t[65000];
    double T_alpha[200];
    double delta;
    double delta2;
    double N;
    double gamma;
    double x[200];
    double f_del2[200];
    double A[200];
    double alpha;
    double n;
    double h1;
    cout << "m=";
    cin >> m;
    cout << endl;
    cout << "delta=";
    cin >> delta;
    cout << endl;
    cout << "delta2=";
    cin >> delta2;
    cout << endl;
    h = 1 / m;
    alpha=pow(0.75, 0.25)*pow(delta, 0.25);
    n=m;
    h1=(2*alpha)/m;
    cout<<"h1="<<h1<<endl;
    cout << "h=" << h << endl;
    cout << "alpha=" << alpha << endl;
    for (int i = 0; i < m; i++)
    {
        x[i] = h*i;
        cout << "x[" << i << "]=" << x[i] << endl;
    }
    N = sqrt((0, 5 * (m*pow(delta / delta2, 2) - 1)) + 1);
    cout << "N=" << N << endl;
    for (int i = 0; i < m; i++)
    {
        A[i] = 1 + rand() % 100;
        cout << "A[" << i << "]=" << A[i] << endl;
    }
    for (int i = 0; i < m; i++)
        f_del2[i] = f(x[i]) + pow(1, i) * A[i] * delta2;
    for (int i = 4; i < m; i++)
        f_del2[i] = f(x[i]) + pow(-1, i) * A[i] * delta2;
    f_del2[2] = f(x[2]) + N*A[2] * delta2;
    f_del2[3] = f(x[3]) - N*A[3] * delta2;  
    for (int i = 0; i < m; i++)
        cout << "f(x[" << i << "])=" << f(x[i]) << endl;
    for (int i = 0; i < m; i++)
        cout << "f_del2[" << i << "]=" << f_del2[i] << endl;
    for(int i = 0; i < n - 1; i++) {
            if (( ((x[i] - alpha) >= 0) && ((x[i] - alpha) <= 1) ) && ( ((x[i] + alpha) >= 0) && ((x[i] + alpha) <= 1) )) 
            {
                  int g=0;
                  double a=x[i]-alpha;                                
                  //cout << "a "<< a << endl;
                  double b=x[i]+alpha;                                
                  //cout << "b "<< b << endl;
                  sum=0;
                  while (a<=b)                                         
                  {                                                   
                  t[g]=a;                                             
                  f_del2[g] = f(t[g]) + pow(1, g) * A[i] * delta2;
                  sum += (3 / (2 * pow(alpha, 3))) * ( (t[g] - x[i]) * f_del2[g] );
                  //cout<<3 / (2 * pow(alpha,3))<<" umnozhenie "<<g<<endl;
                  //cout<<(t[g] - x[i])<<" skobka "<<g<<endl;
                  //cout<<f_del2[g]<<" fdelta "<<g<<endl;
                  //cout<<t[g]<<endl;
                  a+=h1;
                  g++;
                  }
                  T_alpha[i]=sum;
                  cout << "i = " << i << " T_alpha = " << T_alpha[i] << endl;
                } 
            }
    system("PAUSE");
    return 0;
}