Аппроксимация, постройка графика

#pragma hdrstop
#pragma argsused
#define _USE_MATH_DEFINES
 
#ifdef _WIN32
#include <tchar.h>
#include <iostream>
#include <cmath>
#include <iomanip>
#include <math.h>
using namespace std;
#else
  typedef char _TCHAR;
  #define _tmain main
#endif
 
#include <stdio.h>
void filter(double gk[][2])
{
    double stdv = 1.0;
    double r, s = 2.0 * stdv * stdv;  // Assigning standard deviation to 1.0
    double sum = 0.0;   // Initialization of sun for normalization
    for (int x = -2; x <= 2; x++) // Loop to generate 5x5 kernel
    {
        for (int y = -2; y <= 2; y++)
        {
            r = sqrt(x*x + y*y);
            gk[x + 2][y + 2] = (exp(-(r*r) / s)) / (M_PI * s);
            sum += gk[x + 2][y + 2];
        }
    }
 
    for (int i = 0; i < 10; ++i) // Loop to normalize the kernel
        for (int j = 0; j < 2; ++j)
            gk[i][j] /= sum;
}
 
 
 int _tmain(int argc, _TCHAR* argv[])
{
double gk[10][2];
    for (int i = 0; i < 10; i++){
for (int j = 0; j < 2; j++)
 
        {
            cout << "vedit GK[" << i << "," << j << "]:";
            cin >> gk[i][j];
        }
    }
 
 
    filter(gk); // Function call to create a filter
 
    for (int i = 0; i < 10; ++i)
    {
        for (int j = 0; j < 2; ++j)  // loop to display the generated 5 x 5 Gaussian filter
        {
            cout << gk[i][j] << "\t";
            cout << endl;
        }
    }
return 0;
}