Ef = [0;0.1;0.2;0.5];
A = matrix(length(Ef),100);

for (i = 1:length(Ef)){

switchtolayout;
setmaterial("C (graphene) - broadband","chemical potential (eV)", Ef(i));
run;
A(i,1:100) = 1-transmission('R');

}

f = getdata('R','f');
plot(f/1e12,A(1,1:100),A(2,1:100),A(3,1:100),A(4,1:100),'frequency (THz)', 'Absorbance');
legend('0 eV','0.1 eV','0.2 eV','0.5 eV');