1 subroutine bvalus(n,knot,coef,nk,x,s,order)
2C Args
3 integer n, nk, order
4 double precision knot(*),coef(*),x(*),s(*)
5C Local
6 double precision bvalue
7 integer i
8
9 do 10 i=1,n
10 s(i)=bvalue(knot,coef,nk,4,x(i),order)
11 10 continue
12 return
13 end
14