本程序的功能是用牛顿拉夫逊法进行潮流计算