  We give a simple and practical algorithm to compute the link polynomials,
which are defined according to the skein relations. Our method is based on a
new total order on the set of all braid representatives. As by-product a new
complete link invariant are obtained.
