这个方程怎样用matlab编程实现

lhjm

dx/dt=(-du(x)/x+Asin(2pi*f*t))+sqrtD*n(t ）
其中U（x）=-a*x^2/2+b*x^4/4; n(t)为高斯白噪声；
还得怎么画出曲线来。哪位好心人帮忙看看，谢谢！