北航数值分析A课程的第一个大作业,程序有效,结果正确