A possible solution is as follows
clf;
figure
hold on
tic
% Using spherical coordinates to plot the cubes
for r = linspace(0,1,20) % define the range for radius of sphere
for phi = linspace(0,pi,10) % define the range for phi angle
for theta =linspace(0,2*pi,50) % define the range for theta angle
origin = [(r)*cos(theta)*sin(phi), (r)*sin(theta)*sin(phi), (r)*cos(phi)];
edges = [0.03,0.03,0.03];
plotcube(edges,origin,0.8,[0.5 0.5 0.5]);
end
end
end
toc
hold off
Note - Since it's a nested for loop, it is computationally expensive
