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'); % Symmetric y = -(1/(2*k^2) + R*cos(fi)); plot(x,y,'k'); end end axis([0, 2.2,-1.1, 1.1 ]);axis off;