c言語における行列成分の計算

前提・実現したいこと

c言語で以下条件で行列を表示させるプログラムを作成しています。

次元数をNとし、????(????,????) = (???? + 1) ∙ (???? + 1)

次元数N=5の場合には
1.0 2.0 3.0 4.0 5.0
2.0 4.0 6.0 8.0 10.0
3.0 6.0 9.0 12.0 15.0
4.0 8.0 12.0 16.0 20.0
5.0 10.0 15.0 20.0 25.0
という表示を得たいのです。

求めてる結果を得るにはどのようなコードに書き換えれば良いのでしょうか？

発生している問題・エラーメッセージ

以下のコードでは、得られる表示が欲しい結果と異なります。実際にN=5の場合に実行すると
1.0 2.0 3.0 4.0 5.0
0.0 2.0 4.0 6.0 8.0
10.0 0.0 3.0 6.0 9.0
12.0 15.0 0.0 4.0 8.0
12.0 16.0 20.0 0.0 5.0
と表示されます。

該当のソースコード

c
1#include <stdio.h>
2#include <stdlib.h>
3
4void PrintMatrix(double **A, int n){
5  int i, j;
6  double *ai;
7
8  ai = *A;
9
10  for( i = 0 ; i < n ; i++){
11    for( j = 0 ; j < n ; j++){
12      printf("%6.1f", *(ai+j));
13    }
14    ai+=n;
15    printf("\n");
16  }
17}
18
19int main(void){
20  int i = 0, j = 0, N;
21  double **A;
22
23  //次元Nをキーボードから入力
24  printf("N=");
25  scanf("%d", &N);
26
27  //行列A用のメモリを確保
28  A = malloc(N*sizeof(double *));
29  for(i=0;i<N;i++){
30    A[i] = malloc(N*sizeof(double));
31  }
32  if(A==NULL){
33    printf("Cannot allocate memory.\n");
34    exit(1);
35  }
36
37
38  //行列Aの要素
39  for( i = 0 ; i < N ; i++ ){
40    for( j = 0 ; j < N ; j++ ){
41      A[i][j] = (i+1)*(j+1);
42    }
43  }
44 
45  //行列Aを表示
46  printf("A = \n");
47  PrintMatrix(A,N);
48
49
50  //行列Aメモリを開放
51  free(A);
52  return 0;
53}