隐式微分方程求解。隐式微分方程就是不能转换成显式常微分方程组的微分方程,在Matlab中提供专门

的函数odel5i直接求解隐式微分方程。若隐式微分方程的形式为给定初始条件x(t0)=x0,(to)=x,则可以编写函数描述该隐式微分方程,然后调用命令就可以求解该隐式微分方程。其中,fun为Matlab函数描述隐式微分方程,[t0,tn]为微分方程的求解区间;x0为x(t0)的初始值,xp0为&(t)的初始值。但是隐式微分方程不同于-般的显式微分方程,求解之前,除了给定x(1)的初始值，还需要i(1)的初始值,xi(1)的初始值不能任意赋值,必须满足微分方程的相容性条件,否则将可能出现矛盾的初始值。通常使用函数decic求出这些未完全定义的初值条件,函数decie的使用格式为

其中x0是给定的x(t)的初始值，xp0是任意给定的x(1)的初始值，fixed_:x0和fixed_xp0是与xp0同维数的列向量，其分量为1表示需要保留的初值，为0表示需要求解的初始值。若fixed_x0和fixed_xp0等于空矩阵[]，表示允许所有的初值分量可以发生变化。分别用显式和隐式解法求下列微分方程的数值解