Supported by the National Natural Science Foundation of China(11074184);the Foundation for University Key Young Teacher of Henan Province(2009GGJS-163)
Renormalization group recursions are obtained by virtue of the variational cumulant expansion method. Good qualitative estimates are obtained for the d=2 square Ising system.