3 REAL*8 FUNCTION r8seval3(N, U, X, Y, B, C, D, deriv, deriv2)
7 INTEGER,
PARAMETER :: R8=selected_real_kind(12,100)
11 real*8 u, x(n), y(n), b(n), c(n), d(n),
deriv
48 IF ( u .LT. x(i) )
GO TO 10
49 IF ( u .LE. x(i+1) )
GO TO 30
56 IF ( u .LT. x(k) ) j = k
57 IF ( u .GE. x(k) ) i = k
58 IF ( j .GT. i+1 )
GO TO 20
64 r8seval3 = y(i) + dx*(b(i) + dx*(c(i) + dx*d(i)))
65 deriv = b(i) + dx*(2.0_r8*c(i) + dx*3.0_r8*d(i))
66 deriv2 = 2.0_r8*c(i) + dx*6.0_r8*d(i)
69 zslop=(y(n)-y(n-1))/(x(n)-x(n-1))
71 zslop=(y(i+1)-y(i))/(x(i+1)-x(i))
73 r8seval3=y(i)+dx*zslop