Больше полугода ушло на написание материала, где постарался описать принципы расчета динамических режимов нелинейных цепей с применением BDF и метода Ньютона–Рафсона. Ну вот один из вариантов готов.