void print( int n,int right, int x,int y )
{
int i,j;
if( n == 3 )
{
for( i = 0; i < 9; i++ )
for( j = 0; j < 9; j++ )
fractal[x+i][y+j] = out3[right][i][j];
}
else
{
print(n-1,0,x,y);//0
fill_b(n-1,x,y+width[n-1]);//1
// for( i = 0; i < power[n-2]; i++ )
// fractal[x+i][y+width[n-1]*2] = ' ';
print(n-1,right,x,y+width[n-1]*2);//2
fill_b(n-1,x+power[n-2],y);//3
print(n-1,right,x+power[n-2],y+width[n-1]);//4
if( !right )//5
{
// for( i = power[n-2]; i < power[n-2]*2; i++ )
// fractal[x+i][y+width[n-1]*2] = ' ';
fill_b(n-1,x+power[n-2],y+width[n-1]*2);
}
print(n-1,0,x+power[n-2]*2,y);//6
fill_b(n-1,x+power[n-2]*2,y+width[n-1]);//7
//for( i = power[n-2]*2; i < power[n-1];i++ )
// fractal[x+i][y+width[n-1]*2] = ' ';
print(n-1,right,x+power[n-2]*2,y+width[n-1]*2);//8
}
}
int main(int argc, char *argv[])
{
int i, j;
matrix_double A;
double B[A.nrows];
A = alloc_matrix_double(3, 3);
srand(100);
for (i = 0; i < A.nrows; i++) {
for (j = 0; j < A.ncols; j++) {
A.data[i][j] = (double) (rand() % 100);
}
}
fill_b(A, B);
printf("********************************************\n");
printf("Trying a random system with solution 1, 2, 3\n");
printf("********************************************\n");
solve(A, B);
A.data[0][0] = 1;
A.data[0][1] = 2;
A.data[0][2] = 3;
A.data[1][0] = 4;
A.data[1][1] = 5;
A.data[1][2] = 6;
A.data[2][0] = 7;
A.data[2][1] = 8;
A.data[2][2] = 9;
fill_b(A, B);
printf("\n*****************************\n");
printf("Trying with a singular matrix\n");
printf("*****************************\n");
solve(A, B);
return 0;
}
