Cramers' Rule

#include

int i,j;
char display(i)
{
if(j==0) return 'x';
if(j==1) return 'y';
if(j==2) return 'z';
return 0;
}

int main()
{
float a[3][4],d,x,y,z;
system("clear");

for(i=0;i<3;i++)
{
printf("\n\tEnter the values of equation %d\n",i+1);
for(j=0;j<4;j++)
{
if(j!=3) printf("\tEnter the coefficent of %c = ",display(j));
else printf("\tEnter the value of equat = ");
scanf("%f",&a[i][j]);
}
}

d=((a[0][0]*(a[1][1]*a[2][2]-a[1][2]*a[2][1]))-(a[0][1]*(a[1][0]*a[2][2]-a[2][0]*a[1][2]))+(a[0][2]*(a[1][0]*a[2][1]-a[1][1]*a[2][0])));
x=((((a[0][3]*(a[1][1]*a[2][2]-a[1][2]*a[2][1]))-(a[0][1]*(a[1][3]*a[2][2]-a[2][3]*a[1][2]))+(a[0][2]*(a[1][3]*a[2][1]-a[1][1]*a[2][3]))))/d);
y=((((a[0][0]*(a[1][3]*a[2][2]-a[1][2]*a[2][3]))-(a[0][3]*(a[1][0]*a[2][2]-a[2][0]*a[1][2]))+(a[0][2]*(a[1][0]*a[2][3]-a[1][3]*a[2][0]))))/d);
z=((((a[0][0]*(a[1][1]*a[2][3]-a[1][3]*a[2][1]))-(a[0][1]*(a[1][0]*a[2][3]-a[2][0]*a[1][3]))+(a[0][3]*(a[1][0]*a[2][1]-a[1][1]*a[2][0]))))/d);
printf("\tThe solution is\n\tDelta = %f,\n\tX= %f,\n\tY= %f,\n\tZ= %f\n",d,x,y,z);

return 0;
}