Найти значения функции F(x) в указанных точках x(j) и сравнить с точными значениями f(x(j))-F(x(j))

// InterpolationDlg.cpp : implementation file
//
 
#include "stdafx.h"
#include "Interpolation.h"
#include "InterpolationDlg.h"
#include "math.h"
 
#ifdef _DEBUG
#define new DEBUG_NEW
#endif
 
 
 
// CInterpolationDlg dialog
 
#define POINT(x,y) (int)(((x)+x0)*xMash),(int)(((y)+y0)*yMash)// Директива #define определяет подстановку в тексте программы
//Используется для определения макросов:
//#define имя( параметры ) текст_подстановки
// пересчитывает старые координаты в новые
 
CInterpolationDlg::CInterpolationDlg(CWnd* pParent /*=NULL*/)
    : CDialog(CInterpolationDlg::IDD, pParent)
    , t(0)
{
    m_hIcon = AfxGetApp()->LoadIcon(IDR_MAINFRAME);
}
 
void CInterpolationDlg::DoDataExchange(CDataExchange* pDX)
{
    CDialog::DoDataExchange(pDX);
    DDX_Text(pDX, IDC_EDIT1, t);
}
double f(double x)
{
    //return x*x*x*x*x;
return 1/x;
}
double fac (int z)
{   
    double fn=1.;
    for ( int i=1; i<=z; i++)
    {
        fn*=i;
    }
    double f=1/fn;
    return f;
}
 
 
double F (double x, int n) //функция подсчета приближенного значения функции в точке.
 
 
//                      построение таблицы разностей
 
{
    double a=0.1, b=1.1;
    // аргумент первого и последнего узла интерполяции
    // n -колличество узлов интерполяции
    double h=(b-a)/(2*n);
    //шаг между узлами интерполяции
    double * x_arg = new double [n+n+10+1];     //массив аргументов функции в узлах интерполяции
    //cout<<endl<<"h:   "<<h<<endl;
 
    int s=0;//тут введена х чтоб умножать h на положит
    for (int i=-n; i<=n+1; ++i) 
    {
        x_arg[i]=a+s*h; s++;// заполнение массива аргументов от -n до n
    //  cout<<x_arg[i]<<"   ";
    }
    //cout<<endl<<"___________________"<<endl;
    //cout<<endl<<endl;
 
 
    double ** razn= new double *[n+n+10];   //создание динамического массива для таблицы конечных разностей
    // [n+n+2]столько у нас строк в таблице
    int k=0;
    for (int i=-n;i<=(n+1);++i)
    {
        razn[i] = new double [n+n+10+k];
        --k;// выделяем память для каждой строки определенное кол столбцов [n+n+1+k] 
        //к- уменьшает число столбцов в строчке на 1, т.к.
        //таблица разностей имеет диагональный вид
    }
    for ( int i=-n;i<=(n+1);++i)// заполняем нулевой столбец массива значениями функции в узлах
    {
        razn [i][0]=f(x_arg[i]);
    //  cout<<razn[i][0]<<"  "<<endl;
    }
    //cout<<endl<<"----------------"<<endl;
    int r=n+1;
    for (int j=1;j<=(n+n+1);++j)
    {
        for(int i=-n;i<r;++i)
        {
            int s=-i;
            razn[i][j]=razn[-(s-1)][j-1]-razn[i][j-1];//получаем таблицу разностей
    //      cout<<razn[i][j]<<"  "<<endl;
        }
    //  cout<<endl<<"----------------"<<endl;
        --r;//это чтоб уменьшать число столбцов, в строке, приводить к диагональному виду таблицу разности
        //cout<<endl<<"----------------"<<endl;
    }
 
 
    //                      Считаем приближенное значение функции в точке
 
 
 
 
    double q,w1,w2,w;
    w=0.;
    q=(x-x_arg[0])/h; //cout<<"q :"<<q<<endl;
    w2=(q-0.5);w1=1.;
    for(int i=1;i<=n;i++)
    {
        w1*=(q+i-1)*(q-i);//cout<<"w1:  "<<w1<<endl;
        w2*=(q-i)*(q+i-1);//cout<<"w2:  "<<w2<<endl;//умножение с присваиванием 
        w+=w1*(razn[-i][2*i]+razn[-i+1][2*i])*0.5*fac(2*i)+w2*razn[-i][2*i+1]*fac(2*i+1);//сложение с присваиванием
        //cout<<endl<<"w:  "<<w<<"    "<<endl;
    }
    //приближенное значение функции в точке
    //cout<<"  fn net w:  "<<fn<<endl;
    double fn=(razn[0][0]+razn[1][0])*0.5+(q-0.5)*razn[0][1]+w;
    return fn;  
    delete[] x_arg;
    delete [] razn;
}
 
 
double pogr (double x, int n)// считает погрешность
{
    double fn=fabs(F(x,n)-f(x));
    return fn;
}
 
BEGIN_MESSAGE_MAP(CInterpolationDlg, CDialog)
    ON_WM_PAINT()
    ON_WM_QUERYDRAGICON()
    //}}AFX_MSG_MAP
    ON_BN_CLICKED(IDC_BUTTON1, &CInterpolationDlg::OnBnClickedButton1)
END_MESSAGE_MAP()
 
 
// CInterpolationDlg message handlers
 
