examp7_19.m

《高等应用数学问题MATLAB求解》的PDF版及习题解答

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)  % 绘制各个状态的时间响应曲线