function [width, minAbsorption, detunings, absorptions] = loadSimulation(dataFile) data = load(dataFile); detunings = data(:,1); absorptions = data(:,2); [width, minAbsorption] = getPeak(detunings, absorptions); plot(detunings, absorptions); end