Татушник, я вспомнил, мы в курсаче по электротехнике численным интергрированием рассчитаывали переходные процессы в цепи постоянного тока с кондером и индуктивностью. Входные значения всех элементов вводились. А вот функции расчета разными методами - были в программе. Писали на с++. Методом Прямоугольников, трапеций, и вроде еще методом Ньютона