Аппроксимация методом Ньютона

@genrix · Регистрация: 05.01.2016

Author24 — интернет-сервис помощи студентам

Доброго времени суток, интересно найдутся среди множества программистов на сайте способные решить данную задачу:

Требуется аппроксимировать функцию X^2-10*sinквадратx на итервале [0,3]. Кол-во неизвестных параметров 4, кол-во точек 4. Таблица исходной функции yi=f(xi) вычисляется в точках xi=0+(i-1)3/3 , i=1,4. Используя полученную таблицу (хi,yi) требуется вычислить значения функций f(xi), фи(xj,c) и погрешность d(xj)=f(xj)-фи(xj,c) в точках xj=0+(j-1)3/20; j=1,21. Вот мой код программы, но мне кажется не совсем верный, да и ошибки Invalid indirection выдаёт.

C

#include<stdio.h>
#include<conio.h>
#include<math.h>
float func(float x)
{
return pow(x,3)+10*pow(x,2);
}
void main()
{
int m=5,n=5,i,j,k,del,N;
double x[5], y[5],Xj[21],Yj[21], Fj[21], a=-8, b=2,p,e;
clrscr();
getch();
for(i=0;i<m;i++)
{
x[i]=a+i*(b-a)/(m-1);
y[i]=func(x[i]);
}
for(i=0;i<21;i++)
{
Xj[i]=a+i*(b-a)/(m-1);
Yj[i]=func(Xj[i]);
}
for(j=0;j<21;j++)
{
for(i=0;i<n;i++)
{
del[i]=y[i];
}
p=1;
N=y[i];
for(k=1;k<n-1;k++)
{
for(i=1;i<n-k;i++)
{
del[i]=(del[i]-del[i+1])/(x[i]-x[i+k]);
}
p=p*(Xj[j]-x[k]);
N=N+p*pow(del[i],k);
}
Fj[j]=p;
}
printf("X f(x)\n");
for(i=0;i<5;i++)
printf("%f %f\n",x[i],y[i]);
getch();
clrscr();
printf("x f(x) g(x)\n");
for(i=0;i<21;i++)
printf("%f\t %f\t %f\t Pogr: %f\n",Xj[i],Yj[i],Fj[i],fabs(Fj[i]-Yj[i]));
getch();
}

	28.12.2011, 19:31