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();
}

Tidak ada komentar