% parameters gleak = [0.3] Eleak = [-54.4] % functions Ileak(X) = -gleak.*(X-Eleak) % linkers @current += Ileak(X)