--- /dev/null
+
+// How to emulate Faust's seq, par, sum.
+// x(k) is assumed to yield the kth signal.
+
+xseq(1,x) = x(0);
+xseq(n,x) = xseq(n-1,x) : x(n-1);
+
+xpar(1,x) = x(0);
+xpar(n,x) = xpar(n-1,x) , x(n-1);
+
+xsum(1,x) = x(0);
+xsum(n,x) = xsum(n-1,x) + x(n-1);
+
+// These are all very similar. Abstracting
+// on the binary "accumulator" function, we
+// get the familiar fold(-left) function:
+
+fold(1,f,x) = x(0);
+fold(n,f,x) = f(fold(n-1,f,x),x(n-1));
+
+// Now seq, par, sum can be defined as:
+
+fseq(n) = fold(n,\(x,y).(x:y));
+fpar(n) = fold(n,\(x,y).(x,y));
+fsum(n) = fold(n,+);
+
+// Often it is more convenient to specify
+// parameters as a Faust tuple. We can match
+// against the (x,y) pattern to decompose
+// these.
+
+vfold(f,(x,y)) = f(vfold(f,x),y);
+vfold(f,x) = x;
+
+// Tuple version of seq, par, sum:
+
+vseq = vfold(\(x,y).(x:y));
+vpar = vfold(\(x,y).(x,y));
+vsum = vfold(+);
+
+// Example: sum of sinusoids.
+
+import("music.lib");
+
+f0 = 440;
+a(0) = 1;
+a(1) = 0.5;
+a(2) = 0.3;
+h(i) = a(i)*osc((i+1)*f0);
+
+v = hslider("vol", 0.3, 0, 1, 0.01);
+
+process = v*fsum(3,h);
+//process = v*xsum(3,h);
+//process = v*vsum((h(0),h(1),h(2)));