13.5 编制计算三质量扭振系统的固有振型的程序