Theoretical Calculation of Nonlinear Optical Susceptibilities