How can I create a periodic boundary that moves the integrin in my code to the opposite side of my domain when it crosses the boundary

N=10; % number of integrins
L=100; % size of the domain 100nm
D=.01; % diffusion coefficient
dt=1; % time step
F=randi([1, 10],1); % force on the bond with a random number of 1 to 10
P_a=randi([0,1],1); % binding rate of a random number from 0 to 1
P_ub=1/(1.5*F); % unbinding rate based upon our previous values
Tmax=500; % max time
% Intialize our integrins
x=randi([0,360],N,1); % x-coordinates generated randomily for each integrin
y=randi([0,50],N,1); % y-coordinates generated for each 
state=ones(N,1); % Each integrin is set to an inactive state at 0
% our figure
figure;
axis([0 L 0 L]);
set(gca,'nextplot','replacechildren');
% position update
for t = 0:dt:Tmax
% colors for our integrin states
colors = repmat([1 0 0], N,1); % red for inactive state
colors(state == 2,:,:) = repmat([0 1 0],sum(state == 2),1);
% plotting our integrins
 scatter(x,y,50,colors,'filled');
 title(sprintf('Time=%.2f',t));
 drawnow;
     for i=1:N
         w=randi([0,1],1); 
         if state(i)==1
         dx=randi([-360,360],1);
         dy=randi([-50,50],1);
         x(i)=x(i)+dx;
         y(i)=y(i)+dy; 
         end 
     end
     x(x<0)=0;
     x(x>L)=L;
     y(y<0)=0;
     y(y>L)=L;

     for i=1:N
         w_1=randi([0,1],1);
         if w_1 < P_a && state(i)==1
             state(i)=2;
             
         end
     end
     for i=1:N
         w1=randi([0,1],1);
         if state(i)==2 && w1<P_ub
            state(i)=1;
        
         end
     end
    
end

Answer 1

[Walter Roberson]

Replace

     x(x<0)=0;
     x(x>L)=L;
     y(y<0)=0;
     y(y>L)=L;

with

     x = mod(x, L);
     y = mod(y, L);

[BAILEY MCMASTER]

Awesome, I didn't knw anything about this mod application. Thank you so much.