BOOL CInterpolationDlg::OnInitDialog()//инициализация объектов, которые будут передаваться в начале запуска программы
{   
    //setlocale(LC_ALL, ".1251");
    CDialog::OnInitDialog();
 
    // Set the icon for this dialog.  The framework does this automatically
    //  when the application's main window is not a dialog
    SetIcon(m_hIcon, TRUE);         // Set big icon
    SetIcon(m_hIcon, FALSE);        // Set small icon
 
    // TODO: Add extra initialization here
 
 
//*************************************************************************************************************************
 
    t = 2;
    UpdateData(false);  //из памяти в окно, обновляет данные, это для n
//*************************************************************************************************************************
 
    // пишем во всех программах связанных с картинками, а вначале мы их объявили в h.файле
 
    picWnd=GetDlgItem(IDC_PIC);                             // передали ID картинки 
    picDC=picWnd->GetDC();                                  // получили контекст устройства 
    picWnd->GetClientRect(&pic);                            // получаем размер картинки
 
    //карандаши которыми будем рисовать
 
    penOs.CreatePen(PS_SOLID, 2, RGB(0,0,0));               // проинициализировали карандаш который рисует черные оси 
    penSetka.CreatePen(PS_SOLID,1,RGB(100,100,100));        //сетка серая
    penGrafik.CreatePen(PS_SOLID,2,RGB(255,0,0));           //график красный
    penInterpolation.CreatePen(PS_SOLID,1,RGB(0,255,0));
    penkrug.CreatePen(PS_SOLID, 5, RGB (255,0,0));
 
    xMin=-3; 
    xMax=3; 
    yMin=-2;
    yMax=2; 
    x0=-xMin;
    y0=-yMax;
    xMash=pic.Width()/(xMax-xMin); yMash=-pic.Height()/(yMax-yMin);//масштабирующий коэффициент 
    //int k=pic.Width();
 
    return TRUE;  // return TRUE  unless you set the focus to a control
}
 
// If you add a minimize button to your dialog, you will need the code below
//  to draw the icon.  For MFC applications using the document/view model,
//  this is automatically done for you by the framework.
 
void CInterpolationDlg::OnPaint()
{
    //setlocale (LC_ALL,"rus");
 
    if (IsIconic())
    {
        CPaintDC dc(this); // device context for painting
 
        SendMessage(WM_ICONERASEBKGND, reinterpret_cast<WPARAM>(dc.GetSafeHdc()), 0);
 
        // Center icon in client rectangle
        int cxIcon = GetSystemMetrics(SM_CXICON);
        int cyIcon = GetSystemMetrics(SM_CYICON);
        CRect rect;
        GetClientRect(&rect);
        int x = (rect.Width() - cxIcon + 1) / 2;
        int y = (rect.Height() - cyIcon + 1) / 2;
 
        // Draw the icon
        dc.DrawIcon(x, y, m_hIcon);
    }
    else
    {
        CDialog::OnPaint();
        picDC->FillSolidRect(&pic,RGB(255, 255, 255));//нарисовать закрашенный прямоугольник
        picDC->SelectObject(&penOs);//контекстному устройству передаём ссылку  
        //picDC->MoveTo((int)((5+x0)*xMash),(int)((5+y0)*yMash));
        
        //picDC->LineTo((int)((7+x0)*xMash),(int)((9+y0)*yMash));
        //координатные оси
        //ось ох
        picDC->MoveTo(POINT(0,yMax));
        picDC->LineTo(POINT(0,yMin));
        // ось оу
        picDC->MoveTo(POINT(xMax,0));
        picDC->LineTo(POINT(xMin,0));
        //picDC->MoveTo(POINT(0,0));
        //picDC->LineTo(POINT(7,0));
 
        // Рисуем сетку
 
        picDC->SelectObject(&penSetka);//контекстному устройству передаём ссылку на penSetka
        for (double x=xMin; x<xMax; x+=1.)
        {
            picDC->MoveTo(POINT(x,yMin));
            picDC->LineTo(POINT(x,yMax));
        }
        for (double y=yMin; y<yMax; y+=1.)
        {
            picDC->MoveTo(POINT(xMin,y));
            picDC->LineTo(POINT(xMax,y));
        }
 
        // рисуем окружности, узлы интегрирования
        char *str = new char [10];//это для перевода цифр в текст
        picDC->SelectObject(&penkrug);
        for (double x=xMin; x<xMax; x+=1.)
        {
            picDC->Ellipse(POINT (x,0),POINT((x+0.01),(0+0.01)));
            picDC->TextOutW(POINT (x,0),(CString)itoa(x,str,10));//передаём нашу координату в строку, макс размер 10.
        }
        for (double y=yMin; y<yMax; y+=1.)
        {
            picDC->Ellipse(POINT (0,y),POINT((0+0.01),(y+0.01)));
            picDC->TextOutW(POINT (0,y),(CString)itoa(y,str,10));//itoa переводит числа в строку
        }
        picDC->SelectObject(&penGrafik);//контекстному устройству передаём ссылку на penGrafik
        picDC->MoveTo(POINT(xMin, f(xMin)));
        for(double x=xMin; x<xMax; x+=0.001 )
        { 
            picDC->LineTo(POINT(x,f(x)));
        }
        
 
    }
}
 
// The system calls this function to obtain the cursor to display while the user drags
//  the minimized window.
HCURSOR CInterpolationDlg::OnQueryDragIcon()
{
    return static_cast<HCURSOR>(m_hIcon);
}
 
 
void CInterpolationDlg::OnBnClickedButton1()
{
    // TODO: Add your control notification handler code here
    // Рисуем график интерполяции
//**********************************************************************************************
 
    UpdateData(true);
    
//  *************************************************************************************************
    OnPaint();
    int n=t;
    picDC->SelectObject(&penInterpolation);
    picDC->MoveTo(POINT(0, F(0,n)));
    for(double x=0.1; x<=1.1;x+=0.04)
    { 
        picDC->LineTo(POINT(x,F(x,n)));
    }
    FILE *File;
    File=fopen("pogr.txt", "w");
    char str[]="Результаты вычислений \n";
    fprintf (File,"%s%s%s%s","x   \t\t", "f(x) \t \t", "F(x) \t \t"," pogr \t \t\n");
    for(double x=0.1; x<=1.1;x+=0.25)
    {
        fprintf(File,"%f\t%f\t%f\t%f\n",x,f(x),F(x,n),pogr(x,n));
    }
    fclose(File);
}

double ** razn= new double *[n+n+10];
for (int i=-n;i<=(n+1);++i)
{
   razn[i] = new double [n+n+10+k];
}
...
delete [] razn;

for (int i=-n;i<=(n+1);++i)
{
delete razn[i]
}
delete [] razn

 
// InterpolationDlg.cpp : implementation file
//
 
