mseAM = mean(errAM.^2);
mseAM0 = mean((thetaAM_B(:) - a.^2).^2);
mseWP = mean((log2(dgammaML(:)) - log2(dgamma)).^2);
mseWP0 = mean((thetaWP_B(:) - log2(dgamma)).^2);
fprintf(<span class="string">'\nMean square errors: \n\n AM estimation:\n Baseline method: %.2e \n Algorithm estimation: %.2e \n\n Time Warping estimation:\n Baseline estimation: %.2e \n Algorithm estimation: %.2e\n'</span>,[mseAM0,mseAM,mseWP0,mseWP])
</pre><pre class="codeoutput">
Mean square errors: 

 AM estimation:
 Baseline method: 2.01e