回路シミュレーション ロマンある
「t=0で単位ステップ波を入力したときの過渡応答」以外を解く必要がなければ初期値さえ確定しておけば手計算で式が出せるので、これをアルゴリズムに展開すればよいはずです。
あとは、実際に使われている回路シミュレータの原理あたりを持ち込むのも手かも。概要をつかめる記事を紹介します:
SPICEの内側を探る――節点法とは:SPICEの仕組みとその活用設計(1)(1/3 ページ) - MONOist (
https://monoist.itmedia.co.jp/mn/articles/1304/19/news010.html )