Mercurial > hg > octave-nkf

diff scripts/polynomial/residue.m @ 5158:955ed0dcdc5f

Find changesets by keywords (author, files, the commit message), revision number or hash, or revset expression.
[project @ 2005-02-23 00:40:57 by jwe]
author jwe
date Wed, 23 Feb 2005 00:40:57 +0000
parents c08cb1098afc
children 4c8a2e4e0717
line wrap: on
line diff
--- a/scripts/polynomial/residue.m
+++ b/scripts/polynomial/residue.m
@@ -294,7 +294,7 @@
       else
         cp = p (1:mpi-1);
       endif
-      rhs (1, rhi:rhi+lp-1) = prepad (poly (cp), lp);
+      rhs (1, rhi:rhi+lp-1) = prepad (poly (cp), lp, 0, 2);
       rhi = rhi + lp;
     endfor
     mpi = mpi + M (dpi);
@@ -302,11 +302,11 @@
   endwhile
   if (MM > 1)
     for index = 2:MM
-      lhs (index, :) = prepad (polyderiv (lhs (index-1, :)), lb);
+      lhs (index, :) = prepad (polyderiv (lhs (index-1, :)), lb, 0, 2);
       ind = 1;
       for rhi = 1:lp
         cp = rhs (index-1, ind:ind+lp-1);
-        rhs (index, ind:ind+lp-1) = prepad (polyderiv (cp), lp);
+        rhs (index, ind:ind+lp-1) = prepad (polyderiv (cp), lp, 0, 2);
         ind = ind + lp;
       endfor
     endfor