Matlab implementation of Mie theory