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

[heaad]

用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);
}

[星羽]

你讲讲你的算法吧

[heaad]

我已经自己解决了。谢谢各位

[redjuly]

关键位置log出去或者trace一下看看

[HelloDan]

#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);
}
//楼主还是重贴一次代码吧。用格式贴现率

[yuncmt]

你的问题　是怎么解决的　告诉下