Разбить код на функции

#include <conio.h>
#include <time.h>
#include <stdlib.h>
#include <math.h>
#include <vector>
#include <iostream>
#pragma argsused
int main(int argc, char* argv[])
{
    using namespace std;
    srand((unsigned)time(NULL));
    int n = 10;
    int m = 10;
    int r_limit = 10;
    vector<vector<double> > matrix;
    int i, j;
    for (i = 0; i < n; i++) {
        vector<double> v;
 
        for (j = 0; j < m; j++)
            v.push_back((rand() % 2) ? (static_cast<double>(rand() % r_limit) + 1. / static_cast<double>(rand() % r_limit + 1.)) :
                -(static_cast<double>(rand() % r_limit) + 1. / static_cast<double>(rand() % r_limit + 1.)));
        matrix.push_back(v);
    }
    cout << "Initial matrix:" << endl;
    for (i = 0; i < n; i++) {
        for (j = 0; j < m; j++)
            cout << matrix[i][j] << " ";
        cout << endl;
    }
    vector<vector<double> > smooth_matrix;
    for (i = 0; i < n; i++) {
        vector<double> v;
        for (j = 0; j < m; j++) {
            double avg = 0.;
            int elements_count = 0;
            if (j - 1 >= 0) {
                if (i - 1 >= 0) {
                    avg += matrix[i - 1][j - 1];
                    elements_count++;
                }
                avg += matrix[i][j - 1];
                elements_count++;
                if (i + 1 < n) {
                    avg += matrix[i + 1][j - 1];
                    elements_count++;
                }
            }
            if (i - 1 >= 0) {
                avg += matrix[i - 1][j];
                elements_count++;
            }
            if (i + 1 < n) {
                avg += matrix[i + 1][j];
                elements_count++;
            }
            if (j + 1 < m) {
                if (i - 1 >= 0) {
                    avg += matrix[i - 1][j + 1];
                    elements_count++;
                }
                avg += matrix[i][j + 1];
                elements_count++;
                if (i + 1 < n) {
                    avg += matrix[i + 1][j + 1];
                    elements_count++;
                }
            }
            if (elements_count != 0)
                avg /= static_cast<double>(elements_count);
 
            v.push_back(avg);
        }
        smooth_matrix.push_back(v);
    }
    cout << endl << "Smooth matrix:" << endl;
    for (i = 0; i < n; i++) {
        for (j = 0; j < m; j++)
            cout << smooth_matrix[i][j] << " ";
        cout << endl; }
    double summ = 0;
    for (i = 0; i < n; i++)
        for (j = 0; j < i; j++)
            summ += fabs(smooth_matrix[i][j]);
    cout << endl << "Summ of elements modules at main diagonal in smooth matrix is: " << summ << endl;
    getch();
    return 0;}