Modelirovanie i sglazhivanie

#include <iostream>
#include <ctime>
#include <math.h>
#include <cmath>
#include <cstdlib>
using namespace std;

double f(double x)
{
    return x*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 t[65000];
    double T_alpha;
    double delta;
    double delta2;
    double N;
    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 = delta2;
    n=m;
    h1=(2*alpha)/m;
    cout<<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;
    }
    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) ))
            {
                  //double j = x[i] - alpha;
                  //int tt = 0;
                  int g=0;
                  double a=x[i]-alpha;
                  cout<<"a"<<a<<endl;
                  double b=x[i]+alpha;
                  cout<<"b"<<b<<endl;
                  while (a<=b)
                  {
                  t[g]=a;
                  a+=h1;
                  g++;
                  }
                 /* while(j <= x[i] + alpha) {
                          t[tt] = j;
                          tt++;
                          }*/
            } else {
                   cout << "Ne vhodit" << 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() % 10;
        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; i++) {
            T_alpha+= 3 / (2 * pow(alpha,3)) * ( (t[i] - x[i]) * f_del2[i]);
            cout << "i = " << i << " T_alpha = " << T_alpha << endl;
            }
    system("PAUSE");
    return 0;
}