function B = bldB(K, nbc, bc)
B = zeros(K,1);
for i = 1:nbc % boundary condition enforcement
    m = bc(i,1);
    B(m) = bc(i,2);
end
return;