Metode Simpson

Algoritma :

1.    Definisikan fungsi f(x)

2.    Tentukan batas bawah dan batas atas integrasi

3.    Tentukan jumlah pembagi area N

4.    Hitung h=(b-a)/N

5.    Hitung V = phi.(h/3).( fo + 4 ∑ fi(ganjil) + 2 ∑ fi(genap) + fn )

Program :

#include<stdio.h> #include<math.h> double f(double x, double r){                return (sqrt(r*r - x*x) - 4.5)*(sqrt(r*r - x*x) - 4.5); } main(){                int bil=0;                double a, b, h, N, V, sigma=0, ganjil=0, genap=0, r=12.5, x, phi=3.142857143;                printf("\n Menghitung Volume");                printf("\n Integrasi Simpson\n\n");                printf(" batas bawah : "); scanf("%lf",&b);                printf(" batas atas  : "); scanf("%lf",&a);                printf(" pembagi (N) : "); scanf("%lf",&N);                h=(a-b)/N;                for(x=b;x<=a;x+=h){                                if(bil>0 && bil<N){                                                if(bil%2==0){                                                                genap+=f(x,r);                                                } else{                                                                ganjil+=f(x,r);                                                }                                }                                bil++;                }                V=phi*(h/3)*(f(b,r)+(2*genap)+(4*ganjil)+f(a,r));                printf("\n");                printf(" Volume_simpson = %lf \n",V);                fflush(stdin);                getch(); }