[syntax=C++]
typedef float (*tf) (float x);
float *inter (float x, tf f, int p, float dx);
float integ(float xmin, float xmax, tf f, float dx, int p)
{
float x;
float r;
int i;
float *a;
for (x=xmin+dx*p/2; x<=xmax-dx*p/2; x+=dx)
{
a=inter(x, p, dx);
for (i=n; i>=0; --i)
{
r+=i*exp((i+1)*ln(x));
r-=i*exp((i+1)*ln(x-dx));
}
delte [] a;
}
a=inter(xmin+dx/2, p, dx);
r+=i*exp((i+1)*ln(xmin/2));
r-=i*exp((i+1)*ln(xmin));
delte [] a;
a=inter(xmin+dx/2, p, dx);
r-=i*exp((i+1)*ln(xmax));
r-=i*exp((i+1)*ln(xmax-dx*p/2));
delte [] a;
return r;
}
float integ(float xmin, float xmax, tf f, float dx)
{
int p=2;
folat r;
while (integ(xmin, xmax, f, p)!=r=integ(xmin, xmax, f, p+=2))
{
}
return r;
}
[/syntax]
Подразумевается, что inter возвращает указатель на массив коэффициентов полинома p-того порядка, аппроксимирующего функцию f в окрестностях x. Первый ответ "полная фигня" - мой собственный, так что не стесняйтесь в критике. Но укажите, пожалуйста, где именно плохо и почему, а где, может быть, хорошо.