Pract: Develop a C program to implement Simpsons 1/3rd Rule.

#include<stdio.h>
void main()
{
int i,n;
float x[20], y[20], h, sum, integral, s2, s4;
clrscr();
printf(“Enter number of intervals: “);
scanf(“%d”,&n);
printf(“Enter size(height) of interval: “);
scanf(“%f”,&h);
printf(“Enter %d pair of (x,y)\n”,n+1);

for(i=1; i<=(n+1); i++) {
printf(“x[%d]= “,i);
scanf(“%f”,&x[i]);

printf(“y[%d]= “,i);
scanf(“%f”,&y[i]);
}
sum = y[1]+y[n+1];
s2 = s4 = 0.0;
for(i=2; i<=n; i+=2) {
s4 = s4+y[i];
}
for(i=3; i<=(n-1); i+=2) {
s2 = s2+y[i];
}
integral = (h/3.0)*(sum+2*s2+4*s4);
printf(“\nValue of Integral = %.4f\n”,integral);
getch();
}