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;