-->Ab = [1 -2 1]; dAb = 2; zBb = [0 1 1]; dzBb = 2;
 
-->phi = [1 -1.2 0.52]; dphi = 2;
 
-->[S1,dS1,R1,dR1] = xdync(zBb,dzBb,Ab,dAb,phi,dphi)
 dR1  =
 
    1.  
 R1  =
 
    1.    0.32  
 dS1  =
 
    1.  
 S1  =
 
    0.48  - 0.32  
 
-->diary off