3.1 Implicit Euler Solver qn+1=qn+hvn+1(1) vn+1=vn+hM−1(fint(qn+1+fext)(2) 将(2)代入(1)得 qn+1=qn+hvn+1(1)