Estadísticas de Hp-adaptabilidad automática aplicando reconstrucción del gradiente por mínimos cuadrados en la resolución numérica de leyes de conservación por el método de Rothe y formulaciones en elementos finitos estabilizados con difusividad implícita