salut ,
J'ai besoin un coup de main
Je dois de calculer l'interpolation de Newton de la fonction exponentielle sur l'intervalle [-1,1] on commence par remplir le tableau de différences devisé
mon essai:
Code : Sélectionner tout - Visualiser dans une fenêtre à part
1
2
3
4
5
6
7
8
9
10
11
12
13
14 function [DD]=calc_table_dd(x,y,n); xx=linspace(sizeof(x)); z=diag(DD); y=exp(x); p=y(1); for k=1:1:n-1 [s]=produit(a,xx,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)
Partager