Я получил ошибку деления на ноль во время выполнения при использовании al oop, и я не знаю почему. Не волнуйтесь, код еще не закончен. Я пытаюсь преобразовать массив с помощью функции ulj, если вы можете помочь. Вот код:
import numpy as np
def ulj(x,epsilon = 4*(10**-4),sigma=6):
"""Compute Lennard-Jones potential as a function of inter-atomic distance"""
return epsilon*((((sigma)/(x))**12)-((((sigma)/(x))**6)))
positions=np.array([0., 7.0,5,6,2,6,2,2,4,3]) # Initial position of the two objects (in atomic units)
i = 0
x1 = []
x2 = []
while i <(len(positions)):
x1.append(ulj(positions[i]))
print(x1)
Спасибо за вашу помощь