fid=fopen('uvvis.txt');
k=1; while ~feof(fid) line{k}=fgetl(fid); k=k+1; end
fclose(fid);
line(1:2)=[]; data=[];
for k=1:length(line) data(k,:)=cell2mat(textscan(strrep(line{k},',','.'),'%f %f')); end
bg=load('bg.txt'); bg=bg(1:length(data),:); bg(:,2)=bg(:,2)-bg(end,2);
bg(:,2)=bg(:,2)*(data(data(:,1)==520,2)/bg(bg(:,1)==520,2));
data(:,2)=data(:,2)-bg(:,2);
plot(data(:,1),data(:,2),'k-');
hold on; plot(data(:,1),zeros(length(data(:,2)),1),'k--'); hold off;
xlabel('{\it wavelength} / nm'); ylabel('{\it intensity} / OD');
set(get(gca,'xlabel'),'fontsize',12);set(get(gca,'ylabel'),'fontsize',12);
set(get(gca,'xlabel'),'fontname','Arial');set(get(gca,'ylabel'),'fontname','Arial');
set(gca,'fontsize',12); set(gca,'fontname','Arial');
set(gcf,'paperunits','centimeters');set(gcf,'papersize',[16 10]);
set(gcf,'paperpositionmode','auto');
set(gca,'Units','centimeters');set(gca,'OuterPosition',[0 0 16 10]);
set(gcf,'Units','centimeters');oldpos = get(gcf,'Position');
set(gcf,'Position',[oldpos([1 2]) 16 10]);
print(gcf,'data.pdf','-dpdf');