Gauss-Jacobi method