[Octave-bug-tracker] [bug #49291] residue function gives incompatible result

[Rik]

Follow-up Comment #7, bug #49291 (project octave):

@Marco: Could you try the following code in Matlab?  It is from one of the
BIST tests for residue.  If Matlab is really not applying any filtering for
small results than br will be a 4-element vector where the first entry is
nearly 0.


b = [1, 0, 1];
a = [1, 0, 18, 0, 81];
[r, p, k, e] = residue (b, a);
r1 = [-5i; 12; +5i; 12]/54;
p1 = [+3i; +3i; -3i; -3i];
assert (r, r1, 1e-12);
assert (p, p1, 1e-12);
assert (isempty (k));
assert (e, [1; 2; 1; 2]);
[br, ar] = residue (r, p, k);
assert (br, b, 1e-12);
assert (ar, a, 1e-12);




    _______________________________________________________

Reply to this item at:

  <http://savannah.gnu.org/bugs/?49291>

_______________________________________________
  Message sent via/by Savannah
  http://savannah.gnu.org/