%Function for projecting red and magenta mobile stimulus

function redmagentasector2(sectors,time)

close all force

fig=figure;
axis off
axes('Position', [0 0 1 1])
xlim([-.1 .1])
ylim([-.1 .1])
set(fig,'Position',[0 0 1080 1080])

t=0;
% timeLimit=180;
timestep=0.03;
direction=1;
dirSwitchP=1;           %mean wait time between switches (s)
stimvel=0.05;
thetastep=2*pi()*stimvel*timestep;

theta=rand()*2*pi();
r=1;

while t<time
    cla;
    hold on;
    
    if  rand()<timestep/dirSwitchP;
        direction=direction*-1;
    end
    theta=theta+direction*thetastep;
    
    
    fill([1 1 -1 -1], [-1 1 1 -1], 'm','EdgeColor','none')
    
            X1=[];
        Y1=[];
    
    for i=1:sectors

        X1=[X1 0 cos(theta+2*pi()*2*i/(2*sectors)) cos(theta+2*pi()*(2*i+0.5)/(2*sectors)) cos(theta+2*pi()*(2*i+1)/(2*sectors))];
        Y1=[Y1 0 sin(theta+2*pi()*2*i/(2*sectors)) sin(theta+2*pi()*(2*i+0.5)/(2*sectors)) sin(theta+2*pi()*(2*i+1)/(2*sectors))];

    end
     fill(X1,Y1,'r','EdgeColor','none')
    
    
    pause(timestep)
    t=t+timestep;
end


fill([1 1 -1 -1], [-1 1 1 -1], 'k','EdgeColor','none')