To solve PDE's requiring high order base functions but still using triangular 3-noded elements,

use NEM...