I need to find root for Legendre polynomials in terms of the RadauIIA method
using ClassicalOrthogonalPolynomials, NLsolve function P(x, n) return legendrep(n, x) end function radc(x, s) v = 2 * x - 1 return P(v, s) - P(v, s-1) end s = 3 k = zeros(s) function radc!(fk, k) fk .= radc(k, s) end c = nlsolve(radc!, k)
nlsolve is difficult to use. I used this same method of nlsolve for something else and it worked, but not for the above.
Please help me.