24小时热门版块排行榜

返回列表

zhangjin_7

铜虫 (小有名气)

应助: 1 (幼儿园)
金币: 207.3
散金: 1155
红花: 1
帖子: 257
在线: 99.2小时
虫号: 2087044
注册: 2012-10-26
专业: 化学计量学与化学信息学

[求助] 求助romberg调试

#include
#include
#define jifen(x) pow(3,x)*pow(x,1.4)*(5*x+7)*sin(pow(x,2))//定义待求积分函数
#define E 0.00001//定义控制迭代误差大小
main()
{
double newarray[20]={0},oldarray[20]={0};//定义两个数列分别表示迭代是用到的两组数据
double fistarray(double a,double b,int j);//声明一个函数计算新数组的第一个元素
int i=1,j,k;//定义i为新数组中元素个数,j控制将新数组赋值给旧数组，k控制计算新数组
oldarray[0]=jifen(1)+jifen(3);//初始化旧数组中第一个元素
newarray[0]=(3-1)/2*(jifen(1)+jifen(3)-2*jifen(2));//初始化新数组中第一个元素
newarray[1]=(4*newarray[0]-oldarray[0])/(4-1);//初始化新数组中第二个元素
while(fabs(newarray-oldarray[i-1])>=E)
{
for(j=0;j<=i;j++)//将新数组数据赋给旧数组
oldarray[j]=newarray[j];
i++;//新数组长度增加一位
newarray[0]=0.5*(oldarray[0]+fistarray(1,3,i));
for(k=1;k<=i;k++)//给新数组第一个元素以后的赋值
{
newarray[k]=(pow(4,k)*newarray[k-1]-oldarray[k-1])/(pow(4,k)-1);
}
}
printf("%f",newarray);

}
double fistarray(double a,double b,int j)//帮助计算新数组第一个元素
{
int i;
double s=0,ss=0;
for(i=1;i<=pow(2,j-1);i++)
{
s+=jifen(a+(2*i-1)*((b-a)/pow(2,j)));
}
ss=(b-a)/pow(2,j-1)*s;
return(ss);
}