Static and dynamic beta and gamma hyperpolarizability calculations using closed shell TDDFT or TDHF have already been mature before in GAMESS, and you also can calculate directly for many exciting NLO properties according to you choices during one run with the latter. You can calculate static beta and gamma ones using closed shell mp2 by finite field method, too.
'