for ind_me=1:length(me),
    for ind_sc=1:length(sc),
		[w1,b1,w2,b2,nr,err] = solverb(PN,T,[df me(ind_me) eg sc(ind_sc)]);
        a = simurb(PN,w1,b1,w2,b2);
        wpk(ind_me, ind_sc)=(1-sum(abs(T-a)>=0.5)/length(PN))*100;
        uzyt(ind_me, ind_sc)=nr(length(nr));
    end

end