#include<stdio.h>

void main()
{
  const int DIM=2;

  float a[DIM][DIM], b[DIM][DIM], c[DIM][DIM];

  int i, j, k;

  float r;

  printf("\n Enter first matrix a:\n");

  for(i=0; i<DIM; i++)
    {
      for(j=0; j<DIM; j++)
	{
	  printf("\n Enter a[%d][%d]: ", i,j);
	  scanf("%f", &a[i][j]);

	}

    }

  printf("\n Enter second matrix b:\n");

  for(i=0; i<DIM; i++)
    {
      for(j=0; j<DIM; j++)
	{
	  printf("\n Enter b[%d][%d]: ", i,j);
	  scanf("%f", &b[i][j]);

	}

    }

  for(i=0; i<DIM; i++)
    {
      for(j=0; j<DIM; j++)
	{
	  c[i][j]=0;
	  /* Now, need to compute c[i][j] */
	  for(k=0; k<DIM; k++)
	    c[i][j]+=a[i][k]*b[k][j];

	}

    }


 printf("\n Their product is :\n");

  for(i=0; i<DIM; i++)
    {
      for(j=0; j<DIM; j++)
	{
	  printf("%f ", c[i][j]);

	}
      printf("\n");

    }


}