自洽求解薛定谔-泊松方程的matlab程序