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

[Marco Caliari]

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

To comment #1

I think the code which tries to remove small coefficients should be removed,
since too dangerous. I tried the following in matlab


a=10*[1, -5, 8, -4];
b=[1,0,1];
[r,p,k]=residue(b,a);
[b,a]=residue(r,p,k);


b(2) is about 6.66e-16 and less than max([abs(a),abs(b)])*eps. But it was not
converted to 0 by matlab.

    _______________________________________________________

Reply to this item at:

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

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