function HCPlattice
function HCPlattice
a=1;  %将a归一化。
x0=a/2*[0, 2, 1, -1, -2, -1, 1, 0, -1, 1, 0, 2, 1, -1, -2, -1, 1];
y0=sqrt(3)*a/6*[0, 0, 3, 3, 0, -3, -3, 2, -1, -1, 0, 0, 3, 3, 0, -3, -3];
z0=sqrt(6)*a/3*[0,0,0,0,0,0,0,1,1,1,2,2,2,2,2,2,2];
% x0=[0,a,a/2,-a/2,-a,-a/2,a/2,0,-a/2,a/2,0,a,a/2,-a/2,-a,-a/2,a/2];
%
% y0=[0,0,sqrt(3)*a/2,sqrt(3)*a/2,0,-sqrt(3)*a/2,-sqrt(3)*a/2,sqrt(3)*a/3,-sqrt(3)*a/6,-sqrt(3)*a/6,0,0,sqrt(3)*a/2,sqrt(3)*a/2,0,-sqrt(3)*a/2,-sqrt(3)*a/2];
%
% z0=[0,0,0,0,0,0,0,sqrt(6)*a/3,sqrt(6)*a/3,sqrt(6)*a/3,2*sqrt(6)*a/3,2*sqrt(6)*a/3,2*sqrt(6)*a/3,2*sqrt(6)*a/3,2*sqrt(6)*a/3,2*sqrt(6)*a/3,2*sqrt(6)*a/3];

%(注意，在matlab中运行时，要将x0,y0,z0后面的三个表达式拉到1排里。)


%下面这句可以用来初步查看坐标位置     

figure, plot3(x0,y0,z0,'o')


%下面画原子图

r=a/2;   %可自行调整半径  不同半径的画图效果差别蛮大的

rho=0:0.05*pi:2*pi;

phi=-pi:0.05*pi:pi;

figure

for i=1:17

    x=x0(i)+r*(cos(phi))'*(cos(rho));

    y=y0(i)+r*(cos(phi))'*(sin(rho));

    z=z0(i)+r*(sin(phi))'*(ones(size(rho)));

    surf (x,y,z);

    shading interp   %使图形光滑，无边界阴影

    hold on

end

axis equal