rodft second derivatives are computed numerically

[edoapra]

Describe the bug
rodft second derivatives are computed numerically.
Instead, they should be computed analytically.

To Reproduce

geometry
 he 0 0  0.9498
 he 0 0 -0.9498
end

basis
 * library 6-31g
end

dft
  rodft
  mult 3
end

task dft freq