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