22 <span class="keyword">for</span> i=1:length(solid_red)
0123     jj=solid_red(i);
0124     <span class="keyword">if</span> size(OsteoArticularModel(jj).linear_constraint) ~= [0 0] <span class="comment">%#ok&lt;BDSCA&gt;</span>
0125         Aeq_ik(i,i)=-1;
0126         ind_col = OsteoArticularModel(jj).linear_constraint(1,1);
0127         [~,c]=find(GC.q_map(ind_col,:));
0128         
0129         ind_val = OsteoArticularModel(jj).linear_constraint(2,1);
0130         [~,cc]=find(GC.q_map(ind_val,:));
0131 