%interp_poly
choix_menu = MENU('choisissez une methode d interpolation ','1) polynome de lagrange','2)polynome de newtton en cas ou les points x sont quelconques','3)polynome de newtton dans le cas ou les points sont equidistants');
switch choix_menu
case 1
disp ('vous avez choisi la methode de LAGRANGE ')
a=input('donnezla valuer de a');
x=input('donnez les valeur de x ');
y=input('donnez les valeur de y');
fonc_lagrange(a,x,y) ;
case 2
disp ('vous avez choisi la methode de NEWTTON ')
a=input('donnezla valuer de a');
x=input('donnez les valeur de x ');
y=input('donnez les valeur de y');
n=input('donnez les valeur de n');
[DD]=calc_table_dd(x,y,n);
z=diag(DD);
p=y(1);
for k=1:1:n-1
[s]=produit(a,x,k);
p=p+(z(k)*s);
disp(['Pn(',num2str(a),')=',num2str(p)])
end
d=fopen('newton_quel.m','w');
fprintf(d,'%d',p);
fclose(d)
case 3
disp('vous avez choisi le polynome de NEWTTON en cas ou les points x sont quelconques ')
disp('choisissez lun des deux methodes 1 pour le choix progressif et 2 pour le choix regressif')
choix=input('votre choix ');
x=input('donnez les valeur de x ');
y=input('donnez les valeur de y');
n=input('donnez la valeur de n');
[DF]=calc_table_df(n,x,y,choix);
p=y(1);
for k= 1:1:n-1
[s]=produit_progressif(a,k);
p=p+(s*z(k))/(factorial(k));
disp(['Pn(',num2str(a),')=',num2str(p)])
end
d=fopen('newton_equid.m','w');
fclose(d)
end