updated matlab code

discovered possible weakness in advection algo