/**
***变步长梯形求解微分方程***
获取初值: T1 = h/2[ f(x(k)) + f(x(k+1)) ]
n-1
变步长梯形公式: T2n = 1/2*Tn + h/2 * ∑ f ( x(k+1/2) )
k=0
步长: h=b-a/n
属性: 数值积分法
《数值分析简明教程》-2 Editon -高等教育出版社- page 67 算法流程图
代码维护:2005.6.14 DragonLord
**/
#include<iostream.h>
#include<math.h>
#include<stdio.h>
double f(double x)
{
double f;
if(x==0)f=1; // 1
else f=sin(x)/x; // 举例方程 I = ∫ sin(x)/x dx
return f; // 0
}