code: def calculate_vibr_energy(harmonic, anharmonic, i):
      return numpy.exp(-harmonic * i - anharmonic * (i**2))
  
label_desc: Use image node_id 1 for calling a global function with example usage: calculate_vibr_energy(harmonic, anharmonic, i) and returns: numpy