這個方程無法求出解析解（絕大多數混沌系統應該都是球不出解析解的吧），只能用數值方法求解。以下給出參考代碼：
f=inline('[x(2); -0.1*x(2)-x(1)^5+6*sin(t)]','t','x')
[t,y]=ode45(f,100,[2; 3]);
comet(y(:,1),y(:,2))代碼很短，看看相關函數的參考應該能搞明白，如果還有疑問再追問吧。
捎帶說壹下，妳出的錯是把“function”寫成了“fuction”。