diff --git a/example/daubechies_wavelets/bootstrap_chebyshev.cpp b/example/daubechies_wavelets/bootstrap_chebyshev.cpp
new file mode 100644
index 000000000..b158fabcd
--- /dev/null
+++ b/example/daubechies_wavelets/bootstrap_chebyshev.cpp
@@ -0,0 +1,18 @@
+#include <iostream>
+#include <boost/math/special_functions/daubechies_scaling.hpp>
+#include <boost/math/special_functions/chebyshev_transform.hpp>
+
+template<typename Real, int p>
+void bootstrap()
+{
+    std::cout << "Computing phi. . .\n";
+    auto phi = boost::math::daubechies_scaling<Real, p>();
+    std::cout << "Computing Chebyshev transform of phi.\n";
+    auto cheb = boost::math::chebyshev_transform(phi, phi.support().first, phi.support().second);
+    std::cout << "Number of coefficients = " << cheb.coefficients().size() << "\n";
+}
+
+int main()
+{
+    bootstrap<long double, 9>();
+}
\ No newline at end of file