Integrazione sui nodi di gauss-chebychev

[caty89]

salve a tutti! ho un problema: devo scrivere un codice in matlab per integrare utilizzando i nodi di gauss-chebychev. ma non ho idea di come fare. qualcuno può darmi qualche consiglio su come procedere o qualche documento da leggere? grazie anticipatamente!

[caty89]

"TeM":Scrivere una funzione Matlab che data una funzione \(f\), un intervallo \([a, \; b]\) e un numero naturale \(N\) calcoli:
(i) il polinomio \(P_N\) che interpola \(f\) nei \(N+1\) nodi di Chebyshev-Gauss-Lobatto dell'intervallo \([a, \; b]\);
(ii) l'integrale \(I_N\) tra \(a\) e \(b\) di \(P_N\).

function [ P, I ] = quadcheb( f, a, b, N )
hch = - cos( pi * [ 0 : N ] / N ) ;
xch = ( b + a ) / 2 + ( b - a ) / 2 * hch ;
fxch = feval( f, xch ) ;
P = polyfit( xch, fxch, N ) ;
Q = polyint( P ) ;
I = polyval( Q, b ) - polyval( Q, a ) ;
return

Spero sia abbastanza chiaro

grazie della risposta rapidissima e della disponibilità! è chiarissimo!
io non riuscivo a capire cosa fare dopo aver costruito i nodi e aver calcolato i valori della f (da polyfit in giù non ci ero proprio arrivata). adesso provo a fare come mi hai consigliato. ancora grazie mille!