1 files changed, 31 insertions, 0 deletions
diff --git a/Donjon/src/DETSPLI2.f b/Donjon/src/DETSPLI2.f
new file mode 100644
index 0000000..1d37d25
--- /dev/null
+++ b/Donjon/src/DETSPLI2.f
@@ -0,0 +1,31 @@
+*DECK DETSPLI2
+      SUBROUTINE DETSPLI2(XA,YA,Y2A,N,X,Y)
+*
+*Parameters: 
+* XA       
+* YA       
+* Y2A      
+* N        
+* X         
+* Y    
+*     
+      DIMENSION XA(N),YA(N),Y2A(N)
+      KLO=1
+      KHI=N
+1     IF (KHI-KLO.GT.1) THEN
+        K=(KHI+KLO)/2
+        IF(XA(K).GT.X)THEN
+          KHI=K
+        ELSE
+          KLO=K
+        ENDIF
+      GOTO 1
+      ENDIF
+      H=XA(KHI)-XA(KLO)
+      IF (H.EQ.0.) CALL XABORT('DETSPLI2: BAD XA INPUT.')
+      A=(XA(KHI)-X)/H
+      B=(X-XA(KLO))/H
+      Y=A*YA(KLO)+B*YA(KHI)+
+     *      ((A**3-A)*Y2A(KLO)+(B**3-B)*Y2A(KHI))*(H**2)/6.
+      RETURN
+      END