% octavefile to load response functions from pictures of test chart and
% generates a plot figure

% must load one of frame 100 to 118
[f112 y112]=loadf(112); % one stop under
[f114 y114]=loadf(114); % one third stop under
[f115 y115]=loadf(115); % chosen one
[f116 y116]=loadf(116); % one third stop over
[f118 y118]=loadf(118); % one stop over
% true black then 11 values of steps:
q=[.011; .021; .048; .079; .126; .185; .264; .358; .472; .622; .733; .893];

% now make a nice plot before we lose "y" and return only with "f"

gset term postscript portrait
gset output "bacon.ps"
gset title "AVERAGES OF TEST PATTERN TRANSMISSIONS for v112,4,5,6,8"
%gset size square
gset size 1,0.5

x=linspace(0,11,length(y115));
grid
axis([0 12 0 256]);
xlabel("spatial coordinate, x")
ylabel("response function, y(x)=f(q(x))")
%plot(x,y114,x,y115,x,y116)
plot(x,[y112; y114; y115; y116; y118])

