c++运行程序时遇到问题需要关闭

heaad 2008-04-14 12:11:55
用Romberg法求上限为9,下限为1,根号x的积分
编译可以通过,运行时就出现“4_8.exe遇到问题需要关闭.....",大家帮我看一下是不是越界了,还是其他原因

#include<iostream.h>
#include<math.h>
double f(double);
int k,m;
double a=1,b=9,kexi=1e-6,q,sum=0;
double T[20][20]={0};
void main()
{
T[0][0]=(b-a)*(f(a)+f(b))/2;
k=1;
flag:
q=int(pow(2,k-1));
for(int i=1;i<=q;i++)
{
sum+=f(a+(2*i-1)*(b-a)/(q*2));
}
T[0][k]=0.5*(T[0][k-1]+(b-a)*sum/q);
for(m=1;m<=k;m++)
{
T[m][k-m]=(pow(4,m)*T[m-1][k-m+1]-T[m-1][k-m])/(pow(4,m)-1);
}
if(fabs(T[k][0]-T[k-1][0])<kexi)
cout<<"满足精度的积分近似值为"<<T[k][0]<<endl;
else
{
k++;
goto flag;
}
}
double f(double x)
{
return sqrt(x);
}

...全文
169 5 打赏 收藏 转发到动态 举报
写回复
用AI写文章
5 条回复
切换为时间正序
请发表友善的回复…
发表回复
yuncmt 2008-04-21
  • 打赏
  • 举报
回复
你的问题 是怎么解决的 告诉下
HelloDan 2008-04-14
  • 打赏
  • 举报
回复

#include <iostream.h>
#include <math.h>
double f(double);
int k,m;
double a=1,b=9,kexi=1e-6,q,sum=0;
double T[20][20]={0};
void main()
{
T[0][0]=(b-a)*(f(a)+f(b))/2;
k=1;
flag: //这个是怎样的,实在看不明白。
q=int(pow(2,k-1)); // pow()这个在参数方面也有问题。
for(int i=1;i <=q;i++)
{
sum+=f(a+(2*i-1)*(b-a)/(q*2));
}
T[0][k]=0.5*(T[0][k-1]+(b-a)*sum/q);
for(m=1;m <=k;m++)
{
T[m][k-m]=(pow(4,m)*T[m-1][k-m+1]-T[m-1][k-m])/(pow(4,m)-1);
}
if(fabs(T[k][0]-T[k-1][0]) <kexi)
cout < <"满足精度的积分近似值为" < <T[k][0] < <endl;
else
{
k++;
goto flag;
}
}
double f(double x)
{
return sqrt(x);
}
//楼主还是重贴一次代码吧。用格式贴现率

redjuly 2008-04-14
  • 打赏
  • 举报
回复
关键位置log出去或者trace一下看看
heaad 2008-04-14
  • 打赏
  • 举报
回复
我已经自己解决了。谢谢各位
星羽 2008-04-14
  • 打赏
  • 举报
回复
你讲讲你的算法吧

64,759

社区成员

发帖
与我相关
我的任务
社区描述
C++ 语言相关问题讨论,技术干货分享,前沿动态等
c++ 技术论坛(原bbs)
社区管理员
  • C++ 语言社区
  • encoderlee
  • paschen
加入社区
  • 近7日
  • 近30日
  • 至今
社区公告
  1. 请不要发布与C++技术无关的贴子
  2. 请不要发布与技术无关的招聘、广告的帖子
  3. 请尽可能的描述清楚你的问题,如果涉及到代码请尽可能的格式化一下

试试用AI创作助手写篇文章吧