Task 1:
shells = {p};
info = Initialize[shells, "n", {}, s];
wt = Coulomb[info, 1, 1, 1, 1];
Task 2:
w = wt[[5]] - wt[[6]];
es = EigensystemH[w];
es[[1]]
Task 3:
shells = {d};
states = {};
info = Initialize[shells, jvalues, states, sym];
wt = Coulomb[info, 1, 1, 1, 1];
w = wt[[1]] F0 + (wt[[5]] - wt[[6]]) F2 + (wt[[9]] - wt[[10]]) F4;
en = Expand[Eigenvalues[w]]
The above is 3 tasks in computational physics, kindly run the code as it has some error. will be grateful
Attachments: