Matrix Insertion, Display, Addition, Subtraction and Multiplication

#include

int main(void)
{
        int a[3][3],b[3][3],c[3][3];
        int i,j;
        printf("\a \a \a \a enter element of first matrix\n");
        for(i=0;i<3;i++)
        {
                for(j=0;j<3;j++)
                {
                        printf("\nenter a[%d][%d] element=",i,j);
                        scanf("%d",&a[i][j]);
                }
        }
        printf("\nenter element of second matrix");
        for(i=0;i<3;i++)
        {
                for(j=0;j<3;j++)
                {
                        printf("\nenter b[%d][%d] element=",i,j);
                        scanf("%d",&b[i][j]);
                }
        }
        printf("\n\nmatrix a is....\n");
        for(i=0;i<3;i++)
        {
                for(j=0;j<3;j++)
                {
                        printf("%d\t",a[i][j]);
                }
                printf("\n");
        }
        printf("content of matrix b is=\n");
        for(i=0;i<3;i++)
        {
                for(j=0;j<3;j++)
                {
                        printf("%d\t",b[i][j]);
                }
                printf("\n");
        }
        printf("\n welcom to addition\n");
        for(i=0;i<3;i++)

  {

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

                {

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

                }

        }

        printf("\n\nmatrix after addition is \n=");

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

        {

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

                {

                        printf("%d\t",c[i][j]);

                }

                printf("\n");

        }

        printf("\n welcom to substraction\n");

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

        {

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

                {

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

                }

        }

        printf("\n\nmatrix after substraction is \n=");

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

        {

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

                {

                        printf("%d\t",c[i][j]);

                }

                printf("\n");

        }

        int k;

        printf("\n welcom to multiplication\n");

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

        {

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

                {

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

                        {

                                c[i][j]=c[i][j]+a[i][k]*b[k][j];

                        }

                }

 }

        printf("\n\nmatrix after multi is \n=");

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

        {

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

                {

                        printf("%d\t",c[i][j]);

                }

                printf("\n");

        }

}