Proporcionamos un paquete numérico para el cálculo de un diseño de regresión polinómica óptima cercana de grado en un espacio de diseño finito, mediante pocas iteraciones de un algoritmo multiplicativo básico seguido de una compresión tipo Tchakaloff de la medida discreta manteniendo la eficiencia G alcanzada, a través de una versión acelerada del algoritmo de Lawson-Hanson para problemas de Mínimos Cuadrados No Negativos (NNLS). Este paquete puede resolver en una computadora personal problemas a gran escala donde es de hasta -, siendo . Se presentan varios tests numéricos en formas complejas en y en hipercubos en .