#include "stdafx.h"
#include "Interpolation.h"
#include "InterpolationDlg.h"
#include "math.h"
 
#ifdef _DEBUG
#define new DEBUG_NEW
#endif
 
 
 
// CInterpolationDlg dialog
 
#define POINT(x,y) (int)(((x)+x0)*xMash),(int)(((y)+y0)*yMash)// Директива #define определяет подстановку в тексте программы
//Используется для определения макросов:
//#define имя( параметры ) текст_подстановки
// пересчитывает старые координаты в новые
 
CInterpolationDlg::CInterpolationDlg(CWnd* pParent /*=NULL*/)
        : CDialog(CInterpolationDlg::IDD, pParent)
        , t(0)
{
        m_hIcon = AfxGetApp()->LoadIcon(IDR_MAINFRAME);
}
 
void CInterpolationDlg::DoDataExchange(CDataExchange* pDX)
{
        CDialog::DoDataExchange(pDX);
        DDX_Text(pDX, IDC_EDIT1, t);
}
double f(double x)
{
        //return x*x*x*x*x;
return 1/x;
}
double fac (int z)
{   
        double fn=1.;
        for ( int i=1; i<=z; i++)
        {
                fn*=i;
        }
        double f=1/fn;
        return f;
}
 
 
double F (double x, int n) //функция подсчета приближенного значения функции в точке.
 
 
//                      построение таблицы разностей
 
{
        double a=0.1, b=1.1;
        // аргумент первого и последнего узла интерполяции
        // n -колличество узлов интерполяции
        double h=(b-a)/(2*n);
        //шаг между узлами интерполяции
        double * x_arg = new double [2*n+2];         //массив аргументов функции в узлах интерполяции
        //cout<<endl<<"h:   "<<h<<endl;
 
        int s=0;//тут введена х чтоб умножать h на положит
        for (int i=0; i<=n+1; ++i) 
        {
                x_arg[i]=a+s*h; s++;// заполнение массива аргументов от -n до n
        //      cout<<x_arg[i]<<"   ";
        }
        //cout<<endl<<"___________________"<<endl;
        //cout<<endl<<endl;
 
 
        double ** razn= new double *[2*n+2];   //создание динамического массива для таблицы конечных разностей
        // [n+n+2]столько у нас строк в таблице
        int k=0;
        for (int i=0;i<=(n+1);++i)
        {
                razn[i] = new double [2*n+2+k];
                --k;// выделяем память для каждой строки определенное кол столбцов [n+n+1+k] 
                //к- уменьшает число столбцов в строчке на 1, т.к.
                //таблица разностей имеет диагональный вид
        }
        for ( int i=0;i<=(n+1);++i)// заполняем нулевой столбец массива значениями функции в узлах
        {
                razn [i][0]=f(x_arg[i]);
        //      cout<<razn[i][0]<<"  "<<endl;
        }
        //cout<<endl<<"----------------"<<endl;
        int r=n+1;
        for (int j=1;j<=(n+n+1);++j)
        {
                for(int i=-n;i<r;++i)
                {
                        int s=-i;
                        razn[i][j]=razn[-(s-1)][j-1]-razn[i][j-1];//получаем таблицу разностей
        //              cout<<razn[i][j]<<"  "<<endl;
                }
        //      cout<<endl<<"----------------"<<endl;
                --r;//это чтоб уменьшать число столбцов, в строке, приводить к диагональному виду таблицу разности
                //cout<<endl<<"----------------"<<endl;
        }
 
 
        //                      Считаем приближенное значение функции в точке
 
 
 
 
        double q,w1,w2,w;
        w=0.;
        q=(x-x_arg[0])/h; //cout<<"q :"<<q<<endl;
        w2=(q-0.5);w1=1.;
        for(int i=1;i<=n;i++)
        {
                w1*=(q+i-1)*(q-i);//cout<<"w1:  "<<w1<<endl;
                w2*=(q-i)*(q+i-1);//cout<<"w2:  "<<w2<<endl;//умножение с присваиванием 
                w+=w1*(razn[-i][2*i]+razn[-i+1][2*i])*0.5*fac(2*i)+w2*razn[-i][2*i+1]*fac(2*i+1);//сложение с присваиванием
                //cout<<endl<<"w:  "<<w<<"    "<<endl;
        }
        //приближенное значение функции в точке
        //cout<<"  fn net w:  "<<fn<<endl;
        double fn=(razn[0][0]+razn[1][0])*0.5+(q-0.5)*razn[0][1]+w;
        return fn;      
        delete[] x_arg;
        delete [] razn;
}
 
 
double pogr (double x, int n)// считает погрешность
{
        double fn=fabs(F(x,n)-f(x));
        return fn;
}
 
BEGIN_MESSAGE_MAP(CInterpolationDlg, CDialog)
        ON_WM_PAINT()
        ON_WM_QUERYDRAGICON()
        //}}AFX_MSG_MAP
        ON_BN_CLICKED(IDC_BUTTON1, &CInterpolationDlg::OnBnClickedButton1)
END_MESSAGE_MAP()
 
 
// CInterpolationDlg message handlers
 
