非线性随机常微分方程的matlab模拟用matlab应该怎么写?微分方程为：dx/dt=x-(x^3

非线性随机常微分方程的matlab模拟用matlab应该怎么写?微分方程为：dx/dt=x-(x^3/3)-ydy/dt=x+a+εw其中w是白噪声,ε是噪声的强度我想用matlab画出这个微分方程的关于初值问题的解的轨线的数值模拟,并根据模拟研究解的收敛性质,应该怎样操作?