RT, 虽然曲线动态和论文一致，但是y的变动的具体大小和论文有小的出入，不太明白为什么。
附运算CODE和运算结果。