?? qromb.c
字號:
#include <math.h>
#define EPS 1.0e-6
#define JMAX 20
#define JMAXP JMAX+1
#define K 5
float qromb(func,a,b)
float a,b;
float (*func)();
{
float ss,dss,trapzd();
float s[JMAXP+1],h[JMAXP+1];
int j;
void polint(),nrerror();
h[1]=1.0;
for (j=1;j<=JMAX;j++) {
s[j]=trapzd(func,a,b,j);
if (j >= K) {
polint(&h[j-K],&s[j-K],K,0.0,&ss,&dss);
if (fabs(dss) < EPS*fabs(ss)) return ss;
}
s[j+1]=s[j];
h[j+1]=0.25*h[j];
}
nrerror("Too many steps in routine QROMB");
}
#undef EPS
#undef JMAX
#undef JMAXP
#undef K
?? 快捷鍵說明
復(fù)制代碼
Ctrl + C
搜索代碼
Ctrl + F
全屏模式
F11
切換主題
Ctrl + Shift + D
顯示快捷鍵
?
增大字號
Ctrl + =
減小字號
Ctrl + -