C语言实现龙贝格求积

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


double f(double x)
{
        return (4*x*x/(1+x*x));
}


double Romberg(double top,double bottom,double precision)
{
        int k=1;
        double S,x,T1,T2,S1,S2,C1,C2,R1,R2,h=bottom-top;


/*
转载请注明出处：去转盘网www.quzhuanpan.com
*/


        //S梯形公式，T梯型变步长，S梯形加速，C幸普森加速，R龙贝格求积
        T1=h*(f(top)+f(bottom))/2;//梯形公式
        while(1)
        {
                S=0;
                x=top+h/2;
                do
                {