q=10*(0:999)/1000;
f=q.^2;
subplot(221); xlabel('q'); ylabel('f'); plot(q,f);
g=4*f;
subplot(222); xlabel('f'); ylabel('g'); plot(f,g);
f=(q+sin(2*pi*log2(q))).^2;
f(1)=0;
g=4*f;
subplot(223); xlabel('q'); ylabel('f'); plot(q,f);
subplot(224); xlabel('f'); ylabel('g'); plot(f,g);

