численное интегрирование, ищу код на С Опции

[TigerSHARC]

ищу исходники взятия интеграла методом прямоугольников и сплайнов. нужны исходники на С. Наверняка такое в институтах на лабах есть.

[Gate]

Рекомендую книгу "numerical recipes in с".
http://astronu.jinr.ru/wiki/upload/d/d6/Nu...lRecipesinC.pdf
Можно найти и исходники, если из текста брать неудобно.

[_pv]

ищу исходники взятия интеграла методом прямоугольников и сплайнов. нужны исходники на С.

Код

double Integral(double * x, double * y, int num){
  double r = 0;
  for (int i = 1; i < num; i++) r += y[i] * (x[i] - x[i-1]);
//  for (int i = 1; i < num-1; i++){r += 0.5L * y[i] * (x[i+1]-x[i-1]);} r += 0.5L * (y[0] * (x[1]-x[0]) + y[num-1] * (x [num-1]-x[num-2]));
//  for (int i = 1; i < num; i++) r += 0.5L * (y[i] + y[i-1]) * (x[i] - x[i-1]);
  return r;
}

есть некие различия для неравномерной сетки по х.

еще есть промежуточный вариант симпсона, где не через две точки проводится прямая, как в коде выше, а через три - парабола, ну или через 4 кубический полином. Радость по сравнению со сплайнами, что не надо по массиву несколько раз бегать, для построения честных кубических сплайнов, можно считать на ходу.

за исходниками сюда: alglib