OPTIMAL SOLUTION OF RICHARDS’ EQUATION FOR SLOPE INSTABILITY ANALYSIS USING AN INTEGRATED ENHANCED VERSION OF BLACK HOLE MECHANICS INTO THE FEM