--- ray/src/gen/genworm.c	1989/02/02 11:16:31	1.1
+++ ray/src/gen/genworm.c	1990/07/09 09:55:50	1.5
@@ -1,8 +1,10 @@
-/*
-
 #ifndef lint
 static char SCCSid[] = "$SunId$ LBL";
 #endif
+
+/* Copyright (c) 1989 Regents of the University of California */
+
+/*
  *  genworm.c - program to generate worms (strings with varying thickness).
  *
  *	The program takes as input the functions of t for x, y,
@@ -28,13 +30,14 @@ static char SCCSid[] = "$SunId$ LBL";
 #define  max(a,b)	((a) > (b) ? (a) : (b))
 
 
-double  funvalue(), l_hermite(), argument();
+double  funvalue(), l_hermite(), l_bezier(), l_bspline(), argument();
 
 
 main(argc, argv)
 int  argc;
 char  *argv[];
 {
+	extern long	eclock;
 	char  stmp[256];
 	double  t, f, lastr, r;
 	FVECT  lastp, p;
@@ -42,6 +45,8 @@ char  *argv[];
 
 	varset("PI", PI, NULL);
 	funset("hermite", 5, l_hermite);
+	funset("bezier", 5, l_bezier);
+	funset("bspline", 5, l_bspline);
 
 	if (argc < 8)
 		goto userror;
@@ -67,6 +72,7 @@ char  *argv[];
 		goto userror;
 
 	printhead(argc, argv);
+	eclock = 0;
 
 	for (i = 0; i <= nseg; i++) {
 		t = (double)i/nseg;
@@ -160,4 +166,30 @@ l_hermite()			
 		argument(2)*(-2.0*t+3.0)*t*t +
 		argument(3)*((t-2.0)*t+1.0)*t +
 		argument(4)*(t-1.0)*t*t );
+}
+
+
+double
+l_bezier()
+{
+	double  t;
+
+	t = argument(5);
+	return(	argument(1) * (1.+t*(-3.+t*(3.-t))) +
+		argument(2) * 3.*t*(1.+t*(-2.+t)) +
+		argument(3) * 3.*t*t*(1.-t) +
+		argument(4) * t*t*t );
+}
+
+
+double
+l_bspline()
+{
+	double  t;
+
+	t = argument(5);
+	return(	argument(1) * (1./6.+t*(-1./2.+t*(1./2.-1./6.*t))) +
+		argument(2) * (2./3.+t*t*(-1.+1./2.*t)) +
+		argument(3) * (1./6.+t*(1./2.+t*(1./2.-1./2.*t))) +
+		argument(4) * (1./6.*t*t*t) );
 }