With polar=enonly you calculate the static polarizability and hyper, with numeric methods you get these derivatives and with CPHF you get the same derivatives but with analityc form, in theory you must get values very similar, the advantage of CPHF is that also you can calculate the dinamic polarizability and hyper this mean that you can get values at different frequency.
you need put the keyword polar i did a example simple with your input with a basis 3-21g, i add this input and log
in the image you can see the results of polarizability the first is the static polarizability and the second the dinamic the iso= isotropic polarizability and the aniso= anisotropic, you can compare with your previus calculation with polar=enonly the results must be similar to the static polarizability (the first part), you can recalculate with the correct basis set