Куда уходят поезда ? Вывод значений матрицы

    
                std::cout << "  [ " << (A[i][j]) << "  ]       ";

// Numerical.cpp: определяет точку входа для консольного приложения.
//
 
#include "stdafx.h"
 
#include <iostream>
#include <fstream>
 
using namespace std;
 
 
void Create(double ** &A,int n);
void Destroy(double ** &A,int n);
 
void Show(double ** A,int n);
 
void SwapString(double ** &A,int i,int k);
void SwapColumn(double ** &A,int n,int k,int l);
 
double ** &E(double ** &A,int n); // make Identity matrix
 
void JS(double ** A,double ** &A_inv,int n)
{
    double max = A[0][0];
    int imax = 0,jmax = 0;
 
    for(int i = 0;i < n;i++)
        for(int j = 0;j < n;j++)
            if(A[i][j] > max)
            {
                imax = i;
                jmax = j;
                max = A[i][j];
            }
 
//  SwapString(A,imax,0);
//  SwapColumn(A,n,jmax,0);
 
    for(int i = 0;i < n;i++)
    {
 
        if(A[i][i] != 0)
        {
            double t = 1./A[i][i];
            for(int j = i;j < n;j++)
            {
                A[i][j] = A[i][j]*t;
            }
 
            for(int j = 0;j < n;j++)
            {
                A_inv[i][j] = A_inv[i][j]*t;
            }
 
//          system("pause");
            for(int k = 0;k < i;k++)
            {
                double t = A[k][i];
                for(int l = i;l < n;l++)
                {
                    A[k][l] = A[k][l] - A[i][l]*t;
                }
 
                for(int l = 0;l < n;l++)
                {
                    A_inv[k][l] = A_inv[k][l] - A_inv[i][l]*t;
                }
            }
 
            for(int k = i + 1;k < n;k++)
            {
                double t = A[k][i];
                for(int l = i;l < n;l++)
                {
                    A[k][l] = A[k][l] - A[i][l]*t;
                }
 
                for(int l = 0;l < n;l++)
                {
                    A_inv[k][l] = A_inv[k][l] - A_inv[i][l]*t;
                }
            }
        }
    }
    
//  SwapString(A,imax,0);
//  SwapColumn(A,n,jmax,0);
    
//  SwapColumn(A_inv,n,imax,0);
//  SwapString(A_inv,jmax,0);
        
}
 
double ** &Fill(double ** &A,int n) // make Identity matrix
{
    for(int i = 0;i < n;i++)
    {
        for(int j = 0;j < n;j++)
            {
                A[i][j] = rand()%32;// = rand()%20 * (1.0/2.0)
            }
//      std::cout << std::endl;
    }
    return A;
}
 
void MultL(double ** A,double ** B,double ** &C,int n)
{
    for(int i = 0;i < n;i++)
        for(int j = 0;j < n;j++)
        {
            double s = 0;
            for(int k = 0;k < n;k++)
            {
                s = s + A[i][k] * B[k][j];
            }
            C[i][j] = s;
        }
    
}
 
void Read(double ** &A,int n,char * fname)
{
    std::ifstream fin(fname);
    if(fin)
    {
        for(int i = 0;i < n;i++)
            for(int j = 0;j < n;j++)
            {
                int t;
                fin >> t;
                A[i][j] = (double)t;
            }
        fin.close();
    }
}
 
int main()
{
    const int n = 4;
    double ** p, **p_inv,**res; 
    Create(p,n);
    Create(res,n);
 
    Create(p_inv,n);
        E(p_inv,n);
 
/*  p[0][0] = 0;
    p[0][1] = -2;
    p[1][0] = 0.5;
    p[1][1] = 0;*/
//  Fill(p,n);
 
        Read(p,n,"input");
 
        Show(p,n);
 
        JS(p,p_inv,n);
//          E(p,n);
//      JS(p_inv,p,n);
 
        cout << endl;
 
        Read(p,n,"input");
 
        MultL(p,p_inv,res,n);
        Show(res,n);
        cout << endl;
        Show(p_inv,n);
 
    Destroy(p,n);
    system("pause");
    return 0;
}
 
 
void Create(double ** &A,int n)
{
    A = new double*[n];
    for(int i = 0;i < n;i++)
        A[i] = new double[n];
}
 
void Destroy(double ** &A,int n)
{
    for(int i = 0;i < n;i++)
        delete [] A[i];
    delete [] A;
}
 
double ** &E(double ** &A,int n)
{
    for(int i = 0;i < n;i++)
        for(int j = i + 1;j < n;j++)
            A[i][j] = 0;
 
    for(int i = 0;i < n;i++)
        A[i][i] = 1;
 
    for(int i = 0;i < n;i++)
        for(int j = 0;j < i;j++)
            A[i][j] = 0;
    return A;
}
 
void SwapString(double ** &A,int i,int k)
{
    if(i == k) return;
    double * p = A[i];
        A[i] = A[k];
        A[k] = p;
}
 
void SwapColumn(double ** &A,int n,int k,int l)
{
    if(k == l) return;
 
    for(int i = 0;i < n;i++)
    {
        double t = A[i][k];
        A[i][k] = A[i][l];
        A[i][l] = t;
    }
}
 
void Show(double ** A,int n)
{
    for(int i = 0;i < n;i++)
    {
        for(int j = 0;j < n;j++)
            {
                std::cout.width(10);
                std::cout << "  [ " << (A[i][j]) << "  ]       ";
            }
        std::cout << std::endl;
    }
}

   Read(p,n,"input");

double t = 1./A[i][i];