Simpson kuralı belirli integrali bulmak için bir yaklaşım yöntemidir. Tıpkı trapezoid kurali gibi veya orta nokta metodu gibi. Bu gibi yöntemlerde örnek sayısı ne kadar arttırılırsa o kadar iyi sonuç alınır.Ancak bu bir yaklaşım olduğuna göre makul sayıdaki örnek ile en doğru sonuca ulaşılması gerekir. Simpson metodu diğer iki metoddan daha doğru sonuç vermektedir.
Bu kod aşşağıdaki başlıkların örneğini içermektedir;
- Dizi
- For döngüsü
- If döngüsü
Program öncelikle kaç tane değer giriceğini soruyor. Sonrasında bu değerlerin x ve fonksiyonda karşılığı na denk gelen y değerlerini alıp sonucu bildiriyor.
Kod;
/*****************************************************
* Numerical Integration by Simpson (1/3) rule
* (x,y) points are given as inputs
*****************************************************/
#include < stdio.h>
#include < math.h>
void main()
{
float x[10],y[10],sum=0,h;
int i,n,k=0;
printf(“\nhow many record you will be enter: “);
scanf(“%d”,&n);
for(i=0; i< n; i++)
{
printf(“\n\nenter the value of x%d: “,i);
scanf(“%f”,&x[i]);
printf(“\n\nenter the value of f(x%d): “,i);
scanf(“%f”,&y[i]);
}
h=x[1]-x[0];
n=n-1;
sum = sum + y[0];
for(i=1;i< n;i++)
{
if(k==0)
{
sum = sum + 4 * y[i];
k=1;
}
else
{
sum = sum + 2 * y[i];
k=0;
}
}
sum = sum + y[i];
sum = sum * (h/3);
printf(“\n\n I = %f “,sum);
system(“PAUSE”);
}
