Mercurial > hg > octave-nkf

comparison scripts/polynomial/deconv.m @ 3754:e63a3a6d7797

Find changesets by keywords (author, files, the commit message), revision number or hash, or revset expression.
[project @ 2000-12-14 21:57:14 by jwe]
author jwe
date Thu, 14 Dec 2000 21:57:15 +0000
parents e031284eea27
children 22bd65326ec1
comparison
equal deleted inserted replaced
3753:f751e43de300 3754:e63a3a6d7797
55 b = filter (y, a, 1); 55 b = filter (y, a, 1);
56 else 56 else
57 b = 0; 57 b = 0;
58 endif 58 endif
59 59
60 b = polyreduce (b);
61
62 lc = la + length (b) - 1; 60 lc = la + length (b) - 1;
63 if (ly == lc) 61 if (ly == lc)
64 r = y - conv (a, b); 62 r = y - conv (a, b);
65 else 63 else
66 r = [(zeros (1, lc - ly)), y] - conv (a, b); 64 r = [(zeros (1, lc - ly)), y] - conv (a, b);
67 endif 65 endif
68 66
69 r = polyreduce (r);
70
71 endfunction 67 endfunction