Интеграл, непонятен один момент Опции

[Zelepuk]

Хочу найти функцию вычисления определённого интеграла для произвольного массива чисел.
В сети полно примеров взятия интеграла, например:

CODE

double integral(a, b, f, eps)
double a, b; /* концы отрезка */
double eps; /* требуемая точность */
double (*f)(); /* подынтегральная функция */
{
register long i;
double fab = (*f)(a) + (*f)(; /* сумма на краях */
double h, h2; /* шаг и удвоенный шаг */
long n, n2; /* число точек разбиения и оно же удвоенное */
double Sodd, Seven; /* сумма значений f в нечетных и в
четных точках */
double S, Sprev;/* значение интеграла на данной
и на предыдущей итерациях */
double x; /* текущая абсцисса */

niter = 0;
n = 10L; /* четное число */
n2 = n * 2;

h = fabs(b - a) / n2; h2 = h * 2.0;

/* Вычисляем первое приближение */
/* Сумма по нечетным точкам: */
for( Sodd = 0.0, x = a+h, i = 0;
i < n;
i++, x += h2 )
Sodd += (*f)(x);

/* Сумма по четным точкам: */
for( Seven = 0.0, x = a+h2, i = 0;
i < n-1;
i++, x += h2 )
Seven += f(x);

/* Предварительное значение интеграла: */
S = h / 3.0 * (fab + 4.0 * Sodd + 2.0 * Seven );
do{
niter++;
Sprev = S;

/* Вычисляем интеграл с половинным шагом */
h2 = h; h /= 2.0;
if( h == 0.0 ) break; /* потеря значимости */
n = n2; n2 *= 2;

Seven = Seven + Sodd;
/* Вычисляем сумму по новым точкам: */
for( Sodd = 0.0, x = a+h, i = 0;
i < n;
i++, x += h2 )
Sodd += (*f)(x);

/* Значение интеграла */
S = h / 3.0 * (fab + 4.0 * Sodd + 2.0 * Seven );

} while( niter < 31 && fabs(S - Sprev) / 15.0 >= eps );
/* Используем условие Рунге для окончания итераций */

return ( 16.0 * S - Sprev ) / 15.0 ;
/* Возвращаем уточненное по Ричардсону значение */
}

Непонятно следущее: f - указывает на подынтегральную функцию. Математически это понятно. Но в Си же это подпрограмма! Я наверное принципиально ен понимаю чего-то но как можно интеграл от подпрограммы взять? или программа должна возвращать указатель на массив?
я же надеялся найти реализацию где в качестве указателя будет указатель на массив, поскольку функция представлена в виде набора отсчётов (массива) от a до b.

Заранее спасибо за разъяснения.

Сообщение отредактировал Zelepuk - Dec 14 2012, 09:38

[Xenia]

А кстати, а чего это народ функцию от указателя так вычурно вычисляет?
S += (*f)(х);
На мой взгляд (а мои взгляды могут расходиться со взглядами компилятора ), такую функцию следует вызывать как обычную
S += f(х);
Тем более что прецедент уже есть - существует два разных способа объявления массива:
int array[100], *parray=array;
И тем не менее, обоими массивами пользуются одинаково, когда достают из них элементы:
S += array[x];
или
S += parray[x];
а вовсе не (*parray)[x];
И все только потому, что array - тоже указатель, только константный. А, значит, и оператор [] будет их разыменовывать правильно в обоих случаях.

По аналогии с этим, и функцию по указателю можно вызывать напрямую, как f(x), поскольку обычная функция (не указатель) тоже является указателем на функцию, только константным. Значит, и оператор () тоже должен работать в обоих случаях одинаково.
Или я не права?

[demiurg_spb]

Или я не права?

не совсем. попробуйте двумерный массив.
Ну или посмотрите что вам вернёт sizeof(array1) и что sizeof(array2) для таких вот случаев

Код

uint8_t array1[10];
uint8_t array2[2][5];

[Xenia]

не совсем. попробуйте двумерный массив.
Ну или посмотрите что вам вернёт sizeof(array1) и что sizeof(array2) для таких вот случаев

Код

uint8_t array1[10];
uint8_t array2[2][5];

Ну так я же не просила вас вызывать двухмерную функцию.
Речь шла о том, что оператор [] приложим к указателю, так же как и оператор (). Оба они требуют, чтобы слева от них стоял указатель, хоть бы и вычисляемый. Т.е. я не за идентичность указателей ратовала, а за единообразие выполнения операторов "скобочки". А именно, в выражении f() скобочки должны интерпретировать f, как указатель (в данном случае на функцию), а не как-то иначе. Зачем же тогда f предварительно разыновывать? Ведь если имя массива (без скобочек) это указатель, то и имя фукции без скобочек тоже должно быть указателем. А чем же ещё-то?