我要解一个联立微分方程的数值解并绘图
题目是
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顺序不对

好习惯 把每个函数都印出来看一下 确认跟手算的都一样再继续下一步