rowEchlon.cpp

#include<iostream>
using namespace std;
int findleading(float m[3][4],int ro,int r){
        int i;
        for(i=ro;i<r;i++){
            if(m[i][ro]!=0)
                return i;
            }
}
int main(){
    int ro,i,r=3,c=4,l,j;
    float d,instead;
    float m[3][4]={
    {1,3,2,4},
    {2,6,0,-1},
    {1,7,-1,4}
    };

    for(ro=0;ro<c;ro++){
        //find leading
        l=findleading(m,ro,r);
        //change leading to 1
        d=m[l][ro];
        for(i=0;i<c;i++){
            m[l][i]/=d;
        }
        //show matrix
//      for(i=0;i<r;i++){
//          for(j=0;j<c;j++){
//              cout<<" "<<m[i][j];
//          }
//          cout<<endl;
//      }
        //change position of 1
        for(i=0;i<c;i++){
            instead=m[ro][i];
            m[ro][i]=m[l][i];
            m[l][i]=instead;
        }
        //show matrix
//      for(i=0;i<r;i++){
//          for(j=0;j<c;j++){
//              cout<<" "<<m[i][j];
//          }
//          cout<<endl;
//      }
        //change other to 0
        for(i=ro;i<r;i++){
            d=m[i][ro];
            if(i!=ro){
                for(j=0;j<c;j++){
                    m[i][j]-=(d*m[ro][j]);
             }
            }
        }
        //show matrix
//      for(i=0;i<r;i++){
//          for(j=0;j<c;j++){
//              cout<<" "<<m[i][j];
//          }
//          cout<<endl;
//      }
    }
        for(i=0;i<r;i++){
            for(j=0;j<c;j++){
                cout<<" "<<m[i][j];
            }
            cout<<endl;
        }
}