diff --git a/Project.toml b/Project.toml
index 8340878..c226aa5 100644
--- a/Project.toml
+++ b/Project.toml
@@ -1,6 +1,6 @@
 name = "FastGaussQuadrature"
 uuid = "442a2c76-b920-505d-bb47-c5924d526838"
-version = "1.0.2"
+version = "1.0.3"
 
 [deps]
 LinearAlgebra = "37e2e46d-f89d-539d-b4ee-838fcccc9c8e"
diff --git a/src/besselroots.jl b/src/besselroots.jl
index 37aff8f..635a07a 100644
--- a/src/besselroots.jl
+++ b/src/besselroots.jl
@@ -137,9 +137,7 @@ function piessens(ν::Float64)
     C = PIESSENS_C
     T = _piessens_chebyshev30(ν)
     y = C'*T
-    _y = collect(y)
-    _y[1] *= sqrt(ν+1)  # Scale the first root.
-    return _y
+    return Base.setindex(y, y[1] * sqrt(ν+1), 1)
 end