Du mal àdémarrer ?
23
✞
☎
function y=fonctj(x,nomf,j);
y=feval(nomf,x).*x.^(j+1/2);
✝
✆
✞
☎
function b=sndmemb(nomf,n);
for i=1:n+1
b(i)=quadl(’fonctj’,0,1,[] ,[ ],nomf,i-1);
end;
✝
✆
✞
☎
>> sndmemb(’f1’,4)
ans =
0.4374
0.2416
0.1521
0.1041
0.0755
✝
✆
4. Une fois la solution a déterminée, le polynôme se calcule àl’aide de polyval (attention à
l’ordre des coefficients). Calculer a et tracer le graphe de f puis celui du polynôme p n .Sauvegarde
dans approx2.m.T est avec n = 4, f1.
✞
☎
>> approx2
Coefficients du polynome
a=
0.0022
3.0765
0.5720
-7.2921
3.6428
✝
✆
5. Faire varier n .Changer la fonction. Que se passe-t-il quand n augmente ( n 20) ?
Exemples f 1(x ) = sin p x ; f 2(x ) = sin 4 p x ; f 3(x ) = | x − 1 / 2 | ; n = 2 , 3 , 5 , 10, 15, 20...
DU MAL ÀDÉMARRER
?
3.1 Premiers calculs
1. Ànoter que A
− 1
est donnée dans l’exercice 2.1.
2. Si U vecteur propre, écrire l U = AU et majorer la norme.
3.2C onditionnement et erreur
Utiliser la forme f ( x 0 ,...,x n ) 0.
 Dunod –Laphotocopie non autorisée est un délit