B = [ 1 3 2 ; 0 -4 -3 ; 3 1 3 ]

for j=2:3,
  for i=j:3,
    B(i,:) = B(i,:) - B(j-1,:)*B(i,j-1)/B(j-1,j-1)
  end
end