% generate comparametric plot from sqrtic data

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

grid
title("COMPARAMETRIC PLOT OF SQRTIC f(q) FOR RATIO k=2");
xlabel("f(q)")
ylabel("g(f)=f(kq)")
axis([0 1 0 1]);
gset size square
plot(f,g)


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)