#include "stdafx.h"
#include <math.h>
// 数值积分
// Newton-Cotes积分公式
void
NewtonCotesIntegral()
{
printf
(
"牛顿-科茨积分公式:\n"
);
int
n;
double
h;
int
i;
double
a=0,b=0.5;
double
I;
double
C[5],f[5],x[5];
for
(n=1;n<=4;n*=2)
{
h = (b-a)/n;