求解：一微分方程数值解法问题

[tony211]

x''-A(1-x^2)*x'+Bx=0,该方程描述了一个非线性振荡器，对任何不为0的初始条件最后都会产生等幅震荡，在A大于2的情况下，方程中的第二项动态范围变得相当大，试采用数值积分法对其进行数值仿真，并据此对微分方程的数值解法进行讨论（如RK算法与步长的关系）。（A=5，B=1，X（0）=1.0）。