《计算方法》 李晓红等 第六章 数值积分 例题 代码

#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;