Mercurial > hg > octave-nkf

diff scripts/polynomial/ppval.m @ 6014:437f9086b967

Find changesets by keywords (author, files, the commit message), revision number or hash, or revset expression.
[project @ 2006-09-30 15:27:06 by jwe]
author jwe
date Sat, 30 Sep 2006 15:27:06 +0000
parents b0d4ff99a0c5
children 34f96dd5441b
line wrap: on
line diff
--- a/scripts/polynomial/ppval.m
+++ b/scripts/polynomial/ppval.m
@@ -41,13 +41,13 @@
     xi = xi(:);
     xn = length (xi);
     idx = lookup (pp.x(2:pp.n), xi) + 1;
-    dx = (xi - pp.x(idx))';
+    dx = (xi - pp.x(idx)).';
     dx = reshape (dx(ones(1,prod(pp.d)),:),[pp.d,xn]);
     c = reshape (pp.P(:,1), pp.n, prod (pp.d));
-    yi = reshape (c(idx,:)', [pp.d, xn]);
+    yi = reshape (c(idx,:).', [pp.d, xn]);
     for i  = 2 : pp.k;
       c = reshape (pp.P(:,i), pp.n, prod (pp.d));
-      yi = yi .* dx + reshape (c(idx,:)', [pp.d, xn]);
+      yi = yi .* dx + reshape (c(idx,:).', [pp.d, xn]);
     endfor
     if (transposed && isscalar (pp.d) && pp.d == 1)
       yi = yi.';