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

[anonymous]

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

                 Summary: residue function gives incorrect result
                 Project: GNU Octave
            Submitted by: None
            Submitted on: Fri 07 Oct 2016 12:12:48 PM UTC
                Category: Octave Function
                Severity: 3 - Normal
                Priority: 5 - Normal
              Item Group: None
                  Status: None
             Assigned to: None
         Originator Name: Frank Wiedmann
        Originator Email: address@hidden
             Open/Closed: Open
         Discussion Lock: Any
                 Release: 4.0.3
        Operating System: Microsoft Windows

    _______________________________________________________

Details:

The Octave result


>> clear all
>> rf = [1e3 2e3 1e3 2e3];
>> cf = [316.2e-9 50e-9 31.6e-9 5e-9];
>> [num,den] = residue(1./cf,-1./(rf.*cf),0)
num =

  2.5481e+008  1.5687e+013  2.1808e+017  6.0049e+020

den =

  1.4481e+005  4.9290e+009  4.5817e+013  1.0008e+017

>>

seems to be incorrect (the first element of den is missing). With MATLAB, I
get what seems to be the correct result:


>> format shorte
>> clear all
>> rf = [1e3 2e3 1e3 2e3];
>> cf = [316.2e-9 50e-9 31.6e-9 5e-9];
>> [num,den] = residue(1./cf,-1./(rf.*cf),0)

num =

   2.5481e+08   1.5687e+13   2.1808e+17   6.0049e+20


den =

   1.0000e+00   1.4481e+05   4.9290e+09   4.5817e+13   1.0008e+17

>> 





    _______________________________________________________

Reply to this item at:

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

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