ammaML(:)) - log2(dgamma)).^2);
mseWP0 = mean((thetaWP_B(:) - log2(dgamma)).^2);
fprintf('\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',[mseAM0,mseAM,mseWP0,mseWP])

##### SOURCE END #####
--></body></html>