高分求数值积分两个C程序

来源：百度知道 编辑：UC知道 时间：2024/06/05 19:44:18
(1) 利用Romborg算法计算一个给定函数定积分的近似值。
(2) 利用复合梯形公式和复合Simpson公式计算一个给定函数定积分的近似值。
给定函数f(x)=ex(x为e指数)cosx ，输入求积区间 和误差精度要求后，利用Romborg算法或复合梯形公式和复合Simpson公式计算积分f(x)(上下限为a.b)的近似值。
输入a，b，tol.
急用.调试通过的.再加100分
补充一下.【输入样例】
0 3.14159 0.0001
【输出样例】
1.0986
三楼的不很符合要求

用Romborg的算法为：
#include <stdio.h>
#include <conio.h>
#include <math.h>
float f(float x)
{
return exp(x)*cos(x);
}
float Romberg(float a,float b,float (*f)(float),float epsilon)
{
int n=1,k;
float h=b-a,x,temp;
float T1,T2,S1,S2,C1,C2,R1,R2;
T1=(b-a)/2*((*f)(a)+(*f)(b));
while(1)
{
temp=0;
for(k=0;k<=n-1;k++)
{
x=a+k*h+h/2;
temp+=(*f)(x);
}
T2=(T1+temp*h)/2;
if(fabs(T2-T1)<epsilon)return T2;
S2=T2+(T2-T1)/3.0;
if(n==1){T1=T2;S1=S2;h/=2;n*=2;continue;}
C2=S2+(S2-S1)/15;
if(n==2){C1=C2;T1=T2;S1=S2;h/=2;n*=2;continue;}
R2=C2+(C2-C1)/63;
if(n==4){R1=R2;C1=C2;T1=T2;S1=S2;h/=2;n*=2;continue;}
if(fabs(R2-R1)<epsilon)return R2;
R1=R2;C1=C2;T1=T2;h/=2;n*=2;
}
}
main()
{
float epsilon=5e-6;
printf("R=%f&qu