Вычислить произведение положительных элементов матрицы, расположенных на главной диагонали и выше нее

p=1;for(i=0;i<N;i++)for(j=i;j<N;j++)if(A[i][j]>0)p*=A[i][j];

#include<stdio.h>
#include<stdlib.h>
#include<time.h>
int **matrix(int n,int m);
void print(int **ar,int n,int m);
void free_mem(int **ar,int n);
int main(void)
{
    int n,m,i,j,x=1;
    
    int **arr;
    srand(time(NULL));
    printf("Enter size of matrix: ");
    scanf("%d%d",&n,&m);
    arr=matrix(n,m);
    printf("Matrix:\n");
    print(arr,n,m);
    printf("Result: ");
    for(i=0;i<n;i++)
    for(j=i;j<m;j++)
    if(arr[i][j]>0)
    {
        printf("%d ",arr[i][j]);
        x*=arr[i][j];
    }
    printf("\n%d\n",x);
    free_mem(arr,n);
    return 0;
}
int **matrix(int n,int m)
{
    int i,j,**ar;
    ar=malloc(n*sizeof(int*));
    for(i=0;i<n;i++)
    ar[i]=malloc(m*sizeof(int));
    for(i=0;i<n;i++)
    for(j=0;j<m;j++)
    ar[i][j]=-50+rand()%100;
    
    return ar;
    
}
void print(int **ar,int n,int m)
{
    int i,j;
    for(i=0;i<n;i++)
    {
        for(j=0;j<m;j++)
        printf("%5d ",ar[i][j]);
        printf("\n");
    }
    printf("\n");
}
void free_mem(int **ar,int n)
{
    int i;
    for(i=0;i<n;i++)
    free(ar[i]);
    free(ar);
}