GAUSS Jordan

/*..................GAUSS JORDAN.............
...........AUTHOR:-SUDIPTA KUMAR PAIK......
......................N4_11................*/
# define eqn 20

main()
{ int q,w,m,i,j,kkl,ll,n,h,f,k,d,alp=97;
float r[eqn][eqn+1],a[eqn][eqn+1];

do
{
clrscr();

printf("\n Enter [ 0 ] for EXIT \n");
printf("\n Enter No. of Equation ");

scanf("%d",&n);
if(n==0)
{ printf("\n\n\n\t\tThank you for use this program\n");
printf("\n\t\t.............Prepared By SUDIPTA.............. ");
getch();
exit(1);
}
m=n+1;

input:
printf("\n\n [1]-Manual Input ");
printf("\n\n [2]-Random Input ");
printf("\n\n Choice ");
scanf("%d",&kkl);

switch(kkl)
{
case 1:

for(q=1; q<=n; q++)
{ printf("\n\n Eqution %d\n",q);
for(w=1; w<=n; w++)
{ printf("\n Enter %c%d=",alp,w);
scanf("%f",&r[q][w]);
}
printf("\nEnter Constant K%d=", q);
scanf("%f",&r[q][w]);
alp++;
}
break;

case 2:

for(q=1; q<=n; q++)
{ for(w=1; w<=m; w++)
r[q][w]=rand()%(9);
}
break;

default:
printf("\n Wrong Selection");
getch();
goto input;
}

printf("\n\n\n");
for(q=1; q<=n; q++)
{ printf("\n\n");
for(w=1; w<=n; w++)
{ if(r[q][w]<0)
printf("%4.2fX%d",r[q][w],w );
else if(w==1)
printf("%4.2fX%d",r[q][w],w );
else
printf("+%4.2fX%d",r[q][w],w );
}
printf("=%4.2f",r[q][w]);
}

printf("\n\n\n Ú");
for(q=1; q<=8*m; q++)
printf(" ");
printf(" ¿");
printf("\n ³");
for(q=1; q<=8*m; q++)
printf(" ");
printf(" ³");

for(q=1; q<=n; q++)
{ printf("\n ³");
for(w=1; w<=n; w++)
{ printf("%8.4f",r[q][w]);
}
printf("%8.4f ",r[q][w] );
printf("³\n ³");

for(ll=1; ll<=8*m; ll++)
printf(" ");
printf(" ³");

}
printf("\n À");
for(ll=1; ll<=8*m; ll++)
printf(" ");
printf(" Ù");
getch();

for(k=1,j=1; j<=n; j++,k++)
{
for(f=1; f<=n; f++)
{ for(h=1; h<=m; h++)
a[f][h]=r[f][h];
}

for(i=1; i<=m; i++)
r[j][i]=r[j][i]/a[k][k];

for(f=1; f<=n; f++)
{ for(h=1; h<=m; h++)
a[f][h]=r[f][h];
}

for(i=1; i<=n; i++)
{ if(i!=j)
{ for(d=1; d<=m; d++)
r[i][d]=r[i][d]-a[i][j]*r[j][d];
}
}

}

printf("\n\n\n\n\tResult Determinant.........\n");
printf("\n\n Ú");
for(q=1; q<=8*m; q++)
printf(" ");
printf(" ¿");
printf("\n ³");
for(q=1; q<=8*m; q++)
printf(" ");
printf(" ³");

for(q=1; q<=n; q++)
{ printf("\n ³");
for(w=1; w<=n; w++)
printf("%8.4f",r[q][w]);
printf("%8.4f ",r[q][w] );
printf("³\n ³");

for(ll=1; ll<=8*m; ll++)
printf(" ");
printf(" ³");

}
printf("\n À");
for(ll=1; ll<=8*m; ll++)
printf(" ");
printf(" Ù");

printf("\n\n Result.........\n");
for(q=1; q<=n; q++)
printf("\n\n X%d=%.5f",q,r[q][m]);

getch();
}while(1);

}

No comments:

Post a Comment