clc;
clear all;
close all;
figure(1)
Rad=1;
iter=0;
N=100;
fi=0:(2*pi/N):2*pi;
cx=Rad*cos(fi);
cy=Rad*sin(fi);
dy=(cy(2:N)-cy(1:N-1));
dx=(cx(2:N)-cx(1:N-1));
clLL=abs(2*pi*Rad-sum(sqrt(((dy).^2)+(dx).^2)))/(2*pi*Rad)*100;

display(clLL)
2*pi*Rad

id1=plot(N,clLL,'o-')
title('cdomileba')
ylabel('%')
xlabel('N')
N1=[];
SclLL=[];
%for N=logspace(log10(200),log10(10000),100);
for N=200:100:10000 %linspace(200,10000,100);
  N=floor(N)
N1=[N1,N];
%for fi=0:(2*pi/N):2*pi
%  iter=iter+1;
%  cx(iter)=Rad*cos(fi);
%  cy(iter)=Rad*sin(fi); log
%end
fi=0:(2*pi/N):2*pi;
%fi=space(0,2*pi,N);
cx=Rad*cos(fi);
cy=Rad*sin(fi);
dy=(cy(2:N)-cy(1:N-1));
dx=(cx(2:N)-cx(1:N-1));
clLL=abs(2*pi*Rad-sum(sqrt(((dy).^2)+(dx).^2)))/(2*pi*Rad)*100;
SclLL=[SclLL,clLL];
%display(clLL)
%2*pi*Rad

set(id1,'XData',N1,'YData',SclLL);
%pause(1.5);

end
figure(2)
plot(cx,cy,'o-')
axis('equal')
grid on
title('tsretsiri')
xlabel('x coordinati')
ylabel('y coordinati')