ez:=proc():print(`G, F(a,b), Cube, Tet`):end:

G:=plot([[0,0],[1,1],[1,2],[1,3],[-1,2],[1,3],[2,3],[1,2],[-1,1],[1,2],[-1,2],[-1,1],[0,0],[0,-1],[-1,1],[0,0]],axes=none):

F:=proc(a,b): [4/(2*(1/a+1/b)-1)/a, 2/(2*(1/a+1/b)-1), 4/(2*(1/a+1/b)-1)/b]:end:

Cube:=plot([[0,0],[1,1],[3,1],[3,3],[4,4],[3,3],[1,3],[1,1],[0,0],[4,0],[3,1],[4,0],[4,4],[0,4],[1,3],[0,4],[0,0]],axes=none);

Tet:=plot([[0,0],[0,1],[-2,2],[0,1],[2,2],[0,1],[0,0],[2,2],[-2,2],[0,0]],axes=none);