Derivation of R-matrix from local Hamiltonian density

A computer algebra algoritm for solving the quantum Yang-Baxter equation is presented. It is based on the Taylor expansion of R-matrix which is developed up to the order \lambda^6. As an example the classification of 4x4 R-matrices is given.