Здравствуйте, помогите исправить ошибки, пожалуйста 
Вот задачка:
Данная n + 1 пара значений (x; y), связанных функциональной зависимостью f (x) вашего варианта. Значение x выбраны в n + 1 равноудаленных точках отрезка [a; b]. Необходимо построить интерполяционный много-член Pn (x). Для этого:
1) Рассчитайте значение функции f (x) в равноудаленных точках отрезка [a; b]. Тем самым будут заполнены два массива x и y, которые можно интерпретировать как экспериментально измеренные значения каких-либо физических величин.
2) Рассчитайте разделены разницы по значениям функции в n + 1 узле.
3) Запрограммируйте вычисления интерполяционного многочлена Pn (x) по схеме Горнера.
4) Постройте в графическом окне графики исходной функции f (x) и интерполяционного многочлена Pn (x) на отрезке [a; b +]. Значение> 0 Пожалуйста, выберите таким, чтобы область построения графика по оси x была в 1: 5 2 раза больше отрезка, на котором выбрали узлы дискретизации.
5) Опыт поведение интерполяционного многочлена при различной количестве узлов дискретизации (4,5,7,10 точек).
6) Смоделируйте погрешность задания значений функции в узлах дискретамизации путем добавления к значениям функции случайных чисел.
| Код |
#include <math.h>
#include <windows.h>
#define e 2.71
int main()
{
// Инициализация графики
int gd, gm;
detectgraph(&gd,&gm);
initgraph(&gd,&gm,"");
// цвет фона - белый
setbkcolor(WHITE);
// очистка окна
cleardevice();
// текущий цвет - черный
setcolor(BLACK);
// координаты центра
int X0=getMaxx()/2;
int Y0=getMaxy()/2;
// оси координат
line(0,Y0,getMaxx(),Y0);
line(X0,0,X0,getMaxy());
// единичные отрезки в пикселях
int edx=20;
int edy=15;
// текущий цвет - зелёный
setcolor(LIGHTBLUE);
// переменная цикла
int k;
// рисуем засечки
// по оси х
for(k=0;edx*k<X0;k++)
{
line(X0+k*edx,Y0-3,X0+k*edx,Y0+3); line(X0-k*edx,Y0-3,X0-k*edx,Y0+3);}
// по оси у for(k=0;edy*k<Y0;k++)
{
line(X0-3,Y0+k*edy,X0+3,Y0+k*edy); line(X0-3,Y0-k*edy,X0+3,Y0-k*edy);}
// текущий цвет - розовый
setcolor(LIGHTRED);
// стиль линии - сплошная толстая
setlinestyle(SOLID_LINE,0,THICK_WIDTH);
// описание переменных
float a,b,x,y,delta;
// интервал дискретизации
delta=0.1;
// границы построения по оси х
a=-2; b=2;
// Вычисление начальных координат
x=a;
y=e^(-x);
// перемещение указателя в начало графика
moveto(X0+edx*x,Y0-edy*y);
// цикл построения графика
for(x=a;x<=b;x=x+delta)
{
y=e^(-x);
lineto(X0+edx*x,Y0-edy*y);};
// установка стиля текста
settextstyle(SANS_SERIF_FONT,0,2);
// вывод заголовка
outtextxy(180,10,"График функции f(x)=e^(-x)");
getch(); // ожидание
// закрытие графического окна
closegraph();
return 0;
}
|
f(x) [a,b]
e^(-x) [-2,2] |
