First commit calculate_optical_density.m