matlab解二阶微分方程 方程:y''=680000(y-60)/(46.05*(1.5-x)), 边界条件:y(0)=70; y'(1.5)=0。可以用数值方法来求解,如变步长的Runge——Kutta法。求解思路:
1、自定义微分方程函数
function f = ode_fun(x,y)
f=[y(2);680000*(y(1)-60)/(46.05*(1.5-x))];
end
2、调用ode45函数
ode45函数使用格式:
[X,Y] = ode45(@ode_fun,tspan,y0)
tspan——自变量初值
y0——初值
运行结果