--- /dev/null	Thu Jan 01 00:00:00 1970 +0000
+++ b/m-toolbox/classes/@rational/fromParfrac.m	Wed Nov 23 19:22:13 2011 +0100
@@ -0,0 +1,59 @@
+%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%
+%
+% FUNCTION:    fromParfrac
+%
+% DESCRIPTION: Construct a rational from parfrac model
+%
+% CALL:        r = Parfrac(a, pl)
+%
+% PARAMETER:   pl   - plist
+%
+%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%
+function r = fromParfrac(r, pli)
+
+  VERSION = '$Id: fromParfrac.m,v 1.10 2011/08/15 12:41:52 hewitson Exp $';
+
+
+  % get pzmodel info
+  ii = rational.getInfo('rational', 'From Parfrac');
+
+  % Set the method version string in the minfo object
+  ii.setMversion([VERSION '-->' ii.mversion]);
+
+  % Combine input plist with default values
+  pl = applyDefaults(ii.plists, pli);
+  
+  % Extrac model
+  
+  pf = find(pl, 'parfrac');
+  [a,b] = residue(pf.res,pf.poles,pf.dir);
+  
+  % Set fields
+  r.num = a;
+  r.den = b;
+
+  % Override model properties from the parfrac object
+  
+  if isempty(pl.find('ounits'))
+    pl.pset('ounits', pf.ounits);
+  end
+  
+  if isempty(pl.find('iunits'))
+    pl.pset('iunits', pf.iunits);
+  end
+    
+  if isempty(pl.find('name'))
+    pl.pset('name', sprintf('rational(%s)', pf.name));
+  end
+  
+  if isempty(pl.find('description'))
+    pl.pset('description', pf.description);
+  end
+  
+  % Add history
+  r.addHistory(ii, pl, [], pf.hist);
+
+  % Set object properties
+  r.setObjectProperties(pl);
+  
+end