定步长复化梯形求解积分方程

/**
***定步长复化梯形求解积分方程***
                                                     n-1                                   n-1
  定步长复化梯形公式： Tn=   ∑  h/2〔 f(a) + f(b) + 2 ∑ f(x(k)) 〕
                                                    k=0                                   k=0
   属性：数值积分法
                                       n-1
   误差(精度)：I - Tn =  ∑  [ - h*h/12 f''(ε(k)) ] ≈ - h*h/12 [ f'(b) - f'(a) ]
                                      k=0 
   步长：h = ( b - a ) / n

  《数值分析简明教程》-2 Editon -高等教育出版社- page 64 无算法流程图
  
  代码维护：2005.6.14  DragonLord

**/

#include<iostream.h>
#include<math.h>
#include<stdio.h>
double f(double x)
{
    double f;
 if(x==0)f=1;
 else f=sin(x)/x;