3 REAL FUNCTION SEVAL(N, U, X, Y, B, C, D)
5 REAL U, X(N), Y(N), B(N), C(N), D(N)
32 IF ( u .LT. x(i) )
GO TO 10
33 IF ( u .LE. x(i+1) )
GO TO 30
40 IF ( u .LT. x(k) ) j = k
41 IF ( u .GE. x(k) ) i = k
42 IF ( j .GT. i+1 )
GO TO 20
47 seval = y(i) + dx*(b(i) + dx*(c(i) + dx*d(i)))