PROBLEMA DE UM TRABALHO

[PROBLEMA DE UM TRABALHO]

escreva um programa que eprmita o usuario, informar os elementos de uma matriz A de numeros inteiros, ed ordem 3x3.
apos a digitação da matriz A o programa deve:
1. gerar a matriz B que corresponde e transposta de A.
2. gerar a matriz C que é a soma da A+B.
3. Gerar a matriz D tal que e a mult de A*B.
4. Mostrar as matrizes A,B,C e D.
5. Dizer quantos numeros primos existem na diagonal principal na matriz D
6. Mostrar o fatorial de todos os elementos da matriz B.
7. Mostrar a Media aritmetica dos elementos da diagonal principal da matriz C.

obs:. 
1.o programa tera as seguintes funcÕes
    a) uma que recebe o numero inteiro e devolve o seu fatorial.
    b) uma que recebe o numero inteiro e devolve 1 se ele for primo, 0 se não.
    c) a função main.
2. somente a função main faz I/O.

 

FIZ ATE AQUI

int main(int argc, char *argv[]) 
// MOSTRAR MATRIZ
    {int matA[3][3], i,j,k,matB[3][3],matC[j],matD[j],somaProd;
    for (i=0; i <3; i++)
        {for(j=0; j <3; j++)
            {printf("\nINFORME ELEMENTO (%d) (%d): ", i+1,j+1);
            scanf ("%d",&matA[j]);
            }
        }
// MATRIZ TRANSPOSTA DA MATRIZ A    
    printf("\n\tVEJA SUA MATRIZ\n");
    for(i=0;i<3;i++)
        {for(j=0;j<3;j++)
            printf("\t%d",matA[j]);
        printf("\n");
        }
    for(i=0; i< 3;i++)
        {for(j=0;j<3;j++)
        matB[j]=matA[j];
        }
    printf("\n\tVEJA MATRIZ B\n");
    for(i=0; i<3;i++)
        {for(j=0;j<3;j++)
            printf("\t%d", matB[j]);
        printf("\n");
        }
//SOMA DAS MATRIZES
    for(i=0;i <3; i++)
        {for(j=0;j<3;j++)
         matC[j] = matA[j] + matB[j];
        }
    printf("\n\tVEJA A SOMA DAS MATRIZES, A+B=C\n");
    for(i=0;i<3;i++)
        {for(j=0;j<3;j++)
            printf("\t%d",matC[j]);
        printf("\n");
    }
// MULT DAS MATRIZES
    for(i=0;i<3;i++)
        {for(j=0;j<3;j++)
            {somaProd = 0;
             for(k=0; k<3;k++)
                 somaProd = somaProd + matA[k] * matB[k][j];
            matD[j] = somaProd;
            }
        }
    printf("\n\tVEJA O PRODUTO A*B=D\n");
    for(i=0; i< 3; i++)
        {for(j=0;j<3;j++)
            printf("\t%d",matD[j]);
        printf("\n");
        }
//QTD DE PRIMOS NA DIAGONAL PRINCIPAL DA MATRIZ D