#include <stdio.h>
#include <stdlib.h>
#include <malloc.h>
#define N = 3
void razn (double *m1, double *m2, double *m3)
{
     int l;
     for (l = 0 ; l<N ;l++)
     m3[l] = m1[l] - m2[l];
     return;
 }

void zvezda(double **m1, double *m2, double *m3)
{
     int i, j;
     double s;
     for (i=0; i<N; i++)
     {
         s=0;
         for (j=0;j<N;j++)
           s = s + m1[i][j]*m2[j];
         m3[i] = s;
     }
     return;
 }
 void shag(double **m1, double *m2,double *m3, double *m4)
 {
  int i, j;
  double s1=0, s2=0;
  for (i=1; i<N; i++)
    s1 = s1 + m1[0][i]*m3[i];
  m2[0] = (m4[0] - s1)/m1[0][0];
  for (i=1; i<(N-1); i++)
  {
      s1=0;
      s2=0;
      for (j=0; j<i; j++)
        s1 = s1+ m1[i][j]*m2[j];
      for (j=i+1; j<N; j++)
        s2 = s2 + m1[i][j]*m3[j];
      if (m1[i][i]!=0)
        m2[i] = (m4[i]- s1 - s2)/m1[i][i];
      else
        printf ("\n error");
    }
  s2 = 0;
  for (i=0; i<(N-1); i++)
    s2 = s2+ m1[N-1][i]*m2[i];
  m2[N-1] = (m4[N-1] - s2)/a[N-1][N-1];
  return;
  }
 int main ()
{

        int i,j;

        double **A;
        double *b, *xk1, *r, *xk, *t, *e;
        double c, sum, help, bam;
A = (double **) malloc (N * sizeof(double));
 for(i = 0; i < N; i++)
 A[i] = (double *) malloc (N * sizeof(double));
b = (double *) malloc (N * sizeof(double));
xk1 = (double *) malloc (N * sizeof(double));
r = (double *) malloc (N * sizeof(double));
xk = (double *) malloc (N * sizeof(double));
t = (double *) malloc (N * sizeof(double));
e = (double *) malloc (N * sizeof(double));

for (i=0; i<N; i++)
{
for (j=0; j<N; j++)
scanf (" %lf ", &A[i][j]);

}

for (i=0; i<N; i++)
scanf (" %lf ", &b[i]);

for (i=0; i<N; i++)
scanf (" %lf ", &xk[i]);



for (i=0; i<N; i++)
{
    e[i] = 0;
    t[i] = 0;
}


zvezda(A, xk, t);
razn (t,b, r);
bam = 1;
 for (i=0; i<N; i++)
 {
     if  (r[i] = e[i])
     bam = bam*1.0;
     else
     bam = bam*0.0;
 }

while (bam != 1.0)
{
      shag(A, xk1, xk, b);
      for (i=0; i<N; i++)
      xk[i] = xk1[i];
      for (i=0; i<N; i++)
      t[i] = e[i];
      zvezda (A, xk1, t);
      razn (t,b, r);
      bam = 1;
      for (i=0; i<N; i++)
 {
     if  (r[i] = e[i])
     bam = bam*1.0;
     else
     bam = bam*0.0;
 }

}
printf ("\n sha budet otvet \n");
for (i=0; i<N; i++)
printf (" %lf ", xk[i]);


return (0);
}