Introduced beta parameter - ratio of immobile to mobile volumes
Refactored code, computations moved out from main
Fixed saving results to mat file