BOOL CInterpolationDlg::OnInitDialog()//инициализация объектов, которые будут передаваться в начале запуска программы
{       
        //setlocale(LC_ALL, ".1251");
        CDialog::OnInitDialog();
 
        // Set the icon for this dialog.  The framework does this automatically
        //  when the application's main window is not a dialog
        SetIcon(m_hIcon, TRUE);                 // Set big icon
        SetIcon(m_hIcon, FALSE);                // Set small icon
 
        // TODO: Add extra initialization here
 
 
//*************************************************************************************************************************
 
        t = 2;
        UpdateData(false);      //из памяти в окно, обновляет данные, это для n
//*************************************************************************************************************************
 
        // пишем во всех программах связанных с картинками, а вначале мы их объявили в h.файле
 
        picWnd=GetDlgItem(IDC_PIC);                                                             // передали ID картинки 
        picDC=picWnd->GetDC();                                                                  // получили контекст устройства 
        picWnd->GetClientRect(&pic);                                                    // получаем размер картинки
 
        //карандаши которыми будем рисовать
 
        penOs.CreatePen(PS_SOLID, 2, RGB(0,0,0));                               // проинициализировали карандаш который рисует черные оси 
        penSetka.CreatePen(PS_SOLID,1,RGB(100,100,100));                //сетка серая
        penGrafik.CreatePen(PS_SOLID,2,RGB(255,0,0));                   //график красный
        penInterpolation.CreatePen(PS_SOLID,1,RGB(0,255,0));
        penkrug.CreatePen(PS_SOLID, 5, RGB (255,0,0));
 
        xMin=-3; 
        xMax=3; 
        yMin=-2;
        yMax=2; 
        x0=-xMin;
        y0=-yMax;
        xMash=pic.Width()/(xMax-xMin); yMash=-pic.Height()/(yMax-yMin);//масштабирующий коэффициент 
        //int k=pic.Width();
 
        return TRUE;  // return TRUE  unless you set the focus to a control
}
 
// If you add a minimize button to your dialog, you will need the code below
//  to draw the icon.  For MFC applications using the document/view model,
//  this is automatically done for you by the framework.
 
void CInterpolationDlg::OnPaint()
{
        //setlocale (LC_ALL,"rus");
 
        if (IsIconic())
        {
                CPaintDC dc(this); // device context for painting
 
                SendMessage(WM_ICONERASEBKGND, reinterpret_cast<WPARAM>(dc.GetSafeHdc()), 0);
 
                // Center icon in client rectangle
                int cxIcon = GetSystemMetrics(SM_CXICON);
                int cyIcon = GetSystemMetrics(SM_CYICON);
                CRect rect;
                GetClientRect(&rect);
                int x = (rect.Width() - cxIcon + 1) / 2;
                int y = (rect.Height() - cyIcon + 1) / 2;
 
                // Draw the icon
                dc.DrawIcon(x, y, m_hIcon);
        }
        else
        {
                CDialog::OnPaint();
                picDC->FillSolidRect(&pic,RGB(255, 255, 255));//нарисовать закрашенный прямоугольник
                picDC->SelectObject(&penOs);//контекстному устройству передаём ссылку  
                //picDC->MoveTo((int)((5+x0)*xMash),(int)((5+y0)*yMash));
                
                //picDC->LineTo((int)((7+x0)*xMash),(int)((9+y0)*yMash));
                //координатные оси
                //ось ох
                picDC->MoveTo(POINT(0,yMax));
                picDC->LineTo(POINT(0,yMin));
                // ось оу
                picDC->MoveTo(POINT(xMax,0));
                picDC->LineTo(POINT(xMin,0));
                //picDC->MoveTo(POINT(0,0));
                //picDC->LineTo(POINT(7,0));
 
                // Рисуем сетку
 
                picDC->SelectObject(&penSetka);//контекстному устройству передаём ссылку на penSetka
                for (double x=xMin; x<xMax; x+=1.)
                {
                        picDC->MoveTo(POINT(x,yMin));
                        picDC->LineTo(POINT(x,yMax));
                }
                for (double y=yMin; y<yMax; y+=1.)
                {
                        picDC->MoveTo(POINT(xMin,y));
                        picDC->LineTo(POINT(xMax,y));
                }
 
                // рисуем окружности, узлы интегрирования
                char *str = new char [10];//это для перевода цифр в текст
                picDC->SelectObject(&penkrug);
                for (double x=xMin; x<xMax; x+=1.)
                {
                        picDC->Ellipse(POINT (x,0),POINT((x+0.01),(0+0.01)));
                        picDC->TextOutW(POINT (x,0),(CString)itoa(x,str,10));//передаём нашу координату в строку, макс размер 10.
                }
                for (double y=yMin; y<yMax; y+=1.)
                {
                        picDC->Ellipse(POINT (0,y),POINT((0+0.01),(y+0.01)));
                        picDC->TextOutW(POINT (0,y),(CString)itoa(y,str,10));//itoa переводит числа в строку
                }
                picDC->SelectObject(&penGrafik);//контекстному устройству передаём ссылку на penGrafik
                picDC->MoveTo(POINT(xMin, f(xMin)));
                for(double x=xMin; x<xMax; x+=0.001 )
                { 
                        picDC->LineTo(POINT(x,f(x)));
                }
                
 
        }
}
 
// The system calls this function to obtain the cursor to display while the user drags
//  the minimized window.
HCURSOR CInterpolationDlg::OnQueryDragIcon()
{
        return static_cast<HCURSOR>(m_hIcon);
}
 
 
void CInterpolationDlg::OnBnClickedButton1()
{
        // TODO: Add your control notification handler code here
        // Рисуем график интерполяции
//**********************************************************************************************
 
        UpdateData(true);
        
//      *************************************************************************************************
        OnPaint();
        int n=t;
        picDC->SelectObject(&penInterpolation);
        picDC->MoveTo(POINT(0, F(0,n)));
        for(double x=0.1; x<=1.1;x+=0.04)
        { 
                picDC->LineTo(POINT(x,F(x,n)));
        }
        FILE *File;
        File=fopen("pogr.txt", "w");
        char str[]="Результаты вычислений \n";
        fprintf (File,"%s%s%s%s","x   \t\t", "f(x) \t \t", "F(x) \t \t"," pogr \t \t\n");
        for(double x=0.1; x<=1.1;x+=0.25)
        {
                fprintf(File,"%f\t%f\t%f\t%f\n",x,f(x),F(x,n),pogr(x,n));
        }
        fclose(File);
}

// InterpolationDlg.cpp : implementation file
//
 
#include "stdafx.h"
#include "Interpolation.h"
#include "InterpolationDlg.h"
#include "math.h"
 
#ifdef _DEBUG
#define new DEBUG_NEW
#endif
 
 
 
// CInterpolationDlg dialog
 
#define POINT(x,y) (int)(((x)+x0)*xMash),(int)(((y)+y0)*yMash)// Директива #define определяет подстановку в тексте программы
//Используется для определения макросов:
//#define имя( параметры ) текст_подстановки
// пересчитывает старые координаты в новые
 
