1.We used to trim nextPh by inv before computing reachable tube.
This could generate under-approximated results.
We change the order to compute tube before trimming
	modified:   ../../HybridAutomata/ha_stateReach.m
2. Consider the special case when ph_canon is empty for type 1.
	modified:   ../../Projectagon/Ph/ph_canon.m
Merge /ubc/cs/home/c/chaoyan/cra