added bwalya's implementation as matlab scripts, working on comparison with f90 code