Préciser des limites en abscisse et en ordonnée

xlim([0 Tsim]) % impose xmin=0 xmax=Tsim

ylim([0 40]) % impose ymin=0 ymax=40

OU

axis([0 Tsim 0 40]) % impose xmin=0 xmax=Tsim ymin=0 ymax=40