, где входные значения:
import numpy as np n = 16 delta = 0.1 W = -(1/delta)**2 t = 0 k1 = 10 k2 = -k1 sigma = 20*delta x0 = 10*delta x = np.linspace(k2,k1,n) PSI_0 = ((1/(np.pi*sigma**2)**0.25)*np.exp(-(x-x0)**2/(2*sigma**2)))