new matlab script

as of May 2014