Daniel@0: function subplot2(nrows, ncols, i, j)
Daniel@0: % function subplot2(nrows, ncols, i, j)
Daniel@0: 
Daniel@0: 
Daniel@0: sz = [nrows ncols];
Daniel@0: %k = sub2ind(sz, i, j)
Daniel@0: k = sub2ind(sz(end:-1:1), j, i);
Daniel@0: subplot(nrows, ncols, k);
Daniel@0: 
Daniel@0: if 0
Daniel@0:   ncols_plot = ceil(sqrt(Nplots));
Daniel@0:   nrows_plot = ceil(Nplots/ncols_plot);
Daniel@0:   Nplots = nrows_plot*ncols_plot;
Daniel@0:   for p=1:Nplots
Daniel@0:     subplot(nrows_plot, ncols_plot, p);
Daniel@0:   end
Daniel@0: end