#define _CRT_SECURE_NO_WARNINGS

#include <stdlib.h>
#include <stdio.h>
#include <time.h>
#include <math.h>

#define N 2048
#define PI 3.14159265

//funcao para obter f(a)
double f(double q);

void main()
{
    int i, j, n;
    float r = 0;
    double M[N][3], x, y, Z, inc, a, b;
    FILE *fp;

    //definindo a primeira coluna da matriz
    srand(time(NULL));
    for (i = 0; i < N; i++)
    {
        M[i][0] = r;
        r = r + 0.05;
    }

    //obtendo segunda coluna da matriz
    for (i = 0; i < N; i++)
    {
        x = M[i][0];
        y = f(x); //executando funcao p/ obter f(a)
        M[i][1] = y;
    }

    //obtendo terceira coluna da matriz
    for (i = 0; i < N; i++)
    {
        b = M[i][0]; //definindo a
        a = -1 * b; //definindo -a
        n = 1000;
        inc = (b - a) / n;
        Z = 0;
        //obtendo I(a) pela regra dos retangulos
        for (j = 0; j < n; j++)
        {
            b = a + inc;
            Z = Z + (b - a)*(f(a) + f(b)) / 2;
            a = b;
        }
        M[i][2] = Z;
    }

    //printando matriz
    for (i = 0; i < N; i++)
    {
        printf("%f\n%f\n%f\n\n\n", M[i][0], M[i][1], M[i][2]);
    }

    //criando arquivo
    fp = fopen("ra14194054.dat", "w");
    for (i = 0; i<N; i++)
    {
        fprintf(fp, "%f ", M[i][0]);
        fprintf(fp, "  %f ", M[i][1]);
        fprintf(fp, "  %f ", M[i][2]);
        fprintf(fp, "\n");
    }
    fclose(fp);
    system("pause");
}

double f(double q)
{
    //obtendo f(a)
    double y;
    y = sin((PI * q)* PI / 180);
    y = y / (PI * q);
    return y;
}