Problem sheet problem 4

/* To find Trace and subtraction of two matrices */
#include<stdio.h>
#include<math.h>
int main()
{
int i, j, n, m;
float a[20][20], b[20][20], c[20][20], d[20][20],sum=0,sum1=0,sum2=0,sum3=0;
printf("Enter order of the matrices:\n");
scanf("%d%d", &m, &n);
printf("Enter the elements of matrix A:\n");
for(i=0;i<m;i++)
{
for(j=0;j<n;j++)
scanf("%f", &a[i][j]);
printf("\n");
}
printf("Entered Matrix A:\n");
for(i=0;i<m;i++)
{
for(j=0;j<n;j++)
printf("%10f",a[i][j]);
printf("\n");
}
printf("\nEnter the elements of matrix B:\n");
for(i=0;i<m;i++)
{
for(j=0;j<n;j++)
scanf("%f", &b[i][j]);
printf("\n");
}
printf("Entered Matrix B:\n");
for(i=0;i<m;i++)
{
for(j=0;j<n;j++)
printf("%10f",b[i][j]);
printf("\n");
}
for(i=0;i<m;i++)
{
for(j=0;j<n;j++)
{
if(i==j)
sum=sum+a[i][j];
}
}
    printf("\n\ntrace of the matrix A=%f",sum);
for(i=0;i<m;i++)
{
for(j=0;j<n;j++)
{
if(i==j)
sum1=sum1+b[i][j];
}
}
    printf("\n\ntrace of the matrix B=%f",sum1);
/* Matrix subtraction */
for(i=0;i<m;i++)
{
for(j=0;j<n;j++)
d[i][j]=2*a[i][j]-b[i][j];
}

for(i=0;i<m;i++)
{
for(j=0;j<n;j++)
{
if(i==j)
sum2=sum2+d[i][j];
}
}
    printf("\n\ntrace of the matrix (2A-B)=%f",sum2);
sum3=2*sum-sum1;
printf("\n\n2*Trace A - Trace B=%f",sum3);
if(sum2==sum3)
{
printf("\n\nSo, Trace(2A-B)=2*Trace(A) - Trace(B)");
}
printf("\n\nMatrix subtraction D is\n");
for(i=0;i<m;i++)
{
for(j=0;j<n;j++)
printf("%10f\t", d[i][j]);
printf("\n");
}
return(0);
}