Dev - C++ - Solución de sistemas lineales con matrices

Hola a todos. Tengo un problema con mi variable j que debería tomar valores 1 y 2. Pero solo toma el valor 1. Me refiero al for(int j=j+i-1; j<n; j++){. Espero sus respuestas por favor. Gracias

include <iostream>
#include <math.h>
using namespace std;
int main()
{
float a[2][2], b[2], Ratio;
int n=2;
for(int i=1; i<=n; i++){
for(int j=1; j<=n; j++){
cout<<"introduzca el valor a["<<i<<"]["<<j<<"] ";cin>>a[i][j];
}
}
for(int i=1; i<=n; i++){
cout<<"introduzca el valor b["<<i<<"]";cin>>b[i];
}
for(int i=2; i>=1; i-- ){
if((a[i][i])=0){
cout<<"no hay solución"<<endl;
}
else{
for( int j=j+i-1; j<n; j++){
Ratio= a[j][i]/a[i][i];
cout<<j<<endl;
for(int k=2; k>=1; k-- ){
a[j][k]= a[j][k]-(Ratio*a[i][k]);
cout<<"a["<<j<<"]["<<k<<"] = "<<a[j][k]<<endl;
}
}
}
}
return 0;
}

Hola, te dejo una pista: los vectores y matrices comienzan en cero:
Ejemplo: