c11 = [1,1,1]
c12 = [-1,1,1]
c13 = [-1,-1,1]
c14 = [1,-1,1]
c21 = [1,1,-1]
c22 = [-1,1,-1]
c23 = [-1,-1,-1]
c24 = [1,-1,-1]
cf1=[c11,c12,c13,c14]
cf2=[c21,c22,c23,c24]
show(sum([polygon(cf1, color="red"),
polygon(cf2),
line([c11,c21], color="black", thickness=3),
line([c12,c22]),
line([c13,c23]),
line([c14,c24]),
point([2,2,2], color="green"),
point([-2,-2,-2], color="green"),
point([0,0,0])
]),
frame=False, center=[0,0,0])