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

reg =

     6


hilbertFunc =

           1
          13
          87
         403
        1454
        4342
        8652