Uzytkownik "Roman W" <r...@g...com> napisal w wiadomosci
grup dyskusyjnych:43ee067c-4fa1-4d06-8030-8b1f97cd6bfa@go
oglegroups.com...
> OK, point taken: sprawdzilem wersje z tablica (kod ponizej). Wynik jest
> nadal
> 0.5022749400837603, czyli praktycznie tyle samo, co analityczny rezultat.

Poprawione, z lepszym sumowaniem jest róznica równa, sic!, "zlemu
epsilonowi" w porównaniu z wynikiem "dokladnym" (sinusy itd. tez nie sa
idealnie dokladne). Lepiej widac przy kroju pisma o stalej szerokosci
czcionki.

NPTS = 50001
a = 0
b = 7
f(x) = sin(x)*exp(-x)

exact 0.499356717949745400
NC1 0.499356716316665280 -3.27e-003 ppm
NC1-K 0.499356716316659790 -3.27e-003 ppm
NC2 0.499356717949743840 -3.11e-009 ppm
NC2-K 0.499356717949745340 -1.11e-010 ppm
SPLI3 0.499356719582675580 +3.27e-003 ppm

NCn - Newton-Cotes n-tego stopnia (1-trapez, 2-Simpson)

-1.11e-10 ppm to jest -1.11e-10 * 1e-6 czyli -1.11e-16, a to jest 1/2 z
DBL_EPSILON - czyli "zly epsilon" - samo w sobie ciekawe.


Cóz, aby dalej potrollowac poszukam jakiejs wrednej funkcji ;)