modelirovanie

#include <iostream>
#include <ctime>
#include <math.h>
#include <cmath>
#include <cstdlib>
using namespace std;
double f(double x)
{
    return x*x;
}
int main()
{
    double m;
    double h;
    double delta;
    double delta2;
    double N;
    double x[200];
    double f_del2[200];
    double A[200];
    cout << "m=";
    cin >> m;
    cout << endl;
    cout << "delta=";
    cin >> delta;
    cout << endl;
    cout << "delta2=";
    cin >> delta2;
    cout << endl;
    h = 1 / m;
    cout << "h=" << h << 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;
    system("PAUSE");
    return 0;
}