optim.diff.norm function

performs the M step for measurement density parameters in multinormal case