clf
// cyclic permutations of (0,pm 1, pm phi)
phi = (1+sqrt(5))/2;
param3d([1, -1],[0,0],[phi,phi]);
p=get("hdl"); p.foreground=5;
param3d([1, 0],[0,phi],[phi,1]);
p=get("hdl"); p.foreground=5;
param3d([1, phi],[0,1],[phi,0]);
p=get("hdl"); p.foreground=5;
param3d([1, phi],[0,-1],[phi,0]);
p=get("hdl"); p.foreground=5;
param3d([1, 0],[0,-phi],[phi,1]);
p=get("hdl"); p.foreground=5;
param3d([-1, 0],[0,phi],[phi,1]);
p=get("hdl"); p.foreground=5;
param3d([phi, 0],[1,phi],[0,1]);
p=get("hdl"); p.foreground=5;
param3d([phi, phi],[1,-1],[0,0]);
p=get("hdl"); p.foreground=5;
param3d([0, phi],[-phi,-1],[1,0]);
p=get("hdl"); p.foreground=5;
param3d([0, -1],[-phi,0],[1,phi]);
p=get("hdl"); p.foreground=5;

param3d([-1, 1],[0,0],[-phi,-phi]);
p=get("hdl"); p.foreground=5;
param3d([-1, 0],[0,-phi],[-phi,-1]);
p=get("hdl"); p.foreground=5;
param3d([-1, -phi],[0,-1],[-phi,0]);
p=get("hdl"); p.foreground=5;
param3d([-1, -phi],[0,1],[-phi,0]);
p=get("hdl"); p.foreground=5;
param3d([-1, 0],[0,phi],[-phi,-1]);
p=get("hdl"); p.foreground=5;
param3d([1, 0],[0,-phi],[-phi,-1]);
p=get("hdl"); p.foreground=5;
param3d([-phi, 0],[-1,-phi],[0,-1]);
p=get("hdl"); p.foreground=5;
param3d([-phi, -phi],[-1,1],[0,0]);
p=get("hdl"); p.foreground=5;
param3d([0, -phi],[phi,1],[-1,0]);
p=get("hdl"); p.foreground=5;
param3d([0, 1],[phi,0],[-1,-phi]);
p=get("hdl"); p.foreground=5;

param3d([-1, -phi],[0,-1],[phi,0]);
p=get("hdl"); p.foreground=5;
param3d([-1, -phi],[0,1],[phi,0]);
p=get("hdl"); p.foreground=5;
param3d([0, -phi],[phi,1],[1,0]);
p=get("hdl"); p.foreground=5;
param3d([0, 0],[phi,phi],[1,-1]);
p=get("hdl"); p.foreground=5;
param3d([phi, 0],[1,phi],[0,-1]);
p=get("hdl"); p.foreground=5;
param3d([phi, 1],[1,0],[0,-phi]);
p=get("hdl"); p.foreground=5;
param3d([phi, 1],[-1,0],[0,-phi]);
p=get("hdl"); p.foreground=5;
param3d([phi, 0],[-1,-phi],[0,-1]);
p=get("hdl"); p.foreground=5;
param3d([0, 0],[-phi,-phi],[1,-1]);
p=get("hdl"); p.foreground=5;
param3d([0, -phi],[-phi,-1],[1,0]);
p=get("hdl"); p.foreground=5;

title('$\mbox{Icosahedron } V=12, E=30, F=20, F_{\mbox{deg}} =3, V_{\mbox{deg}} = 5, \mbox{ Water}$','fontsize',4);
a=gca(); a.box="off"; a.isoview="on";
tmp=a.x_ticks; tmp.labels=[]; tmp.locations=[];
a.x_ticks=tmp; xlabel("");
a.y_ticks=tmp; ylabel("");
a.z_ticks=tmp; zlabel("");
a.axes_visible=["off","off","off"];

xs2pdf(0,'icosa.pdf');
