function errorxy(Data,varargin)
% Plot graphs with error bars in both axes. 
% Package: plot
% Description: Plot graphs with error bars in both axes. 
% Input  : - Data matrix.
%          * an arbitrary pairs of keyword and value parameters, where
%            keyword is one of the followings:
%            (In case that the second parameter is a number then call the
%             old version errorxy_o).
%            'ColX'     - Column number containing X values, default is 1.
%            'ColY'     - Column number containing Y values, default is 2.
%            'ColXe'    - Column number containing error X values.
%                         If one parameter is given then assume upper
%                         and lower errors are the same, if two values
%                         are given then the first is the left and the
%                         second is the right errors,
%                         default is NaN (no X errors).
%            'ColYe'    - Column number containing error Y values.
%                         If one parameter is given then assume right
%                         and left errors are the same, if two values
%                         are given then the first is the lower and the
%                         second is the upper errors,
%                         default is 3.
%                         If one of the two sided error is NaN, then plot
%                         an upper/lower limit arrow.
%            'Marker'   - Marker type, see plot for options, default is 'o'.
%            'MarkSize' - Marker size, default is 6.
%            'EdgeColor'- Marker edge color, default is 'b'.
%            'FaceColor'- Marker face color, default is NaN.
%            'ColorEB'  - Error bars color, default is the same as 'EdgeColor'.
%            'LineEB'   - Error bars line type, default is '-'.
%            'WidthEB'  - Error bars line width, default is 0.5.
%            'EdgeEB'   - Length [ticks] of error bar edge, default is 0.
%            'Con'      - Connect data points {'y' | 'n'}, default is 'n'.
%            'ConColor' - Color of Connecting line, default is the
%                         same as 'EdgeColor'
%            'ConLine'  - Type of connecting line, default is '-'.
%            'ConWidth' - Width of connecting line, default is 0.5.
%            'Hist'     - Plot histogram {'y' | 'n'}, default is 'n'.
%            'HistFaceColor'- Histogram face color, default is the same
%                             as 'EdgeColor'.
%            'HistEdgeColor'- Histogram edge color, default is the same
%                             as 'EdgeColor'.
%            'ScaleX'   - X scale {'linear' | 'log'}, default is 'linear'.
%            'ScaleY'   - Y scale {'linear' | 'log'}, default is 'linear'.
%            'DirX'     - X direction {'normal' | 'reverse'},
%                         default is 'normal'.
%            'DirY'     - Y direction {'normal' | 'reverse'},
%                         default is 'normal'.
% Output : null
% Plot   : errorxy plot
% See also: errorxy_o.m
% Tested : Matlab 5.3
%     By : Eran O. Ofek                    Mar 2004
%    URL : http://weizmann.ac.il/home/eofek/matlab/
% Example: errorxy(Data,'ColYe',[5 6],'Marker','^');
% Reliable: 2
%----------------------------------------------------------------------------
Ndots = 450;    % approximate number of dots in paper

if (nargin>1 && ischar(varargin{1})==0)
   %--- call old erroxy_o ---
   plot.errorxy_o(Data,varargin{:});
