我要解一個聯立微分方程的數值解並繪圖
題目是
X''=-0.0675*X'*sqrt(X'^2+Y'^2)
Y''=-9.81-0.0675*Y'*sqrt(X'^2+Y'^2)
初始條件X'(0)=134 X(0)=Y(0)=Y'(0)=0
我降階後的程式如下
function dx=odedx(t,x)
a=0.06725;
g=9.81;
%x(1)=x
%x(2)=x'
%x(3)=y
%x(4)=y'
dx=[x(2);x(4);-a*x(2)*sqrt((x(2))^2+(x(4))^2);-g-a*x(4)*sqrt((x(2))^2+(x(4))^2
)];
end

dx順序不對

好習慣 把每個函數都印出來看一下 確認跟手算的都一樣再繼續下一步