※ 引述《jurian0101 (Hysterisis)》之銘言：
: 直接貼到Mathematica讀
: line1=L*D[i[L][t],t]==V[g][t]-r[L]*i[i][t]-v[0][t];
: line2=i[i][t]==(L*D[i[L][t],t])/R[S]+i[L][t];
: line3=v[0][t]==c*r[c]*D[v[c][t],t]+v[c];
: line4=L*D[i[L][t],t]==-(((r[c]*R[S])*(c*D[v[c][t],t]))/(r[L]+R[S]))+
: (R[S]*V[g][t])/(r[L]+R[S])-(v[c]*R[S])/(r[L]+R[S])+
: (i[L][t]*((-r[L])*R[S]))/(r[L]+R[S]);
: line5=c D[v[c][t],t]==R[L]/(r[c]+R[L])(1/R[S]*L*D[i[L][t],t]+i[L][t])
: -v[c]/(r[c]+R[L])
(* 由方程組直接求解 *)
myans1=L*Derivative[1][i[L]][t]/.
Solve[{line1,line2,line3,line5},{i[L]'[t],i[i][t],v[0][t],
v[c]'[t]}][[1]]
(* Eq 2, 3 代入 Eq 1, 再將 Eq 5代入 Eq 1，最後用Solve做移項整理 *)
myans2=L*Derivative[1][i[L]][t]/.
Solve[(line1/.({line2,line3}/.Equal->Rule))/.(line5/.
Equal->Rule),Derivative[1][i[L]][t]][[1]]
(* 係數 *)
mycoef=Table[FullSimplify@Coefficient[z,#]&/@{i[L][t],v[c],
V[g][t]},{z,{myans1,myans2}}]
SameQ@mycoef