Нужно чтобы в матрице столбец менялся на столбец свободных членов

#include<iostream.h>
#include<math.h>
#include<conio.h>
int main()
{
clrscr();
float a11,a12,a13,a21,a22,a23,a31,a32,a33,b1,b2,b3,D1,D2,D3,D,x1,x2,x3;
cout<<" Vvedite koefichienti:\n";
cout<<" a11=";
cin>>a11;
cout<<" a12=";
cin>>a12;
cout<<" a13=";
cin>>a13;
cout<<" a21=";
cin>>a21;
cout<<" a22=";
cin>>a22;
cout<<" a23=";
cin>>a23;
cout<<" a31=";
cin>>a31;
cout<<" a32=";
cin>>a32;
cout<<" a33=";
cin>>a33;
cout<<" Vvedite svobodnie chleny:\n";
cout<<" b1=";
cin>>b1;
cout<<" b2=";
cin>>b2;
cout<<" b3=";
cin>>b3;
D=a11*(a22*a33-a23*a32)-a12*(a21*a33-a23*a31)+a13*(a21*a32-a22*a31);
if (D==0)
cout<<"beskonechnoe mnozhestvo resheniy";
else{
D1=b1*(a22*a33-a23*a32)-a12*(b2*a33-a23*b3)+a13*(b2*a32-a22*b3);
D2=a11*(b2*a33-a23*b3)-b1*(a21*a33-a23*a31)+a13*(a21*b3-b2*a31);
D3=a11*(a22*b3-a32*b2)-a12*(a21*b3-b2*a31)+b1*(a21*a32-a22*a31);
x1=D1/D;
x2=D2/D;
x3=D3/D;
cout<<"\n Opredeliteli:\n";
cout<<" D="<<D<<"; ";
cout<<" D1="<<D1<<"; " ;
cout<<" D2="<<D2<<"; ";
cout<<" D3="<<D3<<";\n ";
cout<<"\n Korni sistemy uravneniy:\n";
cout<<" x1="<<x1<<"; ";
cout<<" x2="<<x2<<"; ";
cout<<" x3="<<x3<<"; ";
getch();
return 0;
}
}

j=0;
k=0;
for(i=0;i<2;i++){
a[i][j]=d[k];
k++;}