/* Dasproc C - 2019
William Handi Wijaya
0087

Memperkirakan luas suatu fungsi pada interval [a,b]
*/

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

double func_g(double x){ //fungsi g
    return x*x*sin(x);
}

double func_f(double x){ //fungsi f
    return sqrt(4-(x*x));
}

void trap(double a, double b, double n, double f(double)){ //fungsi pemanggilan trap
    //deklarasi dan inisiasi
    double x,h,hp = a,i,T;
    x=0;
    h=(b-a)/n;
    //penghitungan nilai zigma
    for(i=1;i<n;i++){
        hp +=h;
        x += f(hp);
    }
    //penghitungan nilai T
    T = (h/2) * (f(a) + f(b) + (2*x));
    //menampilkan nilai T
    printf("\tT = %lf\n", T);
}

int main(){
    //deklarasi variabel
    int n;

    //penghitungan dan penampilan untuk fungsi g
    printf("for x^2 sin(x)\n");
    for(n=2;n<=128;n*=2){
        printf("\tfor n = %d,", n);
        trap(0,3.14159,n,func_g); // a = 0, b = 3.14159
    }
    //penghitungan dan penampilan untuk fungsi f
    printf("for sqrt(4-(x^2))\n");
    for(n=2;n<=128;n*=2){
        printf("\tfor n = %d,", n);
        trap(-2,2,n,func_f); // a = -2, b = 2
    }

    return 0;
}