import numpy as np
import matplotlib.pyplot as plt
x = np.arange(-10, 10, 0.1) # x軸
y = np.arange(-10, 10, 0.1) # y軸
X, Y = np.meshgrid(x, y)
Z = 0
for Y_off in np.arange(-7.5,7.5,0.1):
Z = Z + 1/np.sqrt((X-5)**2 + (Y-Y_off)**2+1) - 1/np.sqrt((X+5)**2 + (Y-Y_off)**2+1)
#∞を回避するのに分母に1を足した
cont = plt.contour(X,Y,Z,colors=['r', 'g', 'b'])
cont.clabel(fmt='%1.1f', fontsize=14)
plt.xlabel('X', fontsize=14)
plt.ylabel('Y', fontsize=14)
plt.show()
©K.Tachibana