add time

Author: dlfivefifty

README.md

@@ -63,6 +63,15 @@ julia> J[1000] - besselj(999,z) # matches besselj to high (relative) accuracy
 julia> J[11_000] - besselj(11_000-1, z)
 3.3730094946097293e-143
 ```
+We're even faster than SpecialFunctions.jl for constructing a range of Bessel functions:
+```julia
+julia> @time [besselj(k-1, z) for k=0:11_000-1];
+  0.188690 seconds (77.20 k allocations: 3.295 MiB)
+
+julia> @time J = A \ Vcat([besselj(1,z)], Zeros(∞));
+  0.036701 seconds (406.40 k allocations: 46.552 MiB, 24.43% gc time)
+```
+
 
 ## Infinite-dimensional QL factorization