def poly2col(pol):
    """
    Returns the coefficients (columns) of a polynomial in the form of a list.

    Parameters
    ----------
    pol : NumPy nd-array of shape (rows, num_Lagr_nodes)
        Coefficients, i.e., the ones obtained from interpolate_P1 or interpolate_P2.

    Returns
    -------
    List
        Coefficients stored in column-wise format. 

    """
    rows = pol.shape[0]
    cols = pol.shape[1]
    coeff_list = []
    
    for col in range(cols):
        coeff_list.append(pol[:, col].reshape(rows, 1))

    return coeff_list