https://www.dropbox.com/scl/fi/szbhz372p8v98j5csblfy/output.pdf?rlkey=9v5vk5unuk6xmfqtyjz2znlf1&dl=0
import matplotlib.pyplot as plt
import numpy as np
def ejes(max):
hx=np.array([-max,max])
hy=np.array([0,0])
plt.plot(hx,hy,'-k',linewidth=2)
vx=np.array([0,0])
vy=np.array([-max,max])
plt.plot(vx,vy,'-k',linewidth=2)
def hexagono(cx,cy):
x1 = np.array([cx])
y1 = np.array([cy])
plt.plot(x1,y1,'ob')
x5=np.array([cx-10,cx-5 ,cx+5 ,cx+10,cx+5 ,cx-5 ,cx-10])
y5=np.array([cy+0 ,cy+8.6,cy+8.6,cy+0 ,cy-8.6,cy-8.6,cy+0])
plt.plot(x5,y5,'--ob')
ejes(40)
hexagono(0,0)
a_x = 15
a_y = 8.6
b_x = 15
b_y = -8.6
hexagono(a_x,a_y)
hexagono(b_x,b_y)
hexagono(-a_x,-a_y)
hexagono(-b_x,-b_y)
hexagono(2*a_x,2*a_y)
hexagono(2*b_x,2*b_y)
hexagono(-2*a_x,-2*a_y)
hexagono(-2*b_x,-2*b_y)
hexagono(a_x-b_x,a_y-b_y)
plt.grid()
plt.axis('square')
plt.show()