Log interpretation relies on measured parameters (resistivity, bulk density, travel time, radioactivity, etc) to obtain a quantitative evaluation of petrophysical properties (porosity, saturation, permeability). The interpretation of traditional induction logging data is not always reliable due to the insensitivity to anisotropy. Triaxial measurements are proposed to overcome this shortcoming and have sufficient sensitivity to anisotropy. In this study we will describe an inversion algorithm for retrieving information of anisotropic conductivities, bed boundary dimensions, dip angle and rotation angle. This approach is based on the Gauss-Newton method that employs a quadratic model of the cost function. Cost function is constructed by means of least squares between the field logs and computed inversion logs. The step length is adjusted by line search in order to effectively decrease mismatch between measured and predicted responses. Global constrain method is used to impose bounds on the inverted parameters. The inversion algorithm will be validated against synthetic data with accurate and consistent log results.