clear all
%oloklirosis synartisis
N=2*255;
e=0.1;
E=10.^[-1:-1:-6];
e=10^(-6);
%E=0.01*[0.1,0.05,0.01,0.005,0.001];
ALPHA=linspace(pi/1000,pi-pi/100,100);

for ia=1:length(ALPHA);
    alpha=ALPHA(ia);
    %e=E(ia);
a=-1+e;
b=1-e;

%x=linspace(a,b,N)';
x=[a:e:b]';N=length(x);
delta=e;
%delta=x(2)-x(1);

f=x.^2;



c0=cos(alpha);
c=cos(alpha*x);
f=1./sqrt(c-c0);

%plot(x,f)


D=delta*ones(1,N);
D(1)=D(1)/2;D(N)=D(N)/2;
%I/delta

I=D*f;

T(ia)=sqrt(2)*alpha*I;

end;
figure(2)
plot(ALPHA,T,'b',ALPHA,2*pi*exp(0*ALPHA),'r')
hold on
[T(ia),2*pi]


%plot(-log10(E),I)