【数学建模】微分方程数值解，求海上缉私问题

实验目的

1. 掌握用MATLAB求常微分方程初值问题数值解的方法。

2. 用欧拉方法和龙格-库塔方法解决简化的实际问题。

实验内容（问题介绍）

算法介绍

调试过程

（1）

function dv = fun(v,~)
weight=527.436*0.4536*9.8;%单位N，实际上11bf=4.44822162N，这里11bf=4.44528N
buoyancy=470.327*0.4536*9.8;%单位N	
K=0.08*0.4536*9.8/0.3048;%阻力比例系数
g=9.8;
dv = g*(weight-buoyancy)/weight-K*g*v/weight;
end
function dy= fun2(v,~)
weight=527.436*0.4536*9.8;%单位N，实际上11bf=4.44822162N，这里11bf=4.44528N
buoyancy=470.327*0.4536*9.8;%单位N    
K=0.08*0.4536*9.8/0.3048;%阻力比例系数
g=9.8;
dy = (v*(weight/g))/(weight-buoyancy-K*v);
end
 
[v,t]=ode45('