Mackay-Icosahedron

This function generates so called Mackay-Icosahedron. In this function it's possible to change distance between atoms (bonding length) number of shells and also composition of two different kinds of atoms

Use following piece of code to visualise the particle:
radius = 1; details = 15;
for i=1:size(X,1)
[X_mesh,Y_mesh,Z_mesh] = ellipsoid(X(i),Y(i),Z(i),radius,radius,radius,details);
s(i) = surface(X_mesh,Y_mesh,Z_mesh);
set(s(i),'EdgeColor','none','Tag',strcat('atom_',int2str(i)));
if(type(i)==1)
set(s(i),'FaceColor','red');
else
set(s(i),'FaceColor','blue');
end
end

axis equal
axis off
light
lighting gouraud