Please solve the circuit by using equation 2 and 4 and get graph as same as below. For further explanation please read the journal.
Iijn+1= Iijn + Δt/Lij(Vin+1/2 -Vjn+1/2 - RijIijn + Eijn+1/2). (2)
Vin+1/2 = ((CiVin-1/2)/Δt + Hin - k=1ΣMiIikn)/((Ci/Δt) + Gi) (4)
Attachment:- TRANSACTIONS ON CIRCUITS AND SYSTEMS.pdf