CInterpolationDlg::CInterpolationDlg(CWnd* pParent /*=NULL*/)
        : CDialog(CInterpolationDlg::IDD, pParent)
        , t(0)
{
        m_hIcon = AfxGetApp()->LoadIcon(IDR_MAINFRAME);
}
 
void CInterpolationDlg::DoDataExchange(CDataExchange* pDX)
{
        CDialog::DoDataExchange(pDX);
        DDX_Text(pDX, IDC_EDIT1, t);
}
double f(double x)
{
        //return x*x*x*x*x;
return 1/x;
}
double fac (int z)
{   
        double fn=1.;
        for ( int i=1; i<=z; i++)
        {
                fn*=i;
        }
        double f=1/fn;
        return f;
}
 
 
double F (double x, int n) //функция подсчета приближенного значения функции в точке.
 
 
//                      построение таблицы разностей
 
{
        double a=0.1, b=1.1;
        // аргумент первого и последнего узла интерполяции
        // n -колличество узлов интерполяции
        double h=(b-a)/(2*n);
        //шаг между узлами интерполяции
        double * x_arg = new double [2*n+2];         //массив аргументов функции в узлах интерполяции
        //cout<<endl<<"h:   "<<h<<endl;
 
        int s=0;//тут введена х чтоб умножать h на положит
        for (int i=-n; i<=n+1; ++i) 
        {
                x_arg[i+n]=a+s*h; s++;// заполнение массива аргументов от -n до n
        //      cout<<x_arg[i]<<"   ";
        }
        //cout<<endl<<"___________________"<<endl;
        //cout<<endl<<endl;
 
 
        double ** razn= new double *[2*n+2];   //создание динамического массива для таблицы конечных разностей
        // [n+n+2]столько у нас строк в таблице
        int k=0;
        for (int i=-n;i<=(n+1);++i)
        {
                razn[i+n] = new double [2*n+2+k];
                --k;// выделяем память для каждой строки определенное кол столбцов [n+n+1+k] 
                //к- уменьшает число столбцов в строчке на 1, т.к.
                //таблица разностей имеет диагональный вид
        }
        for ( int i=-n;i<=(n+1);++i)// заполняем нулевой столбец массива значениями функции в узлах
        {
                razn [i+n][0]=f(x_arg[i]);
        //      cout<<razn[i][0]<<"  "<<endl;
        }
        //cout<<endl<<"----------------"<<endl;
        int r=n+1;
        for (int j=1;j<=(n+n+1);++j)
        {
                for(int i=-n;i<r;++i)
                {
                        int s=-i;
                        razn[i+n][j]=razn[-(s-1)][j-1]-razn[i+n][j-1];//получаем таблицу разностей
        //              cout<<razn[i][j]<<"  "<<endl;
                }
        //      cout<<endl<<"----------------"<<endl;
                --r;//это чтоб уменьшать число столбцов, в строке, приводить к диагональному виду таблицу разности
                //cout<<endl<<"----------------"<<endl;
        }
 
 
        //                      Считаем приближенное значение функции в точке
 
 
 
 
        double q,w1,w2,w;
        w=0.;
        q=(x-x_arg[0])/h; //cout<<"q :"<<q<<endl;
        w2=(q-0.5);w1=1.;
        for(int i=1;i<=n;i++)
        {
                w1*=(q+i-1)*(q-i);//cout<<"w1:  "<<w1<<endl;
                w2*=(q-i)*(q+i-1);//cout<<"w2:  "<<w2<<endl;//умножение с присваиванием 
                w+=w1*(razn[-i][2*i]+razn[-i+1][2*i])*0.5*fac(2*i)+w2*razn[-i][2*i+1]*fac(2*i+1);//сложение с присваиванием
                //cout<<endl<<"w:  "<<w<<"    "<<endl;
        }
        //приближенное значение функции в точке
        //cout<<"  fn net w:  "<<fn<<endl;
        double fn=(razn[0][0]+razn[1][0])*0.5+(q-0.5)*razn[0][1]+w;
        return fn;      
        delete[] x_arg;
        delete [] razn;
}
 
 
double pogr (double x, int n)// считает погрешность
{
        double fn=fabs(F(x,n)-f(x));
        return fn;
}
 
BEGIN_MESSAGE_MAP(CInterpolationDlg, CDialog)
        ON_WM_PAINT()
        ON_WM_QUERYDRAGICON()
        //}}AFX_MSG_MAP
        ON_BN_CLICKED(IDC_BUTTON1, &CInterpolationDlg::OnBnClickedButton1)
END_MESSAGE_MAP()
 
 
// CInterpolationDlg message handlers
 
BOOL CInterpolationDlg::OnInitDialog()//инициализация объектов, которые будут передаваться в начале запуска программы
{       
        //setlocale(LC_ALL, ".1251");
        CDialog::OnInitDialog();
 
        // Set the icon for this dialog.  The framework does this automatically
        //  when the application's main window is not a dialog
        SetIcon(m_hIcon, TRUE);                 // Set big icon
        SetIcon(m_hIcon, FALSE);                // Set small icon
 
        // TODO: Add extra initialization here
 
 
//*************************************************************************************************************************
 
        t = 2;
        UpdateData(false);      //из памяти в окно, обновляет данные, это для n
//*************************************************************************************************************************
 
        // пишем во всех программах связанных с картинками, а вначале мы их объявили в h.файле
 
        picWnd=GetDlgItem(IDC_PIC);                                                             // передали ID картинки 
        picDC=picWnd->GetDC();                                                                  // получили контекст устройства 
        picWnd->GetClientRect(&pic);                                                    // получаем размер картинки
 
        //карандаши которыми будем рисовать
 
        penOs.CreatePen(PS_SOLID, 2, RGB(0,0,0));                               // проинициализировали карандаш который рисует черные оси 
        penSetka.CreatePen(PS_SOLID,1,RGB(100,100,100));                //сетка серая
        penGrafik.CreatePen(PS_SOLID,2,RGB(255,0,0));                   //график красный
        penInterpolation.CreatePen(PS_SOLID,1,RGB(0,255,0));
        penkrug.CreatePen(PS_SOLID, 5, RGB (255,0,0));
 
        xMin=-3; 
        xMax=3; 
        yMin=-2;
        yMax=2; 
        x0=-xMin;
        y0=-yMax;
        xMash=pic.Width()/(xMax-xMin); yMash=-pic.Height()/(yMax-yMin);//масштабирующий коэффициент 
        //int k=pic.Width();
 
        return TRUE;  // return TRUE  unless you set the focus to a control
}
 
