#include <stdio.h>
#include <math.h>

double wyznacznik (double a[][2])
{
return a[0][0]*a[1][1]-a[0][1]*a[1][0];
}
double macierz (double a[][2],double b[], double x[])
{
double w=wyznacznik(a);

double t[2][2];

t[0][0]=b[0];
t[1][0]=b[1];
t[0][1]=a[0][1];
t[1][1]=a[1][1];
double wx=wyznacznik(t);

t[0][0]=a[0][0];
t[1][0]=a[1][0];
t[0][1]=b[0];
t[1][1]=b[1];
double wy=wyznacznik(t);

if (w!=0.0)
{
x[0]=wx/w;
x[1]=wy/w;
return 1;
}
else if (wx==0.0 && wy==0.0)
return 1000;
else return 0;
}

int main()
{
double a[2][2],b[2],x[2];
int i=0,j=0;
printf("Podaj uklad postaci: \n");
printf("a11x+a12y=b1 \n");
printf("a21x+a22y=b2 \n");
for(j=0;j<2;j++)
{
 for (i=0;i<2;i++)
 {
 printf("Podaj a%d %d: \n",j+1,i+1);
 scanf("%lg",&a[j][i]);
 }
 printf("Podaj b%d: \n",j+1);
 scanf("%lg",&b[j]);
}
wyznacznik(a);
macierz(a,b,x);
if (macierz(a,b,x)==1)
{
printf("Układ postaci: \n");
printf("%+gx%+gy=%+g \n",a[0][0],a[0][1],b[0]);
printf("%+gx%+gy=%+g \n",a[1][0],a[1][1],b[0]);
printf("posiada jedno rozwiązanie x=%+g i y=%+g",x[0],x[1]);
}
if (macierz(a,b,x)==1000)
{
printf("Układ postaci: \n");
printf("%+gx%+gy=%+g \n",a[0][0],a[0][1],b[0]);
printf("%+gx%+gy=%+g \n",a[1][0],a[1][1],b[0]);
printf("posiada nieskończenie wiele rozwiązań \n");
}
if (macierz(a,b,x)==0)
{
printf("Układ postaci: \n");
printf("%+gx%+gy=%+g \n",a[0][0],a[0][1],b[0]);
printf("%+gx%+gy=%+g \n",a[1][0],a[1][1],b[0]);
printf("nie posiada rozwiązań");
}
return 0;
}