Refactor to easily create new shape function and flux point distribution
Fix ERK schemes