obj = HOKE(caNl, iOrder, ...)
  caNl is the n_l in [m2/W^l]
 
  iOrder = 1   is SPM
  iOrder = 2   n_l is n4
  options ae same as in SPM