PSD-BPA暂态稳定程序,算法:微分方程采用隐式梯形积分法;代数方法采用三角分解法;微分方程和代数方程交替求解