4 REAL*8 FUNCTION r8speval(N, U, X, Y, B, C, D)
6 INTEGER,
PARAMETER :: R8=selected_real_kind(12,100)
8 real*8 u, x(n), y(n), b(n), c(n), d(n)
34 IF ( u .LT. x(i) )
GO TO 10
35 IF ( u .LE. x(i+1) )
GO TO 30
42 IF ( u .LT. x(k) ) j = k
43 IF ( u .GE. x(k) ) i = k
44 IF ( j .GT. i+1 )
GO TO 20
49 r8speval = b(i) + dx*(2._r8*c(i) + 3._r8*dx*d(i))