Правильно ли написан код для метода трапеций на с++?

Question

[tr1ck]

h = (b - a) / n;
	for (int i = 1; i <= (n - 1); i++)
		S += f(a + h*i);
	integ = h * ((f(a) + f(b)) / 2 + S);

Answer 1

[rdev]

Хотя бы описал что у тебя означает f(). А так какой то огрызок кода, который еще надо самому додумать.

Пример с киберфорума.

#include <cmath>
#include <cstdio>
#include <cstdlib>
 
typedef double(*function)(double);
 
double integral(function f, double a, double b, unsigned step_count) {
  double sum = .0, step;
  size_t i;
  if (0 == step_count) return sum;
 
  step = (b - a) / (1.0 * step_count);
  for ( i = 1 ; i < step_count ; ++i ) {
    sum += f (a + i * step);
  }
  sum += (f(a) + f(b)) / 2;
  sum *= step; 
  return sum;
}
 
double f (double x) {
  return 2 * x;
}
 
int main() {
  printf ("\\int_0^10(x) = %f\n", integral(f, 0, 10, 15));
}

Comments

[tr1ck]

#include "stdafx.h"
#include
#include
#include

using namespace std;

double f(double x)
{
double E = 0.000000001, S = 1, R = 0, N, n = 0, y, Y = 1, F = 1;

//Разложение 1/sqrt(1-x) в ряд Маклорена и извлечение корня
for (int i = 1; i < 10; i++){
S += pow(x,i);

}

return sqrt(S)*(2 * x / sqrt ((1+x)*(1+x*x)));
}

int _tmain(int argc, _TCHAR* argv[])
{
setlocale(LC_ALL, "Russian");
cout « "Программа считает интеграл 2x/sqrt(1-x^4) с границами от 0 до 1 методом трапеции";
cout « "с точность E=10^(-9)!" « endl;
double a = 0, b = 1, n = 4, S = 0, h, integ, IntegPrev;
h = (b - a) / n;
for (int i = 1; i <= (n - 1); i++)
S += f(a + h*i);
integ = h * ((f(a) + f(b)) / 2 + S);
const double PI = 3.14159265;
//посчитали первый интеграл для n=4
do {
//IntegPrev = integ;
n += 1;
h = (b - a) / n;
S = 0;
for (int i = 1; i <= (n - 1); i++)
S += f(a + h*i);
integ = h * ((f(a) + f(b)) / 2 + S);
} while (abs(PI/2 - integ) > 0.01);
cout « "\nAnswer: " « integ « endl;
getchar();
return 0;
}

Answer 2

[Вячеслав]

Отвечая на вопрос: Да, этот кусок кода правильный.

Рекомендации на будущее:
1) Задавая вопрос желательно писать чуть подробней.
2) Для проверки кода можно составить формулу которую он реализует и сравнить с теоретической.