
paramS = {scalesS,Nf};

tic;
[aML, dgammaML, Sxest, evol_crit] = estim_altern(y,Dt,dgamma0,a0,paramWAV,paramWP,paramAM,paramS,stop_crit,Nit);
toc;
</pre><pre class="codeoutput"> Iteration 1 
 Relative update WP: Inf % 
 Relative update AM: 65.07 %

 Iteration 2 
 Relative update WP: 28.61 % 
 Relative update AM: 5.67 %

 Iteration 3 
 Relative update WP: 20.46 % 
 Relative update AM: 109.17 %

 Iteration 4 
 Relative update WP: 25.98 % 
 Relative update AM: 78.51 %

 Iteration 5 
 Relative update WP: 2.28 