各位好
目前在利用ode解簡單的力學問題
方程式為 F=-kx
當我的dt設置的範圍再0至10(時間很短的情況下)曲線在+-1震盪
但當我dt設置範圍改為0至100時(時間拉長了)曲線後來就不在+-1震盪了
想詢問是否有方法解決這個問題
感謝
以下是我寫
function dydt=odefcn(t,y,m,b,k,w)
dydt=zeros(2,1);
dydt(1)=y(2);
dydt(2)=-(k/m)*y(1);

依舊在+-1之間阿

又改成 dt=0:0.01:100嗎