Coupled Hartree-Fock computational scheme for nonlinear response tensors in molecules