#include<stdio.h>
int main()
{
float va[]={50.0f,46.0f,40.0f,32.2f,22.5f},temp,a,b,h;
int j,i,n;
a=0;
b=8;
n=2;
h=(b-a)/n;
temp=1/va[0]+1/va[4];
//complex T
for(i=0;i<n-1;++i)
{
j=i+1;
temp=temp+2/va[j];
}
printf("complex T :%f\n",h*temp/2);
//complex simpson
h=4;
temp=1/va[0];
for(i=0;i<=n-1;++i)
{
j=i+1;
temp=temp+4/va[j];
j=j+1;
temp=temp+2/va[j];
}
printf("complex simpson :%f\n",h*temp/6);
return 0;
}
int main()
{
float va[]={50.0f,46.0f,40.0f,32.2f,22.5f},temp,a,b,h;
int j,i,n;
a=0;
b=8;
n=2;
h=(b-a)/n;
temp=1/va[0]+1/va[4];
//complex T
for(i=0;i<n-1;++i)
{
j=i+1;
temp=temp+2/va[j];
}
printf("complex T :%f\n",h*temp/2);
//complex simpson
h=4;
temp=1/va[0];
for(i=0;i<=n-1;++i)
{
j=i+1;
temp=temp+4/va[j];
j=j+1;
temp=temp+2/va[j];
}
printf("complex simpson :%f\n",h*temp/6);
return 0;
}
本文介绍了一段使用C语言实现的数值积分计算代码,通过复化梯形法则和复化辛普森法则来近似计算特定函数的积分。首先定义了待积分的函数形式,并使用两种不同的数值积分方法进行计算。通过对各节点函数值的加权求和,最终得到积分的近似值。
1489

被折叠的 条评论
为什么被折叠?



