void balanc(double a[6][6],int n)
{
	double radix,sqrdx,c,r,f,g,s;
	int last,i,j;
    radix = 2.0;
    sqrdx = radix *radix;
loop1:   last = 1;
	for (i = 1; i<=n; i++)
	{
        c = 0.0;
        r = 0.0;
        for (j = 1; j<=n; j++)
		{
            if (j != i)
			{
                c = c + fabs(a[j][i]);
                r = r + fabs(a[i][j]);
            }
        }
        if ((c != 0.0) && (r != 0.0))
		{
            g = r / radix;
            f = 1.0;
            s = c + r;
loop2:      if (c < g)
			{ 
                f = f * radix;
                c = c * sqrdx;
                goto  loop2;
			}
            g = r * radix;
loop3:      if (c > g)
			{ 
                f = f / radix;
                c = c / sqrdx;
                goto  loop3;
			}
            if (((c + r) / f) < (0.95 * s))
			{
                last = 0;
                g = 1.0 / f;
                for (j = 1 ; j<=n; j++)
				{
                    a[i][j] = a[i][j] * g;
                }
                for (j = 1 ; j<=n; j++)
				{
                    a[j][i] = a[j][i] * f;
                }
            }
        }
    }
    if (last == 0) goto loop1;
}