Write add_bias(matrix, bias) that broadcasts a per-column bias vector across every row of matrix. The bias length equals the number of columns; add bias[i] to column i of each row. Return the new matrix.