Update multivariatelanczos.jl

Author: dlfivefifty

examples/multivariatelanczos.jl

@@ -2,9 +2,14 @@ using MultivariateOrthogonalPolynomials, ClassicalOrthogonalPolynomials, Test
 
 P = Legendre()
 P² = RectPolynomial(P,P)
+p₀ = expand(P², 𝐱 -> 1)
+𝐱 = axes(P²,1)
+x,y = first.(𝐱),last.(𝐱)
+w = P²/P²\ (x-y).^2
+
+w .* p₀
 
 
-p₀ = expand(P², 𝐱 -> 1)
 sum(p₀)
 using LazyBandedMatrices
 KronTrav(sum(P;dims=1), sum(P;dims=1))