Integral equation modeling of three-dimensional magnetotelluric response