feas_s10.m

function y=feas_s10(M,d1,d2,d3,d4,d5,d6,d7,d8,d9,d10,N_r)
y=1;
if (2*(d1*(M-d1)+d2*(M-d2)+d4*(M-d4)+d5*(M-d5)+d6*(M-d6)+d7*(M-d7)+d8*(M-d8)+d9*(M-d9)+d10*(M-d10))+N_r*M^2<=2*...
        (d1*d2+d1*d3+d1*d4+d1*d5+d1*d6+d1*d7+d1*d8+d1*d9+d1*d10+...
               d2*d3+d2*d4+d2*d5+d2*d6+d2*d7+d2*d8+d2*d9+d2*d10+...
                     d3*d4+d3*d5+d3*d6+d3*d7+d3*d8+d3*d9+d3*d10+...
                           d4*d5+d4*d6+d4*d7+d4*d8+d4*d9+d4*d10+...
                                 d5*d6+d5*d7+d5*d8+d5*d9+d5*d10+...
                                       d6*d7+d6*d8+d6*d9+d6*d10+...
                                             d7*d8+d7*d9+d7*d10+...
                                                   d8*d9+d8*d10+...
                                                         d9*d10))
    y=0;
end