N=1000;
x=0:2/N:2;
t=2*abs(cos(pi*x))./sqrt(1.81-1.8*cos(pi*x));
plot(x,t);
N=1000;
x=0:2/N:2; %111111111
t=2*abs(cos(pi*x))./sqrt(1.81-1.8*cos(pi*x));
plot(x,t);
N=1000;
x=0:2/N:2;
t=2*abs(cos(pi*x))./sqrt(1.81-1.8*cos(pi*x));
plot(x,t);
N=1000;
x=0:2/N:2; %111111111
t=2*abs(cos(pi*x))./sqrt(1.81-1.8*cos(pi*x));
plot(x,t);