function y = arc_function(x,A,L)
%
piL=pi/L;
c=(A*piL)^2;
y = sqrt(1.+ c*(cos(piL*x).^2) );