Mercurial > hg > octave-lyh

diff scripts/polynomial/poly.m @ 1336:33d5c2471c09

Find changesets by keywords (author, files, the commit message), revision number or hash, or revset expression.
[project @ 1995-08-25 01:53:42 by jwe]
author jwe
date Fri, 25 Aug 1995 01:53:42 +0000
parents 611d403c7f3d
children 5d29638dd524
line wrap: on
line diff
--- a/scripts/polynomial/poly.m
+++ b/scripts/polynomial/poly.m
@@ -45,7 +45,8 @@
     usage ("poly (x), where x is a vector or a square matrix");
   endif
   
-  y = [1, zeros (1, n)];
+  y = zeros (1, n+1);
+  y(1) = 1;
   for j = 1:n;
     y(2:(j+1)) = y(2:(j+1)) - v(j) .* y(1:j);
   endfor