Romberg积分法计算定积分（C语言）

/*
-------------Romberg积分法计算定积分-------
参考教材：《数值分析》李乃成，梅立泉，科学出版社
    《计算方法教程》第二版 凌永祥，陈明逵*/

//误差界eps%被积函数为f(x)=(x^3+sin(x))/x;积分区间为[0.3,0.8]

#include<stdio.h>
#include<stdlib.h>
#include<math.h>

int main(void)
{
/* int i,j,n;

 float *d=(float *)malloc(sizeof(float)*N);
 float *x=(float *)malloc(sizeof(float)*N);
 float *y=(float *)malloc(sizeof(float)*N);
 float *u=(float *)malloc(sizeof(float)*N);

 free(a);free(b);free(c);free(d);free(x);free(y);free(l);
 
 return 0;
*/

 double b=0.8;
 double a=0.3;
 double h=0.0; //
 double eps=1.0e-5;//误差界eps
 int kmax=20; //最大递推次数
 double T1=0.0,S1=0.0,C1=0.0,R1=0.0,T2=0.0,S2=0.0,C2=0.0,R2=0.0;
 double sum;
 double *x,*fx;
 int i;
 h=b-a;
 T1=h/2*((pow(a,3)+sin(a))