examp7_19.m

来自「高等应用数学的问题的matlab求解」· M 代码 · 共 3 行

x0=[1;0;0;1]; xd0=[0;0;1;1];  % 定义 x_0 和 dx_0
res=ode15i('c7mimp',[0,2],x0,xd0);  % 求解隐式微分方程
plot(res.x,res.y)  % 绘制各个状态的时间响应曲线