jak tą funkcje scałkować poprawnie?

 double f(double x, double y)
 {
 return A*x*x+B*y*y-C;
 }
double c()
{double s;
for (int i=1;i<=1000;i++)
s+=f(x1+i*hx,y1+i*hy);s*=hx*hy; }

tumdo2 napisał(a)

 double f(double x, double y)
 {
 return A*x*x+B*y*y-C;
 }
double c()
{double s;
for (int i=1;i<=1000;i++)
s+=f(x1+i*hx,y1+i*hy);s*=hx*hy; }

Wg mnie wyjątkowo namieszałeś. Porządnie podaj fcję podcałkową i granice całkowania.

algorytm najprostszy prosty (naiwny).

double całkuj2d(double (*f)(double, double), 
    double xmin,
    double xmax,
    double ymin,
    double ymax,
    double mesh)
{
    assert(f);
    double sum = 0;

    for(double x=xmin; x<xmax; x+=mesh)
         for(double y=ymin; y<ymax; y+=mesh)
                sum += (*f)(x,y);
    return sum*mesh*mesh;
}