// If you add a minimize button to your dialog, you will need the code below
//  to draw the icon.  For MFC applications using the document/view model,
//  this is automatically done for you by the framework.
 
void CInterpolationDlg::OnPaint()
{
        //setlocale (LC_ALL,"rus");
 
        if (IsIconic())
        {
                CPaintDC dc(this); // device context for painting
 
                SendMessage(WM_ICONERASEBKGND, reinterpret_cast<WPARAM>(dc.GetSafeHdc()), 0);
 
                // Center icon in client rectangle
                int cxIcon = GetSystemMetrics(SM_CXICON);
                int cyIcon = GetSystemMetrics(SM_CYICON);
                CRect rect;
                GetClientRect(&rect);
                int x = (rect.Width() - cxIcon + 1) / 2;
                int y = (rect.Height() - cyIcon + 1) / 2;
 
                // Draw the icon
                dc.DrawIcon(x, y, m_hIcon);
        }
        else
        {
                CDialog::OnPaint();
                picDC->FillSolidRect(&pic,RGB(255, 255, 255));//нарисовать закрашенный прямоугольник
                picDC->SelectObject(&penOs);//контекстному устройству передаём ссылку  
                //picDC->MoveTo((int)((5+x0)*xMash),(int)((5+y0)*yMash));
                
                //picDC->LineTo((int)((7+x0)*xMash),(int)((9+y0)*yMash));
                //координатные оси
                //ось ох
                picDC->MoveTo(POINT(0,yMax));
                picDC->LineTo(POINT(0,yMin));
                // ось оу
                picDC->MoveTo(POINT(xMax,0));
                picDC->LineTo(POINT(xMin,0));
                //picDC->MoveTo(POINT(0,0));
                //picDC->LineTo(POINT(7,0));
 
                // Рисуем сетку
 
                picDC->SelectObject(&penSetka);//контекстному устройству передаём ссылку на penSetka
                for (double x=xMin; x<xMax; x+=1.)
                {
                        picDC->MoveTo(POINT(x,yMin));
                        picDC->LineTo(POINT(x,yMax));
                }
                for (double y=yMin; y<yMax; y+=1.)
                {
                        picDC->MoveTo(POINT(xMin,y));
                        picDC->LineTo(POINT(xMax,y));
                }
 
                // рисуем окружности, узлы интегрирования
                char *str = new char [10];//это для перевода цифр в текст
                picDC->SelectObject(&penkrug);
                for (double x=xMin; x<xMax; x+=1.)
                {
                        picDC->Ellipse(POINT (x,0),POINT((x+0.01),(0+0.01)));
                        picDC->TextOutW(POINT (x,0),(CString)itoa(x,str,10));//передаём нашу координату в строку, макс размер 10.
                }
                for (double y=yMin; y<yMax; y+=1.)
                {
                        picDC->Ellipse(POINT (0,y),POINT((0+0.01),(y+0.01)));
                        picDC->TextOutW(POINT (0,y),(CString)itoa(y,str,10));//itoa переводит числа в строку
                }
                picDC->SelectObject(&penGrafik);//контекстному устройству передаём ссылку на penGrafik
                picDC->MoveTo(POINT(xMin, f(xMin)));
                for(double x=xMin; x<xMax; x+=0.001 )
                { 
                        picDC->LineTo(POINT(x,f(x)));
                }
                
 
        }
}
 
// The system calls this function to obtain the cursor to display while the user drags
//  the minimized window.
HCURSOR CInterpolationDlg::OnQueryDragIcon()
{
        return static_cast<HCURSOR>(m_hIcon);
}
 
 
void CInterpolationDlg::OnBnClickedButton1()
{
        // TODO: Add your control notification handler code here
        // Рисуем график интерполяции
//**********************************************************************************************
 
        UpdateData(true);
        
//      *************************************************************************************************
        OnPaint();
        int n=t;
        picDC->SelectObject(&penInterpolation);
        picDC->MoveTo(POINT(0, F(0,n)));
        for(double x=0.1; x<=1.1;x+=0.04)
        { 
                picDC->LineTo(POINT(x,F(x,n)));
        }
        FILE *File;
        File=fopen("pogr.txt", "w");
        char str[]="Результаты вычислений \n";
        fprintf (File,"%s%s%s%s","x   \t\t", "f(x) \t \t", "F(x) \t \t"," pogr \t \t\n");
        for(double x=0.1; x<=1.1;x+=0.25)
        {
                fprintf(File,"%f\t%f\t%f\t%f\n",x,f(x),F(x,n),pogr(x,n));
        }
        fclose(File);
}

int arr [n+n+2];
arr[i+n]; //Это ваш -2-ой элемент, просто n плюсуйте везде

// InterpolationDlg.cpp : implementation file
//
 
#include "stdafx.h"
#include "Interpolation.h"
#include "InterpolationDlg.h"
#include "math.h"
 
#ifdef _DEBUG
#define new DEBUG_NEW
#endif
 
 
 
// CInterpolationDlg dialog
 
#define POINT(x,y) (int)(((x)+x0)*xMash),(int)(((y)+y0)*yMash)// Директива #define определяет подстановку в тексте программы
//Используется для определения макросов:
//#define имя( параметры ) текст_подстановки
// пересчитывает старые координаты в новые
 
CInterpolationDlg::CInterpolationDlg(CWnd* pParent /*=NULL*/)
    : CDialog(CInterpolationDlg::IDD, pParent)
    , t(0)
{
    m_hIcon = AfxGetApp()->LoadIcon(IDR_MAINFRAME);
}
 
