function f=f1cycle(n);

q=2aaaaaaaaaaaaaaa
f1=exp(sin(2*pi*(0:n-1)./n))+1;

f=f1;

aq=linspace(0,sqrt(2)-1,n); f=f+aq;

