clear all
close all

load poly_square_T=200_K=50.mat

cop=calccop(reaction);

figure('Position',[50 50 500 600]);
plot(cop(:,1),cop(:,2),'.',[pos(:,1); pos(1,1)],[pos(:,2); pos(1,2)],'linewidth',2)

axis([-40 40 -60 60])
axis equal
grid on

xlabel('CoP x');
ylabel('CoP y');

legend('Trajectória medida','Trajectória prevista')