code: def var_str(name, shape):
      
  
      size = prod(shape)
      ind = (indices(shape) + 1).reshape(-1, size)
      names = ['[' + ','.join(list(map(str, i))) + ']' for i in zip(*ind)]
      names[0] = '%s %s' % (name, names[0])
      return names
docs: Return a sequence of strings naming the element of the tallyable object.
  
      :Example:
      >>> var_str('theta', (4,))
      ['theta[1]', 'theta[2]', 'theta[3]', 'theta[4]']