Lennard-Jones potential simulation
问题 import pygame import random import numpy as np import matplotlib.pyplot as plt import math number_of_particles = 70 my_particles = [] background_colour = (255,255,255) width, height = 500, 500 sigma = 1 e = 1 dt = 0.1 v = 0 a = 0 r = 1 def r(p1,p2): dx = p1.x - p2.x dy = p1.y - p2.y angle = 0.5 * math.pi - math.atan2(dy, dx) dist = np.hypot(dx, dy) return dist def collide(p1, p2): dx = p1.x - p2.x dy = p1.y - p2.y dist = np.hypot(dx, dy) if dist < (p1.size + p2.size): tangent = math.atan2(dy,