收集的C语言算法程序

Simpson公式法:

#include<iostream.h>
//用宏表示2*n
#define N (2*8)
//声明被积分的函数f(x)
inline double f(double);
void main()
{
 //hh表示公式中h的一半
 double a=0,b=1,sum=0,hh=(b-a)/N;
 double x=a;
 for(int i=1;i<N;i++)
 {
  //以hh的大小为步长递增
     x+=hh;
  //对i的单双数进行不同处理
  if(i%2==1)
     sum+=4*f(x);
  else
     sum+=2*f(x);
 }
 sum+=(f(a)+f(b));
 sum*=hh/3;
 //输出运行结果
  cout<<"经用Simpson公式计算知：\n原函数在【0，1】上的积分近似值为："<<sum<<endl;
}
//定义被积分的函数f(x)
inline double f(double x)
{
 return x/(4+x*x);
}


梯形公式法：
#include<iostream.h>
//用宏表示n
#define N 8
//声明被积分的函数f(x)
inline double f(double);
void main()
{
 double a=0,b=1,sum=0,h=(b-a)/N;
 double x=a;
 for(int i=1;i<N;i++)
 {
  //以h的大小为步长递增
  x+=h;
  sum+=f(x);
 }
 sum*=2;
 sum+=(f(a)+f(b));
 sum*=(h/2);
 //输出运行结果
 cout<<"经用梯形公式计算知：\n原函数在【0，1】上的积分近似值为："<<sum<<endl;
}
//定义被积分的函数f(x)  
inline double f(double x)
{
 return x/(4+x*x);
}