function [Gds11 Gdm11 Gdms11]=gmsd11(Gd11)

ztot=zero(Gd11); ptot=pole(Gd11);
z=ztot(find(ztot>0)); p=ptot(find(ptot>0));
nz=length(z); np=length(p);

s=TF('s');
b=1;a=1;     
for i = 1:np
    b=b*(s-p(i))/(s+p(i));
end;
for i = 1:nz
    a11=a*(s+z(i))/(s-z(i));
end    

Gds11=minreal(tf(b*Gd11));Gdm11=minreal(tf(Gd11*a)); Gdms11=minreal(tf(b*Gd11*a));