Updated the Matalb code from default version

The main change includes:
(1) don't use global t, ci, and pho any more
(2) with the addition of expansion rate which can adjust the parameter
space