Update README.md

README.md

@@ -9,7 +9,8 @@ A library of functions for polynomial bases used in spectral element methods usi
 [FastGaussQuadrature.jl](https://github.com/ajt60gaibb/FastGaussQuadrature.jl) for `Float64` and root finding 
 via the Newton algorithm for other scalar types (such as `BigFloat`). The algorithms for interpolation and
 differentiation use barycentric weights as described in the book "Implementing Spectral Methods for PDEs"
-by David Kopriva.
+by David Kopriva. If [SymPy.jl](https://github.com/JuliaPy/SymPy.jl) is loaded, symbolic computations using
+`SymPy.Sym` are supported.
 
 A brief tutorial is given as 
 [notebook](http://nbviewer.ipython.org/github/ranocha/PolynomialBases.jl/blob/master/notebooks/Tutorial.ipynb).