Re: Simpson vs. Niski Cotes - Grupy dyskusyjne w eGospodarka.pl

eGospodarka.pl › Grupy › pl.comp.programming › Simpson vs. Niski Cotes › Re: Simpson vs. Niski Cotes

Path: news-archive.icm.edu.pl!news.icm.edu.pl!newsfeed.neostrada.pl!unt-exc-02.news.n
eostrada.pl!unt-spo-a-02.news.neostrada.pl!news.neostrada.pl.POSTED!not-for-mai
l
From: "slawek" <s...@h...pl>
Newsgroups: pl.comp.programming
References: <509ee300$0$26682$65785112@news.neostrada.pl>
<k7olf5$rpm$1@news.task.gda.pl> <k7rnav$8eq$1@node2.news.atman.pl>
<k7t76g$3f0$1@zeus.man.szczecin.pl>
<4...@g...com>
<k7tbmv$6se$1@zeus.man.szczecin.pl>
<4...@g...com>
In-Reply-To: <4...@g...com>
Subject: Re: Simpson vs. Niski Cotes
Date: Wed, 14 Nov 2012 13:09:06 +0100
MIME-Version: 1.0
Content-Type: text/plain; format=flowed; charset="iso-8859-1"; reply-type=original
Content-Transfer-Encoding: 8bit
X-Priority: 3
X-MSMail-Priority: Normal
Importance: Normal
X-Newsreader: Microsoft Windows Live Mail 14.0.8117.416
X-MimeOLE: Produced By Microsoft MimeOLE V14.0.8117.416
Lines: 33
Message-ID: <50a389de$0$1222$65785112@news.neostrada.pl>
Organization: Telekomunikacja Polska
NNTP-Posting-Host: 62.69.202.124
X-Trace: 1352894942 unt-rea-b-01.news.neostrada.pl 1222 62.69.202.124:60041
X-Complaints-To: a...@n...neostrada.pl
Xref: news-archive.icm.edu.pl pl.comp.programming:200847
[ ukryj nagłówki ]

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 ;)