examp7_19.m

来自「《高等应用数学问题MATLAB求解》的PDF版及习题解答」· 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)  % 绘制各个状态的时间响应曲线