how to write a script about Hexagonal

I am a beginner from china.I want to write a script about   Hexagonal.But I
cannot now,can you help me?The condition is the same as the flowing script
except tha i want to change lattice  from the Triangle to Hexagonal

(set! num-bands 8)



(set! geometry-lattice (make lattice (size 1 1 no-size)

                         (basis1 (/ (sqrt 3) 2) 0.5)

                         (basis2 (/ (sqrt 3) 2) -0.5)))

(set! geometry (list (make cylinder

                       (center 0 0 0) (radius 0.2) (height infinity)

                       (material (make dielectric (epsilon 12))))))



(set! k-points (list (vector3 0 0 0)          ; Gamma

                     (vector3 0 0.5 0)        ; M

                     (vector3 (/ -3) (/ 3) 0) ; K

                     (vector3 0 0 0)))        ; Gamma

(set! k-points (interpolate 4 k-points))



(set! resolution 32)



(run-tm)

(run-te)
Thanks  very much

_______________________________________________
mpb-discuss mailing list
[email protected]
http://ab-initio.mit.edu/cgi-bin/mailman/listinfo/mpb-discuss