#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;
}