Direi che nel caso di rete con un unico bipolo dinamico ti convenga per prima cosa ottenere il circuito equivalente visto dai morsetti del bipolo dinamico
1, per poi risolvere via discontinuità della variabile di stato.
In questo caso particolare,
$i_L(0^+)=i_L(0^-)+\frac{1}{L}\int_{0^-}^{0^+}v_L(t)\ \text{d}t= \frac{1}{L}\int_{0^-}^{0^+}v_{Th}(t)\ \text{d}t$
dove, ovviamente, nell'integrazione, l'unico termine "utile" sarà quello impulsivo.
Determinate le nuove variabili di stato per t=0+, le successive funzioni del tempo
2 saranno semplicemente quelle relative all'evoluzione libera della rete.
Alternativamente puoi risolvere via Laplace.