code: def max_dim(elements, coordinates):
      
      atom_vdw_vertical = np.matrix(
          [[atomic_vdw_radius[i.upper()]] for i in elements])
      atom_vdw_horizontal = np.matrix(
          [atomic_vdw_radius[i.upper()] for i in elements])
      dist_matrix = euclidean_distances(coordinates, coordinates)
      vdw_matrix = atom_vdw_vertical + atom_vdw_horizontal
      re_dist_matrix = dist_matrix + vdw_matrix
      final_matrix = np.triu(re_dist_matrix)
      i1, i2 = np.unravel_index(final_matrix.argmax(), final_matrix.shape)
      maxdim = final_matrix[i1, i2]
      return i1, i2, maxdim
docs: Return the maximum diameter of a molecule.
  
      Parameters
      ----------
      elements : numpy.ndarray
          An array of all elements (type: str) in a molecule.
  
      coordinates : numpy.ndarray
          An array containing molecule's coordinates.
  
      Returns
      -------