plot(TEC,FOM/max(FOM), 'bd')
xlabel('TEC17')
ylabel('Figure of Merit')

titl = ['Determination of Optimal TEC for ' filename];
title(titl)
set(gcf, 'Position', [50 25 400 350])