Initial commit

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