range = linspace(0,2*pi,3000);

Rx = ones(6, 3000);
Ry = ones(6, 3000);

for i = 1:6
    for n = 1:3000
    point = position(range(n));
    Rx(i,n) = point(i,1)*cos(point(i,2));
    Ry(i,n) = point(i,1)*sin(point(i,2));
    end
end

F(3000) = struct('cdata',[],'colormap',[]);

for i = 1:3000
    plot(   [0 Ry(1,i)] , [0 Rx(1,i)], '-o', ...
            [Ry(1,i) Ry(1,i)+Ry(4,i)], [Rx(1,i) Rx(1,i)+Rx(4,i)], '-o', ...
            [0 Ry(2,i)], [0 Rx(2,i)], '-o', ...
            [0 Ry(6,i)], [0 Rx(6,i)], '-o',...
            [Ry(6,i) Ry(6,i)+Ry(5,i)], [Rx(6,i) Rx(6,i)+Rx(5,i)], '-o', 'LineWidth', 1)

    set(gca, 'XDir','reverse')
    ylim([-6 6]);
    xlim([-6 6]);
    grid
    F(i) = getframe(gcf);
end