Numerical computing of molecular electrostatics through boundary integral equations