else
   %--- new version of errorxy ---

   Nvarg = length(varargin);
   if (0.5.*Nvarg~=floor(0.5.*Nvarg))
      error('Illegal number of input arguments');
   end

   %--- set default values ---
   ArrowLength  = 0.2;
   ArrowBase    = 0.5.*ArrowLength;

   ColX         = 1;
   ColY         = 2;
   ColXe        = [NaN, NaN];
   ColYe        = [3, 3];
   Marker       = 'o';
   MarkSize     = 6;
   EdgeColor    = 'b';
   FaceColor    = 'n';
   ColorEB      = NaN;
   LineEB       = '-';
   WidthEB      = 0.5;
   EdgeEB       = 0;
   Con          = 'n';
   ConColor     = NaN;
   ConLine      = '-';
   ConWidth     = 0.5;
   Hist         = 'n';
   HistFaceColor= NaN;
   HistEdgeColor= NaN;
   ScaleX       = 'linear';
   ScaleY       = 'linear';
   DirX         = 'normal';
   DirY         = 'normal';
   
   for I=1:2:Nvarg-1
      switch varargin{I}
       case 'ColX'
          ColX        = varargin{I+1};
       case 'ColY'
          ColY        = varargin{I+1};
       case 'ColXe'
          ColXe       = varargin{I+1};
       case 'ColYe'
          ColYe       = varargin{I+1};
       case 'Marker'
          Marker      = varargin{I+1};
       case 'MarkSize'
          MarkSize    = varargin{I+1};
       case 'EdgeColor'
          EdgeColor   = varargin{I+1};
       case 'FaceColor'
          FaceColor   = varargin{I+1};
       case 'ColorEB'
          ColorEB     = varargin{I+1};
       case 'LineEB'
          LineEB      = varargin{I+1};
       case 'WidthEB'
          WidthEB     = varargin{I+1};
       case 'EdgeEB'
          EdgeEB      = varargin{I+1};
       case 'Con'
          Con         = varargin{I+1};
       case 'ConColor'
          ConColor    = varargin{I+1};
       case 'ConLine'
          ConLine    = varargin{I+1};
       case 'ConWidth'
          ConWidth   = varargin{I+1};
       case 'Hist'
          Hist       = varargin{I+1};
       case 'HistFaceColor'
          HistFaceColor  = varargin{I+1};
       case 'HistEdgeColor'
          HistEdgeColor  = varargin{I+1};
       case 'ScaleX'
          ScaleX     = varargin{I+1};
       case 'ScaleY'
          ScaleY     = varargin{I+1};
       case 'DirX'
          DirX       = varargin{I+1};
       case 'DirY'
          DirY       = varargin{I+1};
       otherwise
          error('Unknown keyword option');
      end
   end
   
   %--- set NaN colors ---
   if (isnan(ColorEB)==1)
      ColorEB = EdgeColor;
   end
   if (isnan(ConColor)==1)
      ConColor = EdgeColor;
   end
   if (isnan(HistFaceColor)==1)
      HistFaceColor = EdgeColor;
   end
   if (isnan(HistEdgeColor)==1)
      HistEdgeColor = EdgeColor;
   end
   
   if (length(ColXe)==1)
      ColXe = [ColXe, ColXe];
   end
   if (length(ColYe)==1)
      ColYe = [ColYe, ColYe];
   end
   
   %---------------------
   %--- Hold handling ---
   %---------------------
   NextPlot = get(gcf,'NextPlot');
   hold on;
   box on;
   
   %-----------------------------------
   %--- Data span for edge plotting ---
   %-----------------------------------
   SpanX = max(Data(:,ColX)) - min(Data(:,ColX));
   SpanY = max(Data(:,ColY)) - min(Data(:,ColY));
   EdgeX = SpanX.*EdgeEB./Ndots;
   EdgeY = SpanY.*EdgeEB./Ndots;
   
   
   %----------------------
   %--- plot histogram ---
   %----------------------
   switch Hist
    case 'n'
       % do nothing
    case 'y'
       Hhist = bar(Data(:,ColX),Data(:,ColY));
       set(Hhist,'FaceColor',HistFaceColor,...
                 'EdgeColor',HistEdgeColor);
    otherwise
       error('Unknown Hist option');
   end
   
   %------------------------
   %--- plot data points ---
   %------------------------
   N = size(Data,1);
   
   for I=1:1:N
      H         = plot(Data(I,ColX),Data(I,ColY),'o');
      Hpoint(I) = H;
      set(Hpoint,'Marker',Marker,...
            'MarkerSize',MarkSize,...
            'MarkerEdgeColor',EdgeColor,...
            'MarkerFaceColor',FaceColor);
      %--- plot Y error bars ---
      if (xor(isnan(Data(I,ColYe(1))),isnan(Data(I,ColYe(2))))==1)
         %--- plot lower/upper limit arrow ---
	 if (isnan(Data(I,ColYe(1)))==1)
	    ArrowStart = Data(I,ColY)+Data(I,ColYe(2));
            ArrowEnd   = Data(I,ColY);
         else
	    ArrowStart = Data(I,ColY)-Data(I,ColYe(1));
            ArrowEnd   = Data(I,ColY);
         end
