added Matlab code and results