>> [reg, hilbertFunc, stdMon] = HilbertFunc(3, 'Points1', 'Points1dual', 1e-10, 1)

reg =

     1


hilbertFunc =

     1
     2


stdMon =

     0     0     0
     1     0     0

>> [reg, hilbertFunc, stdMon] = HilbertFunc(3, 'Points2', 'Points2dual', 1e-10, 1)

reg =

     2


hilbertFunc =

     1
     3
     4


stdMon =

     0     0     0
     1     0     0
     0     0     1
     1     0     1

>> [reg, hilbertFunc, stdMon] = HilbertFunc(3, 'Points3', 'Points3dual', 1e-10, 1)

reg =

     2


hilbertFunc =

     1
     4
     6


stdMon =

     0     0     0
     1     0     0
     0     1     0
     0     0     1
     1     1     0
     1     0     1

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

reg =

     2


hilbertFunc =

     1
     4
     7


stdMon =

     0     0     0
     1     0     0
     0     1     0
     0     0     1
     1     1     0
     1     0     1
     0     1     1

>> [reg, hilbertFunc, stdMon] = HilbertFunc(3, 'Points5', 'Points5dual', 1e-10, 1)

reg =

     3


hilbertFunc =

     1
     4
     7
     8


stdMon =

     0     0     0
     1     0     0
     0     1     0
     0     0     1
     1     1     0
     1     0     1
     0     1     1
     1     1     1