Function apply_bc_1d 

pub fn apply_bc_1d(
    u: &mut Vec<f64>,
    left: BoundaryCondition,
    right: BoundaryCondition,
)

Applies 1D boundary conditions (left and right) to a solution vector.