function [aa,ba,ca,da,adn]=sablockclc(r1,r2,l);
a1=1;
b1=r1;
c1=0;
d1=1;
dn1=1;
a2=[0 1];
b2=[l r2];
c2=[0 0];
d2=[0 1];
dn2=[0 1];
b1pb2=conv(b1,dn2)+conv(b2,dn1);
c1pc2=conv(c1,dn2)+conv(c2,dn1);
a1ma2=conv(a1,dn2)-conv(a2,dn1);
d1md2=conv(d1,dn2)-conv(d2,dn1);
b1b2=conv(b1,b2);
b1d2pd1b2=conv(b1,d2)+conv(d1,b2);
a1b2pb1a2=conv(a1,b2)+conv(b1,a2);
dn12=conv(dn1,dn2);
adn=conv(b1pb2,dn12);
ca=conv(b1pb2,c1pc2)-conv(a1ma2,d1md2);
aa=conv(a1b2pb1a2,dn12);
ba=conv(b1b2,dn12);
da=conv(b1d2pd1b2,dn12);