Merge pull request #17 from cfbaptista/patch-1

devectorise 'poly' for better performance

Author: jverzani

src/Polynomials.jl

@@ -305,7 +305,9 @@ function poly{T}(r::AbstractVector{T}, var=:x)
     c = zeros(T, n+1)
     c[1] = 1
     for j = 1:n
-        c[2:j+1] = c[2:j+1]-r[j]*c[1:j]
+        for i = j:-1:1
+            c[i+1] = c[i+1]-r[j]*c[i]
+        end
     end
     return Poly(reverse(c), var)
 end