From e4921933bad0110256046e99515b2f4bbad0bd21 Mon Sep 17 00:00:00 2001
From: h_vogt <h_vogt>
Date: Sun, 9 Jun 2013 11:06:05 +0200
Subject: [PATCH] nicomcof.c: xmu used to calculate trap integration
 coefficient

---
 src/maths/ni/nicomcof.c | 12 ++++++++----
 1 file changed, 8 insertions(+), 4 deletions(-)

diff --git a/src/maths/ni/nicomcof.c b/src/maths/ni/nicomcof.c
index 7dcd7ed38..f9d2a9305 100644
--- a/src/maths/ni/nicomcof.c
+++ b/src/maths/ni/nicomcof.c
@@ -7,6 +7,12 @@ Author: 1985 Thomas L. Quarles
 #include "ngspice/cktdefs.h"
 #include "ngspice/sperror.h"
 
+/* xmu=0:    Backward Euler
+ * xmu=0.5:  trapezoidal (standard)
+ * xmu=0.48: good damping of current ringing, e.g. in R.O.s.
+ */
+#define xmu 0.5
+
 
 int
 NIcomCof(CKTcircuit *ckt)
@@ -35,10 +41,8 @@ NIcomCof(CKTcircuit *ckt)
             break;
 
         case 2:
-            ckt->CKTag[0]=1.0/ckt->CKTdelta/(1.0-0.5) ;
-            ckt->CKTag[1]=0.5/(1.0 - 0.5) ;
-            /* above lines should have 'xmu' instead of .5 eventually */
-            /* (in all three places) */
+            ckt->CKTag[0] = 1.0 / ckt->CKTdelta / (1.0 - xmu);
+            ckt->CKTag[1] = xmu / (1.0 - xmu);
             break;
 
         default: