根据星系巡天数据,拟合得到光度函数的算法程序