En el Método-P el tamaño del modelo de partida y final es exactamente el mismo, lo único que varía es el orden polinomial de las caras de los elementos. Se ejecuta el análisis estático y se obtienen los siguientes resultados:

	El Método-P resuelve el problema en menos de 1 min. y en cada iteración incrementa el orden polinomial selectivamente en aquellas caras de elementos cuyo error supere el 1% fijado por el usuario. Al final de la 7ª iteración, se obtiene una tensión sx = 83.8 MPa, lo que supone un error del +6.7% frente al valor teórico.
	Con el Método-P se pueden obtener Gráficos de Convergencia, muy útiles para validar los resultados del análisis y conocer la calidad de la solución obtenida y permitir al usuario conocer la evolución en cada iteración de los resultados. Los resultados disponibles para cada Iteración son los siguientes: Tensiones Von Mises Desplazamiento Resultante Energía de Deformación Grados de Libertad (DOF) y Error
	Las ventajas del METODO-P frente al método tradicional de densidad uniforme de malla son evidentes: mayor precisión de resultados, mayor rapidez de análisis y nulo esfuerzo de mallado por parte del usuario. Frente al METODO-H tiene la ventaja que la malla inicial no se modifica en absoluto, por lo que el método es ajeno a problemas de mallado, es mucho más rápido en grandes problemas, y es el método ideal para analizar problemas con mallas creadas en sistemas CAD muchas veces de "dudosa calidad" por contener elementos de elevado "aspect ratio" e incluso elementos muy deformados.
	El METODO-P es más sensible y permite obtener mayor precisión que el METODO-H.
	Otro aspecto que no se ha comentado en favor del METODO-P frente al H es que admite elementos cuadrilátero PLANE2D de 4-8 nodos, cuya precisión es siempre superior a los triángulos de 3-6 nodos.
	En resumen, en algunos problemas el grado de convergencia del METODO-P es generalmente superior al METODO-H; además es más rápido, especialmente en grandes problemas, más robusto por conseguir la mejora de la solución internamente al elemento, la precisión de resultados no es función del tamaño del elemento sino del orden del mismo, y está especialmente indicado para resolver problemas mallados en sistemas CAD.

Modelo Inicial y final con 38 elementos TRIANG de 6-nodos. La figura muestra el orden final polinomial de cada elemento.	Reparto de Tensiones Resultantes sx en el Punto D = 83.8 MPa

Gráfica de Error vs. nº iteraciones	Gráfica de Energía vs. nº iteraciones

Gráfica de Error vs. Desplazamiento resultante	Gráfica de Grados de Libertad (DOF) vs. nº iteraciones

Gráfica de Tensión vs. nº iteraciones