clear all
close all
n=input('¿CUANTAS CELDAS TIENE EL CUBO FCC?: ');
x=0:n; y=0:n; z=0:n;
%inicio cubo blanco
[X1,Y1,Z1]=meshgrid(x([1 end]),y,z);X1 = permute(X1,[2 1 3]);
Y1 = permute(Y1,[2 1 3]); Z1 = permute(Z1,[2 1 3]);
X1(end+1,:,:) = NaN; Y1(end+1,:,:) = NaN; Z1(end+1,:,:)=NaN;
[X2 Y2 Z2] = meshgrid(x,y([1 end]),z);
X2(end+1,:,:) = NaN; Y2(end+1,:,:) = NaN; Z2(end+1,:,:) = NaN;
[X3 Y3 Z3] = meshgrid(x,y,z([1 end]));
X3 = permute(X3,[3 1 2]); Y3 = permute(Y3,[3 1 2]); Z3=permute(Z3,[3 1 2]);
X3(end+1,:,:) = NaN; Y3(end+1,:,:) = NaN; Z3(end+1,:,:)=NaN;
h = line([X1(:);X2(:);X3(:)], [Y1(:);Y2(:);Y3(:)], [Z1(:);Z2(:);Z3(:)]);
set(h, 'Color',[1 1 1], 'LineWidth',1, 'LineStyle','-')
%fin cubo blanco
hold on
%inicio esfera
r=0.5*n*ones(50,50);
[th,phi]=meshgrid(linspace(0,2*pi,50),linspace(-pi,pi,50));
[i,j,k]=sph2cart(th,phi,r);
i=i+0.5*n;
j=j+0.5*n;
k=k+0.5*n;
surface(i,j,k,'FaceColor', 'none')
%inicio atomos
c=[n/2 n/2 n/2];
r=n/2;
for z1=0:n
for y1=0:n
for x1=0:n
if norm((c-[z1,y1,x1]))<=r
plot3(x1,y1,z1,'bo')
end
end
end
end
%
% %for para puntos medios
for c1=0:n-1
c1=c1+0.5;
for b1=0:n-1
b1=b1+0.5;
for a1=0:n
if norm((c-[c1,b1,a1]))<=r
plot3(a1,b1,c1,'ro')
end
end
end
end
%
for c2=0:n
for b2=0:n
for a2=0:n-1
a2=a2+0.5;
if norm((c-[c2,b2,a2]))<=r
plot3(a2,b2,c2,'go')
end
end
end
end
for c3=0:n
for b3=0:n
for a3=0:n-1
a3=a3+0.5;
if norm((c-[c3,b3,a3]))<=r
plot3(a3,b3,c3,'mo')
end
end
end
end
%fin atomos
hold off
axis off
view(3), axis vis3d
camproj perspective, rotate3d on