M = [   r/k    r/k       b;
        r/k*a  r/k       b-eps;
        c*b    c*(b-eps) 0       ];
ss = inv(M) * [r;r;d]