自洽平均场理论的数值方法