aboutsummaryrefslogtreecommitdiff
path: root/math/deboor/progs/prog7.f
diff options
context:
space:
mode:
Diffstat (limited to 'math/deboor/progs/prog7.f')
-rw-r--r--math/deboor/progs/prog7.f17
1 files changed, 17 insertions, 0 deletions
diff --git a/math/deboor/progs/prog7.f b/math/deboor/progs/prog7.f
new file mode 100644
index 00000000..e08ab120
--- /dev/null
+++ b/math/deboor/progs/prog7.f
@@ -0,0 +1,17 @@
+chapter x. example 4. construction of a b-spline via bvalue
+c from * a practical guide to splines * by c. de boor
+calls bvalue(interv)
+ integer ia
+ real bcoef(1),t(5),value,x
+ real bvalue
+c set knot sequence t and b-coeffs for b(1,4,t)
+ data t / 0.,1.,3.,4.,6. / , bcoef / 1. /
+c evaluate b(1,4,t) on a fine mesh. on (0,6), the values should
+c coincide with those obtained in example 3 .
+ do 20 ia=1,40
+ x = float(ia)*.2 - 1.
+ value = bvalue ( t, bcoef, 1, 4, x, 0 )
+ 20 print 620, x, value
+ 620 format(f10.1,f20.8)
+ stop
+ end