a=sqrt(2); % slope of g=af curve;
f=[1../a/a/a*f1cycle_sqrt_logsin(125) 1../a/a*f1cycle_sqrt_logsin(250) 1../a*f1cycle_sqrt_logsin(500) f1cycle_sqrt_logsin(1000) a*f1cycle_sqrt_logsin(2000) a*a*f1cycle_sqrt_logsin(4000)];
q=(125:7999)/1000;
grid
gset nokey
xlabel('q')
ylabel('f(q)')
%plot(q,f)
qintegers=[1/8 1/4 1/2 1 2 4 8];
fintegers=[1/a/a/a 1/a/a 1/a 1 a a*a a*a*a];
qtrue=(0:7999)/1000; % ok for linear plot, but problem for loglog plot
ftrue=sqrt(qtrue);
subplot(221);
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=sqrt(q);
Ftrue=log(ftrue)/log(2);
subplot(222);
plot(Q,F,Qintegers,Fintegers,'+',Qintegers,Fintegers,"*",Qtrue,Ftrue,".")
subplot(111);

