怎样用matlab求解非线性微分方程

我用解析法试了，不行。
>> dsolve('D2x=3*x^2+Dx+3*x^5+3*x^2*x')
到底是哪一个啊？
>> dsolve('D2x=3*x^2+Dx+3*x^5+3*x^(2*x)')

告诉初始条件，用数值解法吧。

>> dsolve('D2y=3*y^2+Dy+3*y^5+3*y^2*y')

ans =

&where(_a,{[diff(_b(_a),_a)*_b(_a)-3*_a^2-_b(_a)-3*_a^5-3*_a^3 = 0, _b(_a) = diff(y(t),t), _a = y(t), t = Int(1/_b(_a),_a)+C1, y(t) = _a]})

>>