Discussion :

[bsangoku]

Bonjour,

Je suis débutant en Matlab (je n'ai qu'une légère expérience d'un mois) et je cherche à faire la méthode de Newton Cotes. Je me suis aidé d'un livre pour le faire:

[phryte]

Bonjour.
Quelques erreurs de frappes :

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 int= newtcot(a,b,n,fun)
%NEWTCOT Formules fermées de Newton
%INT=NEWTCOT(A,B,N,FUN) calcule une approx de l'intégrale de la
% fonction FUN sur [A,B] par la formule fermée de Newton-Cotes à N noeuds
%FUN accepte en entrée un valeur réel x et renvoie un vecteur réel.
h=(b-a)/n;
n2=fix(n/2)
if n>6, error('n vaut au plus 6'); end
a03=1/3; 
a08=1/8; 
a45=1/45; 
a288=1/288; 
a140=1/140;
.....

Faire :

Code :

Sélectionner tout - Visualiser dans une fenêtre à part

1
2
3
4
5
6
clear
a=-1;
b=1;
n=2;
fun='2*x+1';
res=newtcot(a,b,n,fun)

[bsangoku]

SUPER!!!! Merci de votre réponse: la fonction marche à merveille...

Cependant, je voulais vous demander si vous connaissiez un site dans lequel on peut trouver la valeur des poids des formules de newton cotes pour des valeurs n compris entre 7 et 12 svp...

Je suis en train de chercher et je n'arrive pas à en trouver.
Merci d'avance

[FR119492]

Salut!
Tu trouveras une réponse partielle à ta question dans l'ouvrage de M. Abramowitz et I. Stegun: "Handbook of mathematical functions".
Jean-Marc Blanc

[bsangoku]

Merci de votre réponse mais ceci ne m'aide pas vraiment
En réalité, j'ai envie de compléter mon programme pour N=7,8,9,10 et pour cela j'ai besoin de connaitre les poids correspondants...

[FR119492]

j'ai envie

Dans la vie, on ne fait pas toujours ce dont on a envie. Si on arrive à faire ce qu'on doit, c'est déjà pas mal.
Jean-Marc Blanc