change the runCalcDiffusionMatrix.m