nf = getnamed('R','frequency points');
nr = 45; # number of points for raius sweep

r = linspace(100,540,nr)*1e-9;
A = matrix(length(r),nf);

for (i = 1:length(r)){
switchtolayout;
setnamed('::model','r',r(i));
run;
T = -transmission('T');
R = transmission('R');
A(i,1:nf) = 1-R-T;
?"i = "+num2str(i);
}

f = getdata('T','f');
image(r(1:nr)*1e9,c/f*1e6,A,'radius (nm)','wavelength (um)');