void CInterpolationDlg::DoDataExchange(CDataExchange* pDX)
{
    CDialog::DoDataExchange(pDX);
    DDX_Text(pDX, IDC_EDIT1, t);
}
double f(double x)
{
//  return sin(x*x*x);
    //return x*x*x*x*x;
    //return sin(x);
    //return sin(x*x)
return 1/x;
}
double fac (int z)
{   
    double fn=1.;
    for ( int i=1; i<=z; i++)
    {
        fn*=i;
    }
    double f=1/fn;
    return f;
}
 
 
double F (double x, int n) //функция подсчета приближенного значения функции в точке.
 
 
//                      построение таблицы разностей
 
{
    
    double a=0.1, b=1.1;
    // аргумент первого и последнего узла интерполяции
    // n -колличество узлов интерполяции
    double h=(b-a)/(2*n);
    //шаг между узлами интерполяции
    double * x_arg = new double [n+n+2+1];      //массив аргументов функции в узлах интерполяции
    //cout<<endl<<"h:   "<<h<<endl;
 
    int s=0;//тут введена х чтоб умножать h на положит
    for (int i=-n; i<=n+1; ++i) 
    {
        x_arg[i+n]=a+s*h; s++;// заполнение массива аргументов от -n до n
    //  cout<<x_arg[i]<<"   ";
    }
    //cout<<endl<<"___________________"<<endl;
    //cout<<endl<<endl;
 
 
    double ** razn= new double *[n+n+3];    //создание динамического массива для таблицы конечных разностей
    // [n+n+2]столько у нас строк в таблице
    int k=0;
    for (int i=-n;i<=(n+1);++i)
    {
        razn[i+n] = new double [n+n+3+k];
        --k;// выделяем память для каждой строки определенное кол столбцов [n+n+1+k] 
        //к- уменьшает число столбцов в строчке на 1, т.к.
        //таблица разностей имеет диагональный вид
    }
    for ( int i=-n;i<=(n+1);++i)// заполняем нулевой столбец массива значениями функции в узлах
    {
        razn [i+n][0]=f(x_arg[i]);
    //  cout<<razn[i][0]<<"  "<<endl;
    }
    //cout<<endl<<"----------------"<<endl;
    int r=n+1;
    for (int j=1;j<=(n+n+1);++j)
    {
        for(int i=-n;i<r;++i)
        {
            int s=-i;
            razn[i+n][j]=razn[-(s-1)][j-1]-razn[i][j-1];//получаем таблицу разностей
    //      cout<<razn[i][j]<<"  "<<endl;
        }
    //  cout<<endl<<"----------------"<<endl;
        --r;//это чтоб уменьшать число столбцов, в строке, приводить к диагональному виду таблицу разности
        //cout<<endl<<"----------------"<<endl;
    }
 
 
    //                      Считаем приближенное значение функции в точке
 
 
 
    double q,w1,w2,w;
    w=0.;
    //h=x_arg[0]-x_arg[-1];
    q=(x-x_arg[0])/h; //cout<<"q :"<<q<<endl;
    w2=(q-0.5);w1=1.;
    for(int i=1;i<=n;i++)
    {
        w1*=(q+i-1)*(q-i);//cout<<"w1:  "<<w1<<endl;
        w2*=(q-i)*(q+i-1);//cout<<"w2:  "<<w2<<endl;//умножение с присваиванием 
        w+=w1*(razn[-i][2*i]+razn[-i+1][2*i])*0.5*fac(2*i)+w2*razn[-i][2*i+1]*fac(2*i+1);//сложение с присваиванием
        //cout<<endl<<"w:  "<<w<<"    "<<endl;
    }
    //приближенное значение функции в точке
    //cout<<"  fn net w:  "<<fn<<endl;
    double fn=(razn[0][0]+razn[1][0])*0.5+(q-0.5)*razn[0][1]+w;
    delete[] x_arg;
    
    return fn;  
    delete [] razn;
}
 
 
double pogr (double x, int n)// считает погрешность
{
    double fn=fabs(F(x,n)-f(x));
    return fn;
}
 
BEGIN_MESSAGE_MAP(CInterpolationDlg, CDialog)
    ON_WM_PAINT()
    ON_WM_QUERYDRAGICON()
    //}}AFX_MSG_MAP
    ON_BN_CLICKED(IDC_BUTTON1, &CInterpolationDlg::OnBnClickedButton1)
END_MESSAGE_MAP()
 
 
// CInterpolationDlg message handlers
 
BOOL CInterpolationDlg::OnInitDialog()//инициализация объектов, которые будут передаваться в начале запуска программы
{   
    //setlocale(LC_ALL, ".1251");
    CDialog::OnInitDialog();
 
    // Set the icon for this dialog.  The framework does this automatically
    //  when the application's main window is not a dialog
    SetIcon(m_hIcon, TRUE);         // Set big icon
    SetIcon(m_hIcon, FALSE);        // Set small icon
 
    // TODO: Add extra initialization here
 
 
//*************************************************************************************************************************
 
    t = 2;
    UpdateData(false);  //из памяти в окно, обновляет данные, это для n
//*************************************************************************************************************************
 
    // пишем во всех программах связанных с картинками, а вначале мы их объявили в h.файле
 
    picWnd=GetDlgItem(IDC_PIC);                             // передали ID картинки 
    picDC=picWnd->GetDC();                                  // получили контекст устройства 
    picWnd->GetClientRect(&pic);                            // получаем размер картинки
 
    //карандаши которыми будем рисовать
 
    penOs.CreatePen(PS_SOLID, 2, RGB(0,0,0));               // проинициализировали карандаш который рисует черные оси 
    penSetka.CreatePen(PS_SOLID,1,RGB(100,100,100));        //сетка серая
    penGrafik.CreatePen(PS_SOLID,2,RGB(255,0,0));           //график красный
    penInterpolation.CreatePen(PS_SOLID,1,RGB(0,255,0));
    penkrug.CreatePen(PS_SOLID, 5, RGB (255,0,0));
 
    xMin=-3; 
    xMax=3; 
    yMin=-2;
    yMax=2; 
    x0=-xMin;
    y0=-yMax;
    xMash=pic.Width()/(xMax-xMin); yMash=-pic.Height()/(yMax-yMin);//масштабирующий коэффициент 
    //int k=pic.Width();
 
    return TRUE;  // return TRUE  unless you set the focus to a control
}
 
