- ## Compute the Legendre polynomial over [-1,1] and mangle it.
- P = legendre_p(n)
- P_tilde = @(x) P( (2/(b-a))*x + 1 - (2*b)/(b-a) )
+ ## Compute the Legendre polynomial over [-1,1] and mangle it to fit
+ ## the interval [a,b].
+ P = legendre_p(n);
+ P_tilde = @(x) P( (2/(b-a)).*x + 1 - (2*b)/(b-a) );