The paper reviews the history of B-spline methods for atomic structure calculations bound states. It highlights various aspects variational method, particularly with regard to orthogonality requirements, iterative self-consistent eigenvalue problem, and related sphf, dbsr-hf, spmchf programs. B-splines facilitate mapping solutions from one grid another. following describes a two-stage approach ...