function [Qz, Qp, Qzp]=finn_Qz_Qp_Qzp(dist_zeros, dist_poles)
    for ii=1:length(dist_zeros)
        for jj=1:length(dist_zeros)
            Qz(ii,jj)=1/(dist_zeros(ii)+conj(dist_zeros(jj)));
        end
    end

    for ii=1:length(dist_poles)
        for jj=1:length(dist_poles)
            Qp(ii,jj)=1/(conj(dist_poles(ii))+dist_poles(jj));
        end
    end

    for ii=1:length(dist_zeros)
        for jj=length(dist_poles)
            Qzp(ii,jj)=1/(dist_zeros(ii)-dist_poles(jj));
        end
    end
end