>> [reg, hilbertFunc] = HilbertFunc(4, 'Points4', 'Points4dual', 1e-10, 0)

reg =

     5


hilbertFunc =

     1
     5
    15
    31
    45
    48

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

reg =

     5


hilbertFunc =

     1
     5
    15
    28
    36
    37