TrigSimp is introduced

The case when alpha_i = 0 is treated. cos(th_i+th_j) and sine(th_i+th_j)
are replaced by CIJ and SIJ