#include<stdio.h>

#include<math.h>

float determinant(float[][25],float);

void cofactor(float[][25],float);

void transpose(float[][25],float[][25],float);

void inverse();

void matrixmultiply();

void multireadmatrix(int[][10],int[][10],int,int,int,int);

void multishowmatrix(int[][10],int,int);

void multiprocessmatrix(int[][10],int[][10],int[][10],int,int,int,int);

void matrixadd(float[][25],float[][25],int,int);

void matrixsub(float[][25],float[][25],int,int);

void matrix_output(float[][25],int,int);

void matrixaddsub();

void determinant_input();

void transpose_input();

void main()

{

int a;char str;

printf("--------------WELCOME TO MATRIX CALCULATION--------------\n");

printf("::::1.Calulate Addition or Subtraction of matrix::::::::\n");

printf("::::2.Calculate Determinant of matrix::::\n");

printf("::::3.Calculate Inverse of a matrix::::::\n");

printf("::::4.Calculate Multiplication of matrix:\n");

printf("::::5.Calculate Transpose of matrix::::::\n");

printf("Enter any choices from 1 to 5::\t");

scanf("%d",&a);

switch (a)

{

case 1:

matrixaddsub();

break;

case 2:

determinant_input();

break;

case 3:

inverse();

break;

case 4:

matrixmultiply();

break;

case 5:

transpose_input();

break;

default:

printf("Invalid choice\n");

printf("Do you want to Re-enter.\n");

printf("Press Y to continue or Press N to exit.\t");

scanf("%c",&str);

if (str=='y'||str=='Y')

{

main();

}

break;

}

}

void transpose_input()

{

int i,j,n;float a[25][25],b[25][25];

printf("Calculation of Matrix Transpose\n");

printf("Enter the dimension of square Matrix:\t");

scanf("%d",&n);

printf("\nEnter the Matrix element of %dX%d matrix\n",n,n);

for(i=0;i<n;i++)

{

for(j=0;j<n;j++)

scanf("%f",&a[i][j]);

}

for (i=0;i<n;i++)

{

for (j=0;j<n;j++)

{

b[i][j]=a[j][i];

}

}

printf("The transpose of a matrix is:\n");

matrix_output(b,n,n);

}

void determinant_input()

{

float a[25][25],k,d;

int i,j;

printf("\n Calculation of determinant of matrix\n\n");

printf("Enter the order of the Square Matrix : ");

scanf("%f",&k);

printf("Enter the elements of %.0fX%.0f Matrix : \n",k,k);

for (i=0;i<k;i++)

{

for (j=0;j<k;j++)

{

scanf("%f",&a[i][j]);

}

}

d=determinant(a,k);

printf("Determinant of matrix is %f",d);

}

void matrixadd(float a[][25],float b[][25],int m, int n)

{

int i,j;float integral[25][25];

for(i=0;i<m;i++)

{

for(j=0;j<n;j++)

integral[i][j]=a[i][j] +b[i][j];

}

printf("The addition of the two matrix is \n");

matrix_output(integral,m,n);

}

void matrixsub(float a[][25],float b[][25],int m, int n)

{

int i,j;float integral[25][25];

for(i=0;i<m;i++)

{

for(j=0;j<n;j++)

integral[i][j]=a[i][j] - b[i][j];

}

printf("The subtraction of the two matrix is \n");

matrix_output(integral,m,n);

}

void matrix_output(float dis[][25],int m,int n)

{

int i,j;

for(i=0;i<m;i++)

{

for(j=0;j<n;j++)

{

printf("%10f",dis[i][j]);

}

printf("\n");

}

}

void matrixaddsub()

{

int i,j,m,n,cho; float a[25][25],b[25][25];

printf("Enter the dimension of Matrix as MXN");

scanf("%d %d",&m,&n);

printf("Enter the element of 1st Matrix\n");

for(i=0;i<m;i++)

{

for(j=0;j<n;j++)

{

scanf("%f",&a[i][j]);

}

}

printf("Enter the element of 2nd Matrix\n");

for(i=0;i<m;i++)

{

for(j=0;j<n;j++)

{

scanf("%f",&b[i][j]);

}

}

printf("1.Calculate Matrix Addition.\n");

printf("2.Calculate Matrix Subtraction.\n");

scanf("%d",&cho);

switch(cho)

{

case 1:

matrixadd(a,b,m,n);

break;

case 2:

matrixsub(a,b,m,n);

break;

}

}

void inverse()

