int n=3;
float[][] a;
float[] b;
float c;
int i, j, k;
String sd;

void setup()
{
  a = new float [n][n];
  b = new float [n];
  a[0][0]= 1; a[0][1]= 2; a[0][2]= 4; 
  a[1][0]= 2; a[1][1]= 5; a[1][2]= 2;
  a[2][0]= 4; a[2][1]=-1; a[2][2]= 1;
  b[0]=11; b[1]=3; b[2]=8;
  printab();
  gaussjd();
  printab();  
  printsolution();
}

void gaussjd(){
  for(k=0; k<n; k++) {
    for(i=0; i<n; i++) {
      c = a[i][k]/a[k][k];
      for(j=k; j<n; j++) {
        if((i!=k) && (j==k)) a[i][j]=0;
        else if((i!=k) && (j>k)) a[i][j]-=c*a[k][j];
      }
      if(i!=k) b[i]-=c*b[k];
    }
  }
}

void printab() {
  for(i=0; i<n; i++) {
    for(j=0; j<n; j++) { 
      sd = nfs(a[i][j], 3, 2);
      print(sd+" ");
    }  
    sd = nfs(b[i], 3, 2);
    print(sd);  
    println();
  }
}

void printsolution() {
  for(i=0; i<n; i++) {
    sd = nfs(b[i]/a[i][i], 3, 2);
    print(sd);  
    println();
  }
}