add some docs on methods

Author: jverzani

README.md

@@ -5,27 +5,34 @@ Basic arithmetic, integration, differentiation, evaluation, and root finding ove
 [![Build Status](https://travis-ci.org/Keno/Polynomials.jl.png?branch=master)](https://travis-ci.org/Keno/Polynomials.jl)
 
 #### Poly{T<:Number}(a::Vector)
+
 Construct a polynomial from its coefficients, lowest order first.
+
 ```julia
 julia> Poly([1,0,3,4])
 Poly(1 + 3x^2 + 4x^3)
 ```
 
 An optional variable parameter can be added.
+
 ```julia
 julia> Poly([1,2,3], :s)
 Poly(1 + 2s + 3s^2)
 ```
 
 #### poly(r::AbstractVector)
-Construct a polynomial from its roots. This is in contrast to the `Poly` constructor, which constructs a polynomial from its coefficients.
+
+Construct a polynomial from its roots. This is in contrast to the
+`Poly` constructor, which constructs a polynomial from its
+coefficients.
+
 ```julia
 // Represents (x-1)*(x-2)*(x-3)
 julia> poly([1,2,3])
 Poly(-6 + 11x - 6x^2 + x^3)
 ```
 
-#### +, -, *, /, ==
+#### +, -, *, /, div, ==
 
 The usual arithmetic operators are overloaded to work on polynomials, and combinations of polynomials and scalars.
 ```julia
@@ -44,14 +51,18 @@ Poly(3 + 2x)
 julia> p - q
 Poly(2x + x^2)
 
-julia> p*q
+julia> p * q
 Poly(1 + 2x - x^2 - 2x^3)
 
-julia> q/2
+julia> q / 2
 Poly(0.5 - 0.5x^2)
+
+julia> q ÷ p      # `div`, also `rem` and `divrem`
+Poly(0.25 - 0.5x)
 ```
 
 Note that operations involving polynomials with different variables will error.
+
 ```julia
 julia> p = Poly([1, 2, 3], :x)
 julia> q = Poly([1, 2, 3], :s)
@@ -60,6 +71,7 @@ ERROR: Polynomials must have same variable.
 ```
 
 To get the degree of the polynomial use `degree` method
+
 ```
 julia> degree(p)
 1
@@ -72,6 +84,7 @@ julia> degree(p-p)
 ```
 
 #### polyval(p::Poly, x::Number)
+
 Evaluate the polynomial `p` at `x`.
 
 ```julia
@@ -80,7 +93,11 @@ julia> polyval(Poly([1, 0, -1]), 0.1)
 ```
 
 #### polyint(p::Poly, k::Number=0)
-Integrate the polynomial `p` term by term, optionally adding constant term `k`. The order of the resulting polynomial is one higher than the order of `p`.
+
+Integrate the polynomial `p` term by term, optionally adding constant
+term `k`. The order of the resulting polynomial is one higher than the
+order of `p`.
+
 ```julia
 julia> polyint(Poly([1, 0, -1]))
 Poly(x - 0.3333333333333333x^3)
@@ -90,14 +107,20 @@ Poly(2.0 + x - 0.3333333333333333x^3)
 ```
 
 #### polyder(p::Poly)
-Differentiate the polynomial `p` term by term. The order of the resulting polynomial is one lower than the order of `p`.
+
+Differentiate the polynomial `p` term by term. The order of the
+resulting polynomial is one lower than the order of `p`.
+
 ```julia
 julia> polyder(Poly([1, 3, -1]))
 Poly(3 - 2x)
 ```
 
 #### roots(p::Poly)
-Return the roots (zeros) of `p`, with multiplicity. The number of roots returned is equal to the order of `p`. The returned roots may be real or complex.
+
+Return the roots (zeros) of `p`, with multiplicity. The number of
+roots returned is equal to the order of `p`. The returned roots may be
+real or complex.
 
 ```julia
 julia> roots(Poly([1, 0, -1]))
@@ -114,4 +137,41 @@ julia> roots(Poly([0, 0, 1]))
 2-element Array{Float64,1}:
  0.0
  0.0
+ ```
+
+#### Polyfit
+
+* `polyfit`: fits a polynomial of minimal degree fitting the points
+  specified by `x` and `y` using least squares fit.
+
 ```
+julia> xs = 1:4; ys = exp(xs); polyfit(xs, ys)
+Poly(-7.717211620141281 + 17.9146616149694x - 9.77757245502143x^2 + 2.298404288652356x^3)
+```
+
+#### Other methods
+
+Polynomial objects also have other methods:
+
+* 0-based indexing is used to extract the coefficients of $a_0 + a_1
+  x + a_2 x^2 + ...$, coefficients may be changed using indexing
+  notation.
+
+* `coeffs`: returns the entire coefficient vector
+
+* `degree`: returns the polynomial degree, `length` is 1 plus the degree
+
+* `variable`: returns the polynomial symbol as a degree 1 polynomial
+
+* `norm`: find the `p`-norm of a polynomial
+
+* `truncate`: set to 0 small terms in a polynomial; `chop` chops off
+  any small leading values that may arise due to floating point
+  operations.
+
+* `gcd`: greatest common divisor of two polynomials.
+
+
+* `Pade`: Return the
+  [Pade approximant](https://en.wikipedia.org/wiki/Pad%C3%A9_approximant)
+  of order `m/n` for a polynomial as a `Pade` object.