The essential factor in laser frequency conversion involves phase matching within nonlinear optical crystals. To our knowledge, few studies have investigated the noncollinear phase matching calculation for biaxial crystal out of the principal plane. In this paper, we propose an arbitrary direction phase matching model and a computational method based on gradient descent (GD) algorithm, which can be applied to noncollinear in the principal plane, collinear and noncollinear out of the principal plane. In the case of 1053 nm third harmonic generation (THG) in LiB3O5 (LBO) crystal, the phase matching conditions are converted into a system of nonlinear equations with six variables and six equations, which can be solved by iterative optimization search with the GD algorithm and includes type-I (ss-f) and type-II (fs-f). We reveal the relationship of phase matching angles and effective nonlinear coefficients (d(eff)) for various structures. Our method uncovers the existence of many solutions in the non-principal plane with gamma > 8 degrees and the (eff) close to the maximum value 0.66834 pm/V at theta = 90 degrees, phi = 141.84 degrees and gamma = 0. The resolution of the arbitrary direction phase matching problem holds significant importance, as it expands the possibilities for laser frequency conversion, especially for noncollinear structures.