[Octave-bug-tracker] [bug #34266] Four issues with "residue" and residuez"

[Bernard Grung]

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


The main problem is "residuez" -- how do I get connected with Octave-Forge?

Here is one for "residue"
        z1=+7.0372976777e+06;
        p1=-3.1415926536e+09;
        p2=-4.9964813512e+08;
        r1=-(1+z1/p1)/(1-p1/p2)/p2/p1  ;
        r2=-(1+z1/p2)/(1-p2/p1)/p2/p1  ;
        r3=(1+(p2+p1)/p2/p1*z1)/p2/p1  ;
        r4=z1/p2/p1                    ; 
        r=[r1;r2;r3;r4];
        p=[p1;p2;0;0];
        k=[];
        e=[1;1;1;2];
        b=[1,z1];
        a=[1,-(p1+p2),p1*p2,0,0];
        [br, ar] = residue (r, p, k, e);
        assert ((abs (br - b) < 1e-8
              && abs (ar - a) < 1e-8));


    _______________________________________________________

Reply to this item at:

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

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