Numerical application of Szego's method for constructing Liapunov functions