quad3d.cpp

提供了6种数值积分方法的c++代码实现： 1 梯形求积法(TRAPZD->QTRAP)； 2 辛普森(Simpson)求积法(QSIMP)； 3 龙贝格(Romberg)求积法(QROMB)；

double f(double zz)
{
    z = zz;
    return func(x, y, z);
}

double g(double yy)
{
	double b1,b2,ss;
    y = yy;
    b1 = z1(x, y);
    b2 = z2(x, y);
    qgausz(b1, b2, ss);
    return ss;
}

double h(double xx)
{
	double a1,a2,ss;
    x = xx;
    a1 = y1(x);
    a2 = y2(x);
    qgausy(a1, a2, ss);
    return ss;
}

void quad3d(double xx1, double xx2, double& ss)
{
    qgausx(xx1, xx2, ss);
}