lot(t,thetaWP_B,<span class="string">'g'</span>,t,log2(dgammaML),<span class="string">'r'</span>,t,log2(dgamma),<span class="string">'b-.'</span>,th,WPinf,<span class="string">'b--'</span>,th,WPsup,<span class="string">'b--'</span>,<span class="string">'linewidth'</span>,2); grid <span class="string">on</span>; axis <span class="string">tight</span>;
<span class="comment">% title('Warping estimation');</span>
legend({<span class="string">'Baseline estimate'</span>,<span class="string">'JEFAS estimate'</span