chapter x. example 2. plotting the pol,s which make up a b-spline
c from * a practical guide to splines * by c. de boor
calls bsplvb
c
integer ia,left
real biatx(4),t(11),values(4),x
c knot sequence set here....
data t / 4*0.,1.,3.,4.,4*6. /
do 20 ia=1,40
x = float(ia)/5. - 1.
do 10 left=4,7
call bsplvb ( t, 4, 1, x, left, biatx )
c
c according to bsplvb listing, biatx(.) now contains value
c at x of polynomial which agrees on the interval (t(left),
c t(left+1) ) with the b-spline b(left-4 + . ,4,t) . hence,
c biatx(8-left) now contains value of that polynomial for
c b(left-4 +(8-left) ,4,t) = b(4,4,t) . since this b-spline
c has support (t(4),t(8)), it makes sense to run left = 4,
c ...,7, storing the resulting values in values(1),...,
c values(4) for later printing.
c
10 values(left-3) = biatx(8-left)
20 print 620, x, values
620 format(f10.1,4f20.8)
stop
end