{

float a[25][25],k,d;

int i,j;

printf("-------------------------------------------------------------\n");

printf("-------------------------------------------------------------\n");

printf("-------------------------------------------------------------\n");

printf("\n C Program to find inverse of Matrix\n\n");

printf("Enter the order of the Matrix : ");

scanf("%f",&k);

printf("Enter the elements of %.0fX%.0f Matrix : \n",k,k);

for (i=0;i<k;i++)

{

for (j=0;j<k;j++)

{

scanf("%f",&a[i][j]);

}

}

d=determinant(a,k);

printf("Determinant of the Matrix = %f",d);

if (d==0)

printf("\nInverse of Entered Matrix is not possible\n");

else

cofactor(a,k);

printf("\n\n**** Thanks for using the program!!! ****");

}

//For calculating Determinant of the Matrix

float determinant(float a[25][25],float k)

{

float s=1,det=0,b[25][25];

int i,j,m,n,c;

if (k==1)

{

return (a[0][0]);

}

else

{

det=0;

for (c=0;c<k;c++)

{

m=0;

n=0;

for (i=0;i<k;i++)

{

for (j=0;j<k;j++)

{

b[i][j]=0;

if (i != 0 && j != c)

{

b[m][n]=a[i][j];

if (n<(k-2))

n++;

else

{

n=0;

m++;

}

}

}

}

det=det + s * (a[0][c] * determinant(b,k-1));

s=-1 * s;

}

}

return (det);

}

void cofactor(float num[25][25],float f)

{

float b[25][25],fac[25][25];

int p,q,m,n,i,j;

for (q=0;q<f;q++)

{

for (p=0;p<f;p++)

{

m=0;

n=0;

for (i=0;i<f;i++)

{

for (j=0;j<f;j++)

{

if (i != q && j != p)

{

b[m][n]=num[i][j];

if (n<(f-2))

n++;

else

{

n=0;

m++;

}

}

}

}

fac[q][p]=pow(-1,q + p) * determinant(b,f-1);

}

}

transpose(num,fac,f);

}

//Finding transpose of matrix

void transpose(float num[25][25],float fac[25][25],float r)

{

int i,j;

float b[25][25],inverse[25][25],d;

for (i=0;i<r;i++)

{

for (j=0;j<r;j++)

{

b[i][j]=fac[j][i];

}

}

d=determinant(num,r);

for (i=0;i<r;i++)

{

for (j=0;j<r;j++)

{

inverse[i][j]=b[i][j] / d;

}

}

printf("\n\n\nThe inverse of matrix is : \n");

for (i=0;i<r;i++)

{

for (j=0;j<r;j++)

{

printf("\t%f",inverse[i][j]);

}

printf("\n");

}

}

void matrixmultiply()

{

int matrix1[10][10],matrix2[10][10],m,n,p,q,multi[10][10];

printf("enter dimension of matrix1 MXN\n");

scanf("%d %d",&m,&n);

printf("enter dimension of matrix2 PXQ\n");

scanf("%d %d",&p,&q);

if (n==p)

{

printf("\nEnter elements of matrices\n");

multireadmatrix(matrix1,matrix2,m,n,p,q);

multiprocessmatrix(matrix1,matrix2,multi,m,n,p,q);

printf("\nMatrix multiplication is ::::\n");

multishowmatrix(multi,m,q);

}

else

{

printf("Matrix multiplication not possible");

}

}

void multishowmatrix(int multi[][10],int n,int p)

{

int i,j;

for(i=0;i<n;i++)

{

for(j=0;j<p;j++)

{

printf("%10d ",multi[i][j]);

}

printf("\n");

}

}

void multireadmatrix(int matrix1[][10],int matrix2[][10],int m,int n,int p,int q)

{

int i,j;

printf("\nEnter 1st Matrix\n");

for(i=0;i<m;i++)

{

for(j=0;j<n;j++)

scanf("%d",&matrix1[i][j]);

}

printf("\nEnter 2nd Matrix\n");

for(i=0;i<p;i++)

{

for(j=0;j<q;j++)

scanf("%d",&matrix2[i][j]);

}

}

void multiprocessmatrix(int matrix1[][10],int matrix2[][10],int multi[][10],int m, int n, int p, int q)

{

int i,j,k;

for (i=0;i<m;i++)

{

for(j=0;j<q;j++)

{

multi[i][j]=0;

for(k=0;k<n;k++)

multi[i][j]=multi[i][j]+(matrix1[i][k]*matrix2[k][j]);

}

}

}