Compare PI time and then use the fastest

- new PI using matricial evaluation as asked
- comparison with the one using Bellman
- the new one is faster, so I replaced PI by PIbis in mainTP1.m