Updated for new ApproxFun

Author: dlfivefifty

src/InvariantEnsembles.jl

@@ -33,7 +33,7 @@ export InvariantEnsemble
 function orthonormalpolynomials(μ0,α::Vector,β::Vector,d)
     n=length(α)+1
 
-    p=Array(IFun{Float64,Interval{Float64}},n)
+    p=Array(IFun{Float64},n)
     p[1] = Fun([μ0],d)
     p[2] = multiplybyx(p[1])./β[1] - p[1].*α[1]./β[1]
     for k = 3:n
@@ -257,7 +257,7 @@ function samplespectra(q::Array{Float64,2},d,plan::Function,pts)
 end
 
 
-function iekernel(p::Array{IFun{Float64,Interval{Float64}}})
+function iekernel{D}(p::Array{IFun{Float64,D}})
     ret = 0.*p[1]
     for i = 1:length(p)
         ret += fasttimes(p[i],p[i])
@@ -266,7 +266,7 @@ function iekernel(p::Array{IFun{Float64,Interval{Float64}}})
     ret
 end
 
-function samplespectra(p::Array{IFun{Float64,Interval{Float64}}})
+function samplespectra{D}(p::Array{IFun{Float64,D}})
     n = length(p)
     r=sample(iekernel(p)/n)
     if n==1