

# data is in: "fig_synthetic_comparagrams_data.txt"  
#fig_synthetic_comparagrams
#fig_synthetic_response_estimates
#fig_synthetic_response_estimates_reg

load  -force fig_synthetic_comparagrams_data.txt

xlabel('f=f(q)');
ylabel('g=f(kq)');
title(sprintf("Synthetic Comparagrams"));
grid
plot(f,g);
plot(
     f( 1:16:256),g(8, 1:16:256),"-;k=3.25;",
     f(11:16:256),g(7,11:16:256),"-;k=3.00;",
     f(14:16:256),g(6,14:16:256),"-;k=2.75;",
     f( 1:16:256),g(5, 1:16:256),"-;k=2.50;",
     f( 5:16:256),g(4, 5:16:256),"-;k=2.25;",
     f(10:16:256),g(3,10:16:256),"-;k=2.00;",
     f( 1:16:256),g(2, 1:16:256),"-;k=1.75;",
     f( 1:16:256),g(1, 1:16:256),"-@14;k=1.50;"
    )
gset key 0.85,0.43
axis

gset size square; gset terminal X11; gset terminal postscript portrait 22; system("rm synthetic_comparagrams.ps"); gset output "synthetic_comparagrams.ps"; gset key 0.78,0.51; replot; system("ps2epsi synthetic_comparagrams.ps"); system("gv synthetic_comparagrams.epsi");


