Fix error

Author: Leticia-maria

src/optimizebasis.jl

@@ -85,7 +85,8 @@ function optimizebasis(molecule, settings)
 
     for atom in atoms
         for setting in settings
-            for (idx, (n, l, k)) in enumerate(setting[atom.symbol])
+            println("Optimizing basis for atom: ", atom.symbol)
+            for (n, l, k) in enumerate(setting[atom.symbol])
                 β, γ, ζ, alphas, coeffs, rms = optimise_basis_v2(k, n, l)
                 for momentum in _angularmomentum(l)
                     ℓ = momentum[1]