v0=vector( -1.435000,-1.900833,1.573333); v1=vector( -1.435000,2.469167,1.573333); v2=vector( -1.435000,-1.900833,-1.576667); v3=vector( 1.435000,-1.900833,1.573333); v4=vector( 1.435000,2.469167,1.573333); v5=vector( 1.435000,-1.900833,-1.576667); v6=vector( -1.435000,-0.250833,0.003333); v7=vector( -1.435000,2.469167,0.003333); v8=vector( 1.435000,-0.250833,0.003333); v9=vector( -1.435000,-0.250833,-1.576667); v10=vector( 1.435000,-0.250833,-1.576667); v11=vector( 1.435000,1.199167,0.003333); poly0=poly(list(v8,v6,v7,v11),FALSE); poly1=poly(list(v6,v8,v10,v9),FALSE); poly2=poly(list(v0,v1,v7,v6,v9,v2),FALSE); poly3=poly(list(v4,v3,v5,v10,v8,v11),FALSE); poly4=poly(list(v3,v0,v2,v5),FALSE); poly5=poly(list(v1,v4,v7),FALSE); poly6=poly(list(v3,v4,v1,v0),FALSE); poly7=poly(list(v2,v9,v10,v5),FALSE); poly8=poly(list(v4,v11,v7),FALSE); obj=mergepoly(list(poly0,poly1,poly2,poly3,poly4,poly5,poly6,poly7,poly8)); save("grnblk3.dat",convex(obj));