0001 function [] = isotopomerViewer(mdv1, mdv2, names)
0002
0003
0004
0005 bins = round(sqrt(size(mdv1,2)));
0006
0007 for i = 1:size(mdv1,1)
0008 subplot(4,4,mod(i-1,16)+1);
0009 [x1,x2]=hist(mdv1(i,:), bins);
0010 plot(x2,x1);
0011 if max(x2) < .02
0012 plot(x2,x1, 'k');
0013 end
0014 hold on;
0015 [x1,x2]=hist(mdv2(i,:), bins);
0016 plot(x2,x1, 'g');
0017 if max(x2) < .02
0018 plot(x2,x1, 'k');
0019 end
0020 title(names{i});
0021 hold off;
0022 if mod(i,16) == 0
0023 pause;
0024 end
0025 end
0026
0027
0028
0029
0030