как численно интегрировать переменную, которая вычисляется в программе в качестве указателя (например, с использованием трапециевидного правила) на языке c

Я не знаю, откуда взялся этот код, но он очень уродлив и имеет жестко закодированные ограничения (333 точки и приращение на 0,003). Чтобы использовать его, вам нужно правильно "сэмплировать" вашу функцию и генерировать пары (x, f(x))...

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

Давайте рассмотрим вашу функцию и предположим, что она работает (я думаю, что это не так, это действительно неясный код...; например, когда вы интегрируете функцию, вы ожидаете в результате число; где это число? Может быть, INTE5? обратно не возвращается... а если это так, то зачем окончательное обновление массива W?Это бесполезно, а может у нас есть что-то значимое в W?). Как бы вы его использовали?

Прототип

  void WORK(double *d, double *W);

означает, что РАБОТА требует два указателя. Какими должны быть эти указатели, зависит от кода; взгляд на него предполагает, что вам действительно нужны два массива с элементами N_X_POINTS каждый. Код читает и записывает в массив W и читает только из d. N_X_POINTS int равно 333, поэтому вам нужно перейти к массивам функций как минимум из 333 двойников:

   double d[333];
   double W[333];

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

В любом случае, если вы вызовете его с помощью WORK(d, W), вы не получите ошибку seg, так как массивы достаточно велики. Результат будет неверным, но это сложнее отследить (опять же, извините, для этого нет «реверс-инжиниринга»).

Последнее замечание (тоже из комментариев): если у вас double a[N], то a имеет тип double *.

Ошибка сегментации часто возникает в C, когда вы пытаетесь получить доступ к какой-либо части памяти, к которой вы не должны обращаться. Я подозреваю, что виновато выражение d[N_X_POINTS] (поскольку массивы в C имеют нулевой индекс), но, не видя определения d, я не могу быть уверен.

Попробуйте помещать информативные printf операторы отладки до/после каждой строки кода в вашей функции, чтобы вы могли сузить круг возможных источников проблемы.

Вот простая программа, которая интегрирует $f(x) = x^2$ в диапазоне [0..10]. Это должно направить вас в правильном направлении.

#include <stdio.h>
#include <stdlib.h>

double int_trapezium(double f[], double dX, int n)
{
   int i;
   double sum;

   sum = (f[0] + f[n-1])/2.0;
   for (i = 1; i < n-1; i++)
       sum += f[i];
   return dX*sum;
}

#define N 1000
int main()
{
    int i;
    double x;
    double from = 0.0;
    double to = 10.0;
    double dX = (to-from)/(N-1);
    double *f = malloc(N*sizeof(*f));

    for (i=0; i<N; i++)
    {
        x = from + i*dX*(to-from);
        f[i] = x*x;
    }
    printf("%f\n", int_trapezium(f, dX, N));
    free(f);
    return 0;
}