j'ai bricoler cette algorythme avec la méthode des trapezes; f(x) une fonction entrée par l'utilisateur sur un intervalle [a;b] découpé en n partie égale(plus n est grand plus l'integrale de f(x) sur [a;b] est précise)
vous en pensez quoi ?
- name:Méthode integration trapeze#stepline
- set:#vars:a;b;n;m;k;i;h#fct:f(x)#changefxt:stepl
ine;stl#changefxt:scan;sc#stl
- start:#stl
- sc:a#sc:b#sc:n#sc:f(x)#0;m#0;k#0;i#(b-a)/n;h#stl
- while:m<(b-a)#stl
- i+0,5*h(f(a+k*h)+f(a+(k+1)h));i#stl
- k+1;k#stl
- m+h;m#stl
- whileend#stl
- end:#println:i#stl
- stop