Formulation of Finite Element Equation for 4-Node Quadrilateral Element (2-D Field Problem), intrinsic coordinates and isoparametric transformation