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:
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:
Bonjour.
Quelques erreurs de frappes :
Faire :
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; .....
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)
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
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
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...
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.j'ai envie
Jean-Marc Blanc
Vous avez un bloqueur de publicités installé.
Le Club Developpez.com n'affiche que des publicités IT, discrètes et non intrusives.
Afin que nous puissions continuer à vous fournir gratuitement du contenu de qualité, merci de nous soutenir en désactivant votre bloqueur de publicités sur Developpez.com.
Partager