i = 1;
ourbs(:,20) = [0 0 0 0 0 0 0 0 0 0 0 0 0 0 0 0 0 0]';
for a=1:20
  for q=1:20
    ourbs(:,i) = newb - [0 0 0 0 0 0 0 0 q q 0 0 0 a a 0 0 0]';
    i = i + 1;
  end
end