请教有关龙贝格求微积分的问题
题目是用龙贝格(Romberg)求pow(3,x)*(5*x+7)*sin(x*x)*pow(x,1.4)在(1,3)区间上的积分这是我这个菜鸟编的程序 但是不显示任何结果 正确答案是440左右 请大家帮帮忙 在下不胜感激 谢谢。
#include "stdio.h"
#include<math.h>
#define n 20
double f(double x)
{double z;
z=pow(3,x)*(5*x+7)*sin(x*x)*pow(x,1.4);
return(z);
}
main()
{
int i,j;
double T[n][n];
double h,temp;
T[1][0]=f(1)+f(3);
for(i=1;i<n;i++)
{ h=0;
for(j=1;j<=pow(2,i-1);j++)
temp=2*f(1+2*(2*j-i)/pow(2,i))/pow(2,i-1);
h=h+temp;}
T[1][i]=(T[1][i-1]+h)/2;
for(i=1;i<n;i++)
{for(j=1;j<=n-i+1;j++)
T[i+1][j-1]=(pow(4,i)*T[i][j]-T[i][j-1])/(pow(4,i)-1);
}
printf("%f",T[i][0]);
}