Numerical performance of triangle element approximation for solving 2D Poisson equations using 4-point EDGAOR method