// If you add a minimize button to your dialog, you will need the code below
//  to draw the icon.  For MFC applications using the document/view model,
//  this is automatically done for you by the framework.
 
void CInterpolationDlg::OnPaint()
{
    //setlocale (LC_ALL,"rus");
 
    if (IsIconic())
    {
        CPaintDC dc(this); // device context for painting
 
        SendMessage(WM_ICONERASEBKGND, reinterpret_cast<WPARAM>(dc.GetSafeHdc()), 0);
 
        // Center icon in client rectangle
        int cxIcon = GetSystemMetrics(SM_CXICON);
        int cyIcon = GetSystemMetrics(SM_CYICON);
        CRect rect;
        GetClientRect(&rect);
        int x = (rect.Width() - cxIcon + 1) / 2;
        int y = (rect.Height() - cyIcon + 1) / 2;
 
        // Draw the icon
        dc.DrawIcon(x, y, m_hIcon);
    }
    else
    {
        CDialog::OnPaint();
        picDC->FillSolidRect(&pic,RGB(255, 255, 255));//нарисовать закрашенный прямоугольник
        picDC->SelectObject(&penOs);//контекстному устройству передаём ссылку  
        //picDC->MoveTo((int)((5+x0)*xMash),(int)((5+y0)*yMash));
        
        //picDC->LineTo((int)((7+x0)*xMash),(int)((9+y0)*yMash));
        //координатные оси
        //ось ох
        picDC->MoveTo(POINT(0,yMax));
        picDC->LineTo(POINT(0,yMin));
        // ось оу
        picDC->MoveTo(POINT(xMax,0));
        picDC->LineTo(POINT(xMin,0));
        //picDC->MoveTo(POINT(0,0));
        //picDC->LineTo(POINT(7,0));
 
        // Рисуем сетку
 
        picDC->SelectObject(&penSetka);//контекстному устройству передаём ссылку на penSetka
        for (double x=xMin; x<xMax; x+=1.)
        {
            picDC->MoveTo(POINT(x,yMin));
            picDC->LineTo(POINT(x,yMax));
        }
        for (double y=yMin; y<yMax; y+=1.)
        {
            picDC->MoveTo(POINT(xMin,y));
            picDC->LineTo(POINT(xMax,y));
        }
 
        // рисуем окружности, узлы интегрирования
        char *str = new char [10];//это для перевода цифр в текст
        picDC->SelectObject(&penkrug);
        for (double x=xMin; x<xMax; x+=1.)
        {
            picDC->Ellipse(POINT (x,0),POINT((x+0.01),(0+0.01)));
            picDC->TextOutW(POINT (x,0),(CString)itoa(x,str,10));//передаём нашу координату в строку, макс размер 10.
        }
        for (double y=yMin; y<yMax; y+=1.)
        {
            picDC->Ellipse(POINT (0,y),POINT((0+0.01),(y+0.01)));
            picDC->TextOutW(POINT (0,y),(CString)itoa(y,str,10));//itoa переводит числа в строку
        }
        picDC->SelectObject(&penGrafik);//контекстному устройству передаём ссылку на penGrafik
        picDC->MoveTo(POINT(xMin, f(xMin)));
        for(double x=xMin; x<xMax; x+=0.001 )
        { 
            picDC->LineTo(POINT(x,f(x)));
        }
        
 
    }
}
 
// The system calls this function to obtain the cursor to display while the user drags
//  the minimized window.
HCURSOR CInterpolationDlg::OnQueryDragIcon()
{
    return static_cast<HCURSOR>(m_hIcon);
}
 
 
void CInterpolationDlg::OnBnClickedButton1()
{
    // TODO: Add your control notification handler code here
    // Рисуем график интерполяции
//**********************************************************************************************
 
    UpdateData(true);
    
//  *************************************************************************************************
    OnPaint();
    int n=t;
    picDC->SelectObject(&penInterpolation);
    picDC->MoveTo(POINT(0, F(0,n)));
    for(double x=0.01; x<=1.2;x+=0.04)
    { 
        picDC->LineTo(POINT(x,F(x,n)));
    }
    FILE *File;
    File=fopen("pogr.txt", "w");
    char str[]="Результаты вычислений \n";
    fprintf (File,"%s%s%s%s","x   \t\t", "f(x) \t \t", "F(x) \t \t"," pogr \t \t\n");
    for(double x=0.1; x<=1.1;x+=0.25)
    {
        fprintf(File,"%f\t%f\t%f\t%f\n",x,f(x),F(x,n),pogr(x,n));
    }
    fclose(File);
}

double * x_arg = new double [n+n+2];        //массив аргументов функции в узлах интерполяции
    //cout<<endl<<"h:   "<<h<<endl;
 
    int s=0;//тут введена х чтоб умножать h на положит
    for (int i=-n; i<=n+1; ++i) 
    {
        x_arg[i+n]=a+s*h; s++;// заполнение массива аргументов от -n до n
    //  cout<<x_arg[i]<<"   ";
    }

double ** razn= new double *[n+n+2];    //создание динамического массива для таблицы конечных разностей
    // [n+n+2]столько у нас строк в таблице
    int k=0;
    for (int i=-n;i<=(n+1);++i)
    {
        razn[i+n] = new double [n+n+2+k];
        --k;// выделяем память для каждой строки определенное кол столбцов [n+n+1+k] 
        //к- уменьшает число столбцов в строчке на 1, т.к.
        //таблица разностей имеет диагональный вид
    }

for (int i=-n;i<=(n+1);++i)
{
delete razn[i]
}
delete [] razn;