diff --git a/examples/owl.scad b/examples/owl.scad
new file mode 100644
index 00000000..270af9fc
--- /dev/null
+++ b/examples/owl.scad
@@ -0,0 +1,263 @@
+use <polyhedra/polar_zonohedra.scad>;
+use <polyhedra/octahedron.scad>;
+use <polyhedra/icosahedron.scad>;
+use <shear.scad>;
+use <polyline_join.scad>;
+use <bezier_curve.scad>;
+use <rails2sections.scad>;
+use <sweep.scad>;
+use <experimental/worley_sphere.scad>;
+use <curve.scad>;
+use <ptf/ptf_rotate.scad>;
+use <util/dedup.scad>;
+
+detail = 1;
+head_angles = [5, 0, 20];
+
+owl(detail, head_angles);
+rock();
+
+module owl(detail, head_angles) {
+	n = (detail + 1) * 5;
+	$fn = n;
+
+	module head() {
+		module eye_plane_mask() {
+			translate([10, -66, 82.5])
+			rotate([90, 0, 30])
+			linear_extrude(10)
+				circle(12);
+		}
+
+		module eye_ring() {
+			translate([10, -65, 82.5])
+			rotate([270, 0, 30])
+			scale([1, 1, 2])
+			rotate(180 / $fn)
+			rotate_extrude()
+			translate([12, 0, 0])
+				circle(1.75);
+		}
+
+		module eye_hole_mask() {
+			translate([10, -66, 82.5])
+			rotate([270, 0, 30])
+			rotate(180 / $fn)
+			translate([0, 0, -1])
+			linear_extrude(3, scale = .5)
+			circle(10);
+		}
+
+		module eye_ball() {
+			translate([8, -62, 82.5])
+				icosahedron(4, n / 5 - 1);
+		}
+
+		difference() {
+			union() {
+				difference() {
+					translate([0, -49, 82.5])
+					rotate([15, 0, 0])
+					shear(sy = [0, -.05])
+					scale([25, 25, 25])
+						octahedron(1, n / 5);
+
+					eye_plane_mask();
+					mirror([1, 0, 0])
+						eye_plane_mask();
+				}
+				
+				eye_ring();
+				mirror([1, 0, 0])
+					eye_ring();
+			}
+			
+			eye_hole_mask();
+			mirror([1, 0, 0])
+				eye_hole_mask();
+		}
+
+		eye_ball();
+		mirror([1, 0, 0])
+			eye_ball();
+			
+		eyebrow();
+		mirror([1, 0, 0])
+			eyebrow();
+
+		// beak
+		translate([0, -72.25, 80])
+		scale([.85, 1, 3.5])
+		rotate([25, 0, 0])
+			octahedron(2.5, detail);
+
+	}
+
+	module body() {
+		// back
+		translate([0, 8, 22])
+		rotate([68, 0, 0])
+		shear(sy = [0, .3])
+		scale([11, 10, 15])
+		scale(10 / n)
+			polar_zonohedra(n);
+
+
+		// belly
+		translate([0, -35, 22])
+		rotate([5, 0, 0])
+		scale([10, 9, 12])
+		scale(10 / n)
+			polar_zonohedra(n);
+			
+		translate([0, 7.5, 15])
+		rotate([12, 0, 0])
+		shear(sy = [0, -1])
+		scale([9, 9, 10])
+		scale(10 / n)
+			polar_zonohedra(n);
+		
+		// tail
+		translate([0, 20, 9])
+		rotate([12, 0, 0])
+		shear(sz = [0, .1])
+		shear(sy = [0, -1.1])
+		scale([8, 8, 11])
+		scale(10 / n)
+			polar_zonohedra(n);
+
+		translate([4, 20, 11])
+		rotate([12, 0, 2])
+		shear(sz = [0, .1])
+		shear(sy = [0, -1.1])
+		scale([6, 6, 10])
+		scale(10 / n)
+			polar_zonohedra(n);
+	}
+
+	module wing() {
+		translate([2, 12, 25])
+		rotate([68, 5, 0])
+		shear(sz = [0, .3])
+		shear(sy = [0, .2])
+		scale([12, 10, 14])
+		scale(10 / n)
+			polar_zonohedra(n);
+	}	
+
+	module eyebrow() {
+		t_step = 0.2 / detail;
+		points = bezier_curve(t_step, [
+			[0, -74, 83], 
+			[5, -75, 97], 
+			[15, -70, 90], 
+			[25, -55, 105]
+		]);
+			
+		points2 = bezier_curve(t_step, [
+			[0, -74, 83], 
+			[5, -75, 95], 
+			[20, -65, 82], 
+			[25, -55, 105]
+		]);
+
+		points3 = bezier_curve(t_step, [
+			[0, -74, 83], 
+			[5, -60, 95], 
+			[15, -65, 92], 
+			[25, -55, 105]
+		]);
+
+			
+		points4 = bezier_curve(t_step, [
+			[0, -74, 83], 
+			[-5, -72, 100], 
+			[15, -70, 97], 
+			[25, -55, 105]
+		]);
+
+		rails = [
+			points4, points3, points2, points
+		];
+
+		sections = rails2sections(rails);
+
+		sweep(sections);
+	}
+	
+	module claw() {
+		pts = [
+			[0, 0.55], [1, 0.45], [2.5, 0.375], [6, 0.825], [8, -0.375], 
+			[8, -0.375], [6, 1.875], [4, 1.6], [1.8, 2.5], [1.5, 2.8], [1.2, 3.3], [1.05, 3.8], [1, 4], [0, 8]
+		];
+
+		$fn = 16;
+
+		a = 360 / $fn;
+		x = 6.2 * cos(a);
+		y = 6.2 * sin(a);
+		path = [
+			[0, 0], [2.5, 0], [x, y], [x + 1, y + 1]
+		];
+		path2 = [
+			for(i = len(path) - 1; i > -1; i = i - 1) 
+				ptf_rotate([path[i][0], -path[i][1]], a * 2)
+		];
+
+		t_step = 0.25;
+		claw_path_basic = concat(curve(t_step, path), curve(t_step, path2));
+
+		claw_path1 = [for(p = claw_path_basic) ptf_rotate(p * 1.2, a * 2)];
+		claw_path2 = [for(p = claw_path_basic) ptf_rotate(p * 1.15, a * 4)];
+		claw_path3 = [for(p = claw_path_basic) ptf_rotate(p * 1.1, a * 6)];
+		claw_path4 = [for(p = claw_path_basic) ptf_rotate(p, a * 12)];
+
+		rotate(-15)
+		scale([1.15, 1.3, 1]) 
+		rotate(15)
+		intersection() {
+			rotate_extrude($fn = 7)
+				polygon(pts);
+			linear_extrude(5)
+				polygon(dedup(concat(claw_path1, claw_path2, claw_path3, [[-2, -.75], [-1.45, -1.45]], claw_path4, [[1.45, -1.45], [2, -.75]])));
+		}
+	}
+
+	translate([0, -45, 0])
+	translate([0, 0, 80])
+	rotate(head_angles)
+	translate([0, 45, -80])
+		head(); 	
+		
+	body();
+
+	wing();
+	mirror([1, 0, 0])
+		wing();
+
+	translate([10, -30, 18])
+	scale([1.2, 1.2, 1.5])
+	rotate([10, 0, 180])
+	    claw();
+
+	translate([-10, -30, 18])
+	scale([1.2, 1.2, 1.5])
+	rotate([10, 0, 150])
+	    claw();
+}
+
+module rock() {
+    radius = 16.5;
+    detail = detail;
+    amplitude = .25;
+    dist = "border"; 
+	
+	translate([1, -29, -6.9])
+	scale([1, 1.2, 1])
+	rotate([35, 3, 0])
+    translate([0, 0, 12])
+	scale([3, 1, 1])	
+		worley_sphere(radius, detail, amplitude, dist, seed = 4);
+}
+
+