adding MATLAB script
script updated constants