scilab quadrature

xavierdestev

Bonjour,

J'ai codé cette méthode quadrature mais mon problème est qu'il faut que je l'utilise sur des fonctions $\phi$ qui dépendent aussi de $i$ et de $N$. (c'est-à-dire d'une discrétisation régulière)
Ce sont les fonctions chapeaux de la méthode des éléments finis $1D$.

Comment puis-je utiliser cette méthode sur $\phi(x,I,N)$

function I=integrale2(a,b,r)
~~~H1=5/9;
~~~H2=8/9;
~~~H3=5/9;
~~~q=sqrt(3/5);
~~~e1=a+(1-q)*(b-a)/2;
~~~e2=(a+b)/2;
~~~e3=a+(1+q)*(b-a)/2;
~~~I=((b-a)/2)*(H1*r(e1)+H2*r(e2)+H3*r(e3));
endfunction

Merci d'avance.