bgodard

https://github.com/JuliaPy/PyPlot.jl/issues/291

```julia julia> import Base:* julia> *(x::Unitful.AbstractQuantity,y::Num) = Quantity(x.val*y, unit(x)) * (generic function with 735 methods) julia> *(y::Num,x::Unitful.AbstractQuantity) = x*y * (generic function with 736 methods) julia> 10u"km"*x 10x km julia>...

Related, I see that the interpolation used is Akima interpolation. I guess this is: http://www.iue.tuwien.ac.at/phd/rottinger/node60.html Is there a particular reason for this choice? The references in OP use Lagrange polynomial...