>> [reg, hilbertFunc] = HilbertFunc(3, 'Points', 'Pointsdual', 1e-10, 0)

reg =

     4


hilbertFunc =

     1
     4
    10
    15
    16

>> [reg, hilbertFunc] = HilbertFuncVeronese(3, 'Points', 1e-10, 0)

reg =

     2


hilbertFunc =

     1
     4
     5