function []=figure1b()
close all
figure()
hold on
%% Loading config data
data_name=sprintf('Conf.dat');
fid = fopen(data_name,'r'); numbers=fscanf(fid, '%f'); fclose('all');
N=numbers(1); % number of particles
L=numbers(2); % linear size of the system
numbers=numbers(3:end);
numbers=reshape(numbers,5,N);
xs=(numbers(1,:)).'; % x-positions
ys=(numbers(2,:)).'; % y-positions
sigmas=(numbers(5,:)).'; % diameters

%% Scatter disks
xylim=0.2*L; % plot range
for ii=1:N
    x=xs(ii);
    y=ys(ii);
    if x<xylim+1 && y<xylim+1
        sigma=sigmas(ii)/2;
        viscircles([x,y],sigma,'Color',[0,0,0],'EnhanceVisibility',false);
    end
end
xlim([0,xylim])
ylim([0,xylim])
axis off
set(gca,'Xtick',[],'Ytick',[])
axis square
print(gcf, '-dpdf', 'figure1b.pdf');
end