r = [1;2];
r1 = [2;3];
maxfac = 1.15;

phi = [0 0.1] * pi;
phi = linspace(phi(1), phi(2), 10);
x = r*cos(phi);
y = r*sin(phi);

% xdiff = diff(r*cos(minmax(phi)));


phi1 = [max(phi) min(phi)] + pi;
phi1 = linspace(phi1(1), phi1(2), 10);
x1 = r1*cos(phi1) - 2*cos(phi1(1)) - 2*cos(phi1(1));
y1 = r1*sin(phi1) - 2*sin(phi1(1)) - 2*sin(phi1(1));

% phi2 = [mean(phi) min(phi)] + pi;
% phi2 = linspace(phi2(1), phi2(2), 5);
% x2 = r*cos(phi2) - 2 * r*cos(phi2(end));
% y2 = r*sin(phi2) - 2 * r*sin(phi2(end));

plot(x, y, x1, y1, 0, 0, 'o')
axis(1.1*repmat([-1 1]*max(abs([r; r1])),1,2))