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

[Bernard Grung]

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


Sorry, I forget to remove the "bg_" from the code in the last comment.
Here is a better way -- just run the following in your simulator:

        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                    ; 
        fi= 1.05e+08     ;
        T =1/fi          ;
        p1=exp(p1*T);
        p2=exp(p2*T);
        p3=1;
        p4=1;
        p =[p1;p2;p3;p4];
        r4=r4*T;
        r3=r3-r4;
        r =[r1;r2;r3;r4];
        m=[       1     ,     1     ,     1         ,     1         ;
               -(2+p2)  ,  -(2+p1)  , -(1+p1+p2)    ,  -(p1+p2)     ;
              +(1+2*p2) , +(1+2*p1) , (p1*p2+p1+p2) ,   +p1*p2      ;
                 -p2    ,   -p1     , -p1*p2        ,     0         ];
        a1=1.00;
        a2=-(2+p1+p2);
        a3=1+2*(p1+p2)+p1*p2;
        a4=-(p1+p2+2*p1*p2);
        a5=p2*p1;
        a =[a1,a2,a3,a4,a5];
        b=[m*r]';
        k=[];
        e=[1;1;1;2];
        [rx,px,kx,ex] = bg_residuez(b,a);
        assert ((abs (rx - r) < 1e-8
              && abs (px - p) < 1e-8
              && isempty (kx)
              && ex == e));



    _______________________________________________________

Reply to this item at:

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

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