简单矩形法算法问题

Question

我被要求编写一个矩形方法算法来计算函数的估计面积，但是我的讲师将我的代码发回并要求我更正它，因为结束条件错误？

double prostokaty(double( *f)(double))
{
    double krok, p, suma = 0;
    krok = (c_do - c_od) / lp;
    for (p = c_od + krok; p < c_do; p += krok) {
        suma += (*f)(p);
    }
    return suma * krok;      
}

我花了将近两个小时，找出问题所在，但我失败了。

Answer 1

最好使用矩形的中间c_od + krok/2，而不是右边缘c_od + krok。

double prostokaty(double( *f)(double)) {
  double p, suma = 0;
  double krok = (c_do - c_od) / lp;
  // for (p = c_od + krok; p < c_do; p += krok) {
  for (p = c_od + krok/2; p < c_do; p += krok) {
    suma += (*f)(p);
  }
  return suma * krok;      
}

如果lp是整数类型，最好避免累积增量错误。

  double suma = 0;
  double krok = (c_do - c_od) / lp;
  for (int i = 0; i < lp; i++) {
    double p = c_od + (i + 0.5)*krok;
    suma += f(p);
  }
  return suma * krok;