January 3, 2014 3:44 pm / Leave a comment
clc; close all; fi = 0:0.01:2*pi;
figure('Position',[0 0 500 500])
hold on;
for h =1:15
for k = 1:15
R = 1/(2*k^2);
x = h/k + R*sin(fi);
y = 1/(2*k^2) + R*cos(fi);
plot(x,y,'k');
y = -(1/(2*k^2) + R*cos(fi));
plot(x,y,'k');
end
end
axis([0, 2.2,-1.1, 1.1 ]);axis off;

Like this:
Like Loading...