        Aeq_ik(i,c)=cc;
0132     <span class="keyword">end</span>
0133 <span class="keyword">end</span>
0134 <span class="comment">% Aeq_ik=zeros(nb_solid);</span>
0135 <span class="comment">% beq_ik=zeros(nb_solid,1);</span>
0136 <span class="comment">% for i=1:nb_solid</span>
0137 <span class="comment">%     if size(OsteoArticularModel(i).linear_constraint) ~= [0 0] %#ok&lt;BDSCA&gt;</span>
0138 <span class="comment">%         Aeq_ik(i,i)=-1;</span>
0139 <span class="comment">%         Aeq_ik(i,OsteoArtic