% generate comparametric plot from quadratic data

%q=linspace(0,0.893000); % q on the DSC testpattern goes up to .893
q=linspace(0,0.893000/2);
f=(quadraticf(q));
g=(quadraticf(2*q));

system("rm bacon.ps") % remove else it appends???
gset term postscript portrait
gset output "bacon.ps"
grid
title("COMPARAPLOT OF QUADRATIC f(q) FOR RATIO k=2");
xlabel("f(q)")
ylabel("g(f(q))=f(kq)")
axis([0 1 0 1]);
gset size square
plot(f,g)