これを整理すると
(dX/dT)=(-(R1+R4+R6)*X+R4*Y+IT*W*COS(W*T+THETA)*(-L6)-R6*IT*SIN(W*T+THETA))/(L1+L6)
(dY/dT)=(-(R2+R4+R5+R7)*Y+R4*X+R5*Z+IT*W*COS(W*T+THETA)*(-L7)-R7*IT*SIN(W*T+THETA))/(L2+L7)
(dZ/dT)=(-(R3+R5+R8)*Z+R5*Y+IT*W*COS(W*T+THETA)*(-L8)-R8*IT*SIN(W*T+THETA))/(L3+L8)

という上の式がでます。
さっき言い忘れましたけどＴというのが時間です。
電源の式が
IT*sin(W*T+THETA)で与えられているので、
(d It/dT)=IT*W*COS(W*T*THETA)　になります。
Ｗは角速度です。ＴＨＥＴＡはＴ＝０の時の位相です。