3 REAL FUNCTION SEVAL2(N, U, X, Y, B, C, D, deriv)
5 REAL U, X(N), Y(N), B(N), C(N), D(N), deriv
38 IF ( u .LT. x(i) )
GO TO 10
39 IF ( u .LE. x(i+1) )
GO TO 30
46 IF ( u .LT. x(k) ) j = k
47 IF ( u .GE. x(k) ) i = k
48 IF ( j .GT. i+1 )
GO TO 20
54 seval2 = y(i) + dx*(b(i) + dx*(c(i) + dx*d(i)))
55 deriv = b(i) + dx*(2.0*c(i) + dx*3.0*d(i))
58 zslop=(y(n)-y(n-1))/(x(n)-x(n-1))
60 zslop=(y(i+1)-y(i))/(x(i+1)-x(i))