Re: failure to Integrate in orthogonal polynomials
- To: mathgroup at smc.vnet.net
- Subject: [mg67918] Re: failure to Integrate in orthogonal polynomials
- From: Roger Bagula <rlbagula at sbcglobal.net>
- Date: Fri, 14 Jul 2006 02:11:09 -0400 (EDT)
- References: <e8nt4n$k7c$1@smc.vnet.net>
- Sender: owner-wri-mathgroup at wolfram.com
Based on my experience that a unit sphere in an so(3) matrix gives a square that gives a projective plane I did in the off diagonal elements: M = {{0, LegendreP[n, x], -LegendreP[n + 1, x]}, {-LegendreP[n, x], 0, LegendreP[n + 2, x]}, {LegendreP[n + 1, x], -LegendreP[n + 2, x], 0}} M.M Unnormalized the functions tend to grow larger as they are projected. When presented in this way the result won't integrate: norm[n_]=Sqrt[1/Integrate[LegendreP[n,x]^2,{x,-1,1}]] g=Table[ParametricPlot3D[{norm[n+1]*norm[n+2]*LegendreP[n+1,x]*LegendreP[n+2,\ x],norm[n]*norm[n+2]*LegendreP[n,x]*LegendreP[n+2,x],norm[n]*norm[n+1]*\ LegendreP[n,x]*LegendreP[n+1,x]},{x,-1,1},PlotPoints\[Rule]1000],{n,0,10}] Show[g] What works to give an idea of what these functions look like is: Clear[norm, p0, p, g] p0 = Table[LegendreP[n, x], {n, 0, 12}] norm = Table[(1/Integrate[p0[[n]]*p0[[n]], {x, -1, 1}])^(1/2), {n, 1, 12}] p = Table[p0[[n]]*norm[[n]], {n, 1, 12}] g = Table[ParametricPlot3D[{p[[n + 1]]*p[[n + 2]], p[[n]]*p[[n + 2]], p[[n]]*p[[n + 1]]}, {x, -1, 1}, PlotPoints -> 1000], {n, 1, 10}] Show[g] Show[g, ViewPoint -> {3.354, -0.048, 0.447}] Show[g, ViewPoint -> {0.045, 0.000, 3.383}] The outcome seems to be two incoming and outgoing cyclinders. That meet in a crosscap like figure? They resemble Lerenz and Rössler chaotic trajectories. I think these functions give an interesting set of new functions in 3D. > >