n</span>
[fn,gn,wn] = fh(mn);

<span class="comment">% write results</span>
rsf_write_all(<span class="string">'wn.rsf'</span>,{},wn,[dd(3) 0],[od(3) 0]);
rsf_write_all(<span class="string">'mn.rsf'</span>,{},reshape(mn,n),d,o);

<span class="comment">% end</span>
cd(curdir);
</pre><p class="footer"><br><br></p>
      </div>
      <!--
##### SOURCE BEGIN #####
%%
% The |mbase|, |mlsls|, |mstls| and |mstst| scripts reproduce the results
% from A.Y. Aravkin, T. van Leeuwen and F.J. Herrmann - Source estimatio