Sinxdx

#include <iostream>
#include <cstdlib>
#include <cmath>

#define M_PI 3.141592

/* calkowanie numeryczne funkcji sin(x) metoda Simpsona */

int main( void )
{
    double zakres = 3.0;
    double h = 0.1;

    int n = 30;

    double* x = new double[n];
    double* f = new double[n];

    std::cout << n << std::endl;

    for (int i = 0; i <= n; i++)
    {
        x[i] = i*h;
    }

    for (int i = 0; i <= n; i++)
    {
        f[i] = sin(x[i]);
    }

    for (int i = 0; i <= n; i++)
    {
        std::cout << "x=" << x[i] << "  f= " << f[i] << std::endl;
    }

    double calka = 0;

    double c1 = f[0] + f[n];
    double c2 = 0;

    for (int i = 0; i <= n - 2; i += 2)
    {
        c2 += f[i];
    }
    c2 *= 2;

    double c3 = 0;
    for (int i = 1; i <= n - 1; i += 2)
    {
        c3 += f[i];
    }
    c3 *= 4;

    calka = (c1 + c2 + c3) * 0.0333333333333;

    std::cout << "calka sinxdx od 0 do 3 = " << calka << std::endl;

    system("pause");


    return 0;
}