这是个含有一个电容一个电感的二阶电路动态方程吧,怎么会有四阶.
这个方程不用龙格库塔法解啊,龙格库塔法等方法属于数值解法,针对微分方程无法求出解析形式时,只能求得近似解.
但这个方程很简单,在给出初始条件iL(0)和Uc(0)的情况下直接能求得解析形式
不妨记方程组为x'=Ax+bUs,这里x=[iL Uc]’,那么这个方程的解析解就是x=exp(A*t)x(0)+与u有关的积分项（百度打公式不方便,又不让给链接,你自己搜一下"线性定常系统的解"）
这个解析解是众所周知的,如果你非要用matlab算也很简单,用fsolve就可以了