The 4-EGAOR method for solving triangle element approximation of 2D poisson equations