% octavefile to load response function from picture v115 of test chart % generate a plot figure % must load one of frame 100 to 118 [f115 y115]=loadf(115); % chosen one % 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 "AVERAGE OF TEST PATTERN TRANSMISSION for v115" %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,y115)