Вывод транспонированной матрицы

#include<iostream.h>
#include<math.h>
int SLAU(double **matrica_a,int n,double *massiv_b,
double *x)
{
int i,j,k,r;
double c,M,max,s, **a, *b;
a=new double *[n];
for(i=0;i<n;i++)
a[i]=new double[n];
b=new double [n];
for(i=0;i<n;i++)
for(j=0;j<n;j++)
a[i][j]=matrica_a[i][j];
for(i=0;i<n;i++)
b[i]=massiv_b[i];
for(k=0;k<n;k++)
{
max=fabs(a[k][k]);
r=k;
for(i=k+1;i<n;i++)
if (fabs(a[i][k])>max)
{
max=fabs(a[i][k]);
r=i;
}
for(j=0;j<n;j++)
{
c=a[k][j]; a[k][j]=a[r][j];
a[r][j]=c;
}
c=b[k];b[k]=b[r];b[r]=c;
for(i=k+1;i<n;i++)
{
for(M=a[i][k]/a[k][k],j=k;j<n;j++)
a[i][j]-=M*a[k][j];
b[i]-=M*b[k];
} }
if (a[n-1][n-1]==0)
if(b[n-1]==0)
return -1;
else return -2;
else
{
for(i=n-1;i>=0;i--)
{
for(s=0,j=i+1;j<n;j++)
s+=a[i][j]*x[j];
x[i]=(b[i]-s)/a[i][i];
}
return 0;}
for(i=0;i<n;i++)
delete [] a[i];
delete [] a;
delete [] b;
}
int main()
{int result,i,j,N;
double **a, *b, *x;
cout<<"N="; cin>>N;
a=new double *[N];
for(i=0;i<N;i++)
a[i]=new double[N];
b=new double [N];
x=new double [N];
cout<<"Input Matrix A"<<endl;
for(i=0;i<N;i++)
for(j=0;j<N;j++)
cin>>a[i][j];
cout<<"Input massiv B"<<endl;
for(i=0;i<N;i++)
cin>>b[i];
result=SLAU(a,N,b,x);
if (result==0)
{ cout<<"Massiv X"<<endl;
for(i=0;i<N;i++)
cout<<x[i]<<"\t";
cout<<endl;
}
else if (result==-1)
cout<<"Great number of Solution";
else if (result==-2)
cout<<"No solution";
for(i=0;i<N;i++)
delete [] a[i];
delete [] a; delete [] b; delete [] x; 
}

#include <iostream>
#include <cmath>
 
using namespace std;
 
int main()
{
const int N = 3;
const int M = 4;
int A[N][M], B[M][N], p;
cout << "Введите матрицу:" << endl;
for (int i = 0; i < N; i++)
{
for (int j = 0; j < M; j++)
{
cin >> A[i][j];
}
} 
cout << "Матрица после транспонирования:" << endl;
p = 0;
while (p < M)
{
for (int i = 0; i < N; i++)
{
B[p][i] = A[i][p];
cout << B[p][i] << " "; 
}
cout << endl;
p++;
}
system("pause");
return 0;
}