ガウスフィルタのc++,matlab実現

c++バージョン

#include 
using namespace std;
int gaussian(double **gaus, const int size, const double sigma);
int main()
{
	int size = 4;
	double **gaus = new double *[size];
	for (int i = 0; i < size; i++) {
		gaus[i] = new double[size];  //      
	}
	gaussian(gaus,4, 1);
	return 0;
	 
}
int gaussian(double **gaus, const int size, const double sigma)
{
	const double PI = 4*atan(1.0);
	double radius;
	double center = (size+1) / 2;
	double sum = 0;
	for (int i = 0; i < size; i++)
	{
		for (int j = 0; j < size; j++) {
			radius = ((i + 1 - center)*(i + 1 - center) + (j + 1 - center)*(j + 1 - center));
			gaus[i][j] = (1/(2 * PI *sigma*sigma))* exp(-(radius / (2 * sigma*sigma)));
			sum += gaus[i][j];
		}
	}
	cout << sum << endl;
	for (int i = 0; i < size; i++) {
		for (int j = 0; j < size; j++) {
			gaus[i][j] /= sum;
			cout << gaus[i][j] << "	  ";
		}
		cout << endl;
	}
	return 0;
}

matlabバージョン

#include 
using namespace std;
int gaussian(double **gaus, const int size, const double sigma);
int main()
{
	int size = 4;
	double **gaus = new double *[size];
	for (int i = 0; i < size; i++) {
		gaus[i] = new double[size];  //      
	}
	gaussian(gaus,4, 1);
	return 0;
	 
}
int gaussian(double **gaus, const int size, const double sigma)
{
	const double PI = 4*atan(1.0);
	double radius;
	double center = (size+1) / 2;
	double sum = 0;
	for (int i = 0; i < size; i++)
	{
		for (int j = 0; j < size; j++) {
			radius = ((i + 1 - center)*(i + 1 - center) + (j + 1 - center)*(j + 1 - center));
			gaus[i][j] = (1/(2 * PI *sigma*sigma))* exp(-(radius / (2 * sigma*sigma)));
			sum += gaus[i][j];
		}
	}
	cout << sum << endl;
	for (int i = 0; i < size; i++) {
		for (int j = 0; j < size; j++) {
			gaus[i][j] /= sum;
			cout << gaus[i][j] << "	  ";
		}
		cout << endl;
	}
	return 0;
}