%  Convert to English units 
function [varargout] = EU(varargin)
  J=length(varargin);
  units1=varargin{J};
  for j=1:J-1
    val1=varargin{j};
    switch units1 % check for simple temperature conversion
    case 'K'
      varargout{j}=val1*9/5;
      varargout{J}='R';
      continue;
    case 'C'
      varargout{j}=(val1+273.15)*9/5;
      varargout{J}='R';
      continue;
    case 'F'
      varargout{j}=val1+459.67;
      varargout{J}='R';
      continue;
    end

    SI2EU={'rad',180/pi,'deg';
           'J',778.169297,'ft*lbf';
           'N',1/4.4482216,'lbf';
           'm',1/0.3048,'ft';
           'kg',1/0.45359237,'lbm';
           'Pa',1/47.880259,'lbf/ft^2';
           'W',0.73756215,'lbf*ft/s';
           'K',9/5,'F';
           'deg',1,'deg';
           'lbf',1,'lbf';
           'ft',1,'ft';
           'lbm',1,'lbm';
           's',1,'s'};
    [varargout{j} varargout{J}] = ConvertUnits(SI2EU,val1,units1);
  end
end
