June 29
Simpson Formula
#include <iostream>
#include <stdlib.h>
#include <math.h>
using namespace std;
double c[100];
double Poly(double x,int n){
int i;
double ans=0;
for(i=n;i>=0;i--)
if (c[i]==0);else ans+=pow(x,i)*c[i];
return ans;
}
int main()
{
int n;
int i,j;
double k,a,b;
double s=0;
cin>>n;
for(i=n;i>=0;i--)cin>>c[i];
cin>>a>>b;
k=(b-a)/10000;
for(i=1;i<10000;i++)s+=Poly(a+double(i)*k,n);
s*=2;
s+=Poly(a,n);s+=Poly(b,n);
s=s*k*0.5;
cout<<s;
cin>>a;
}