--- /dev/null	Thu Jan 01 00:00:00 1970 +0000
+++ b/m-toolbox/classes/+utils/@plottools/msuptitle.m	Wed Nov 23 19:22:13 2011 +0100
@@ -0,0 +1,129 @@
+% MSUPTITLE Puts a title above all subplots.
+%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%
+%
+% DESCRIPTION: MSUPTITLE Puts a title above all subplots.
+%	             MSUPTITLE('text') adds text to the top of the figure
+%            	 above all subplots (a "super title"). Use this function
+%              after all subplot commands.
+%
+% CALL:        xaxis(x1,x2)
+%
+% REMARK:      If the figure or axis units are non-default, this will break.
+%
+% VERSION:     $Id: msuptitle.m,v 1.1 2008/08/05 17:51:32 ingo Exp $
+%
+% HISTORY:     15-06-1995 Drea Thomas (drea@mathworks.com)
+%                 Creation.
+%              13-12-2000 John Cristion
+%                 Modified.
+%              13-03-2004 Mark Histed (histed@mit.edu)
+%                 Fix disappearing legend on last plot
+%              xx-xx-xxxx M Hewitson
+%                 Just a copy and rename with some small adjustements
+%
+%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%
+
+function hout = msuptitle(str)
+
+  % Parameters used to position the supertitle.
+
+  % Amount of the figure window devoted to subplots
+  plotregion = .92;
+
+  % Y position of title in normalized coordinates
+  titleypos  = .95;
+
+  % Fontsize for supertitle
+  %fs = get(gcf,'defaultaxesfontsize')+4;
+
+  fs = get(gcf,'defaultaxesfontsize');
+
+  % Fudge factor to adjust y spacing between subplots
+  fudge=1;
+
+  haold = gca;
+  figunits = get(gcf,'units');
+
+  % Get the (approximate) difference between full height (plot + title
+  % + xlabel) and bounding rectangle.
+
+  if (~strcmp(figunits,'pixels')),
+    set(gcf,'units','pixels');
+    pos = get(gcf,'position');
+    set(gcf,'units',figunits);
+  else
+    pos = get(gcf,'position');
+  end
+  ff = (fs-4)*1.27*5/pos(4)*fudge;
+
+  % The 5 here reflects about 3 characters of height below
+  % an axis and 2 above. 1.27 is pixels per point.
+
+  % Determine the bounding rectange for all the plots
+
+  % h = findobj('Type','axes');
+
+  % findobj is a 4.2 thing.. if you don't have 4.2 comment out
+  % the next line and uncomment the following block.
+
+  % h = findobj(gcf,'Type','axes');  % Change suggested by Stacy J. Hills
+
+  % If you don't have 4.2, use this code instead
+  ch = get(gcf,'children');
+  h=[];
+  for i=1:length(ch),
+    if strcmp(get(ch(i),'type'),'axes'),
+      h=[h,ch(i)];
+    end
+  end
+
+
+  max_y=0;
+  min_y=1;
+
+  oldtitle =0;
+  for i=1:length(h),
+    if (~strcmp(get(h(i),'Tag'),'suptitle')),
+      pos=get(h(i),'pos');
+      if (pos(2) < min_y), min_y=pos(2)-ff/5*3;end;
+      if (pos(4)+pos(2) > max_y), max_y=pos(4)+pos(2)+ff/5*2;end;
+    else
+      oldtitle = h(i);
+    end
+  end
+
+  if max_y > plotregion,
+    scale = (plotregion-min_y)/(max_y-min_y);
+    for i=1:length(h),
+      get(h(i),'Tag')
+      if ~strcmp(get(h(i),'Tag'),'legend')
+        pos = get(h(i),'position');
+        pos(2) = (pos(2)-min_y)*scale+min_y;
+        pos(4) = pos(4)*scale-(1-scale)*ff/5*3;
+        set(h(i),'position',pos);
+      end
+    end
+  end
+
+  np = get(gcf,'nextplot');
+  set(gcf,'nextplot','add');
+  if (oldtitle),
+    delete(oldtitle);
+  end
+  ha=axes('pos',[0 1 1 1],'visible','off','Tag','suptitle');
+  ht=text(.5,titleypos-1,str);set(ht,'horizontalalignment','center','fontsize',fs);
+  set(gcf,'nextplot',np);
+  axes(haold);
+
+  % fix legend if one exists
+  legH = legend;
+  if ~isempty(legH)
+    axes(legH);
+  end
+
+  if nargout,
+    hout=ht;
+  end
+
+end
+