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