f=[1/64*f1cycle_sqr(125) 1/16*f1cycle_sqr(250) 1/4*f1cycle_sqr(500) f1cycle_sqr(1000) 4*f1cycle_sqr(2000) 16*f1cycle_sqr(4000)];
q=(125:7999)/1000;
grid
gset nokey
xlabel('q')
ylabel('f(q)')
qintegers=[1/8 1/4 1/2 1 2 4 8];
fintegers=[1/64 1/16 1/4 1 4 16 64];
plot(q,f,qintegers,fintegers,'+',qintegers,fintegers,"*")

qtrue=(0:7999)/1000; % ok for linear plot, but problem for loglog plot
ftrue=qtrue.^2;

subplot(121);
plot(q,f,qintegers,fintegers,'+',qintegers,fintegers,"*",qtrue,ftrue,".")
Q=log(q)/log(2);
F=log(f)/log(2);
Qintegers=log(qintegers)/log(2);
Fintegers=log(fintegers)/log(2);
Qtrue=log(q)/log(2); % log plot problem because it goes to zero
ftrue=q.^2;
Ftrue=log(ftrue)/log(2);
subplot(122);
xlabel('Q')
ylabel('F(Q)')
title('Comparametric Periodicity (''Comperiodicity'')')
plot(Q,F,Qintegers,Fintegers,'+',Qintegers,Fintegers,"*",Qtrue,Ftrue,".")
subplot(111);

