code: def pairwise(X):
      M = len(X)
      N = len(X[0])
      D = [[0 for x in range(M)] for y in range(M)]
      "omp parallel for"
      for i in range(M):
          for j in range(M):
              d = 0.0
              for k in range(N):
                  tmp = X[i][k] - X[j][k]
                  d += tmp * tmp
              D[i][j] = math.sqrt(d)
      return D
  
label_desc: Use image node_id 1 for calling a global function with example usage: pairwise(X) and returns: D