clear all

load('shifted_beta_vs_edge_angle.mat');

flake_num=size(xdata,1);

figure; hold on;
for i=1:flake_num
    plot(xdata{i},ydata{i},'.','MarkerSize',16);
end

box on
xlim([0 200])
ylim([-30 30])
xticks(0:50:200)
yticks(-30:10:30)
xlabel(['Edge Angle \beta (' char(176) ')']);
ylabel(['<\delta\beta> (' char(176) ')']);