%[	    ArrowStart, ArrowEnd]
	 [HarrowLine(I), HarrowHead(I)] = arrow2([Data(I,ColX);Data(I,ColX)],[ArrowStart;ArrowEnd],...
                                                  ArrowLength.*abs(ArrowStart-ArrowEnd),...
						  ArrowBase.*abs(ArrowStart-ArrowEnd));
         set(Hyerr,'Color',ColorEB,...
                   'LineStyle',LineEB,...
      	           'LineWidth',WidthEB);
      else
         %--- plot normal error bars ---

         if (isnan(ColYe(1))==0)
            H        = plot([Data(I,ColX); Data(I,ColX)],...
                            [Data(I,ColY)-Data(I,ColYe(1)); Data(I,ColY)+Data(I,ColYe(2))],...
                            '-');    
            Hyerr(I) = H;
            set(Hyerr,'Color',ColorEB,...
        	        'LineStyle',LineEB,...
      	        'LineWidth',WidthEB);
      
            %--- Error bar edge ---
            if (EdgeEB>0)
               %--- Y error bar lower edge ---
               H        = plot([Data(I,ColX)-EdgeX; Data(I,ColX)+EdgeX],...
                               [Data(I,ColY)-Data(I,ColYe(1)); Data(I,ColY)-Data(I,ColYe(1))],...
   		               '-');
               Heyel(I) = H;
      	       set(Heyel,'Color',ColorEB,...
                         'LineStyle',LineEB,...
      	                 'LineWidth',WidthEB);
               %--- Y error bar upper edge ---
       	       H        = plot([Data(I,ColX)-EdgeX; Data(I,ColX)+EdgeX],...
                               [Data(I,ColY)+Data(I,ColYe(2)); Data(I,ColY)+Data(I,ColYe(2))],...
      		         '-');
               Heyeu(I) = H;
      	       set(Heyeu,'Color',ColorEB,...
        	           'LineStyle',LineEB,...
      	           'LineWidth',WidthEB);
            end
         end
      end
      %--- plot X error bars ---
      if (isnan(ColXe(1))==0)
         H        = plot([Data(I,ColX)-Data(I,ColXe(1)); Data(I,ColX)+Data(I,ColXe(2))],...
                         [Data(I,ColY); Data(I,ColY)],...
                         '-');    
         Hxerr(I) = H;
         set(Hxerr,'Color',ColorEB,...
     	        'LineStyle',LineEB,...
   	        'LineWidth',WidthEB);
   
         %--- Error bar edge ---
         if (EdgeEB>0)
            %--- X error bar lower edge ---
    	 H        = plot([Data(I,ColX)-Data(I,ColXe(1)); Data(I,ColX)-Data(I,ColXe(1))],...
                            [Data(I,ColY)-EdgeY; Data(I,ColY)+EdgeY],...
   		         '-');
            Hexel(I) = H;
   	 set(Hexel,'Color',ColorEB,...
     	           'LineStyle',LineEB,...
   	           'LineWidth',WidthEB);
            %--- X error bar upper edge ---
    	 H        = plot([Data(I,ColX)+Data(I,ColXe(2)); Data(I,ColX)+Data(I,ColXe(2))],...
                            [Data(I,ColY)-EdgeY; Data(I,ColY)+EdgeY],...
   		         '-');
            Hexer(I) = H;
   	 set(Hexer,'Color',ColorEB,...
     	           'LineStyle',LineEB,...
   	           'LineWidth',WidthEB);
         end
      end
   
   
   end
   
   %---------------------------
   %--- connect data points ---
   %---------------------------
   switch Con
    case 'n'
       % do nothing
    case 'y'
       Hcon = plot(Data(:,ColX),Data(:,ColY),'-');
       set(Hcon,'Color',ConColor,...
    	     'LineStyle',ConLine,...
      	     'LineWidth',ConWidth);
    otherwise
       error('Unknown Con option');
   end
   
   %---------------------------
   %--- set scale/direction ---
   %---------------------------
   set(gca,'XScale',ScaleX,...
           'YScale',ScaleY,...
           'XDir',DirX,...
           'YDir',DirY);
   
   %--- return hold to riginal state ---
   set(gcf,'NextPlot',NextPlot);
end
