-
Path: news-archive.icm.edu.pl!agh.edu.pl!news.agh.edu.pl!newsfeed2.atman.pl!newsfeed.
atman.pl!newsfeed.neostrada.pl!unt-exc-01.news.neostrada.pl!unt-spo-a-01.news.n
eostrada.pl!news.neostrada.pl.POSTED!not-for-mail
From: "slawek" <s...@h...pl>
Newsgroups: pl.comp.programming
Subject: Simpson vs. Niski Cotes
Date: Sun, 11 Nov 2012 00:28:01 +0100
MIME-Version: 1.0
Content-Type: text/plain; format=flowed; charset="iso-8859-2"; 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: 79
Message-ID: <509ee300$0$26682$65785112@news.neostrada.pl>
Organization: Telekomunikacja Polska
NNTP-Posting-Host: 62.69.202.124
X-Trace: 1352590080 unt-rea-a-01.news.neostrada.pl 26682 62.69.202.124:65367
X-Complaints-To: a...@n...neostrada.pl
Xref: news-archive.icm.edu.pl pl.comp.programming:200706
[ ukryj nagłówki ]W kontekście rozmaitych forumowych "specjalistów" od metod numerycznych -
zabawne, jak tym razem będą wyjaśniać trywialny w istocie rezultat...
exact 0.5022749400837604
tapez 0.5022749194207212
simpson 0.5022738634614975
/*
* tq.c
* 2012 (C) slawek
*/
#include <stdio.h>
#include <math.h>
#define NPTS 10000 /* NPTS must be even! */
double x[NPTS+1];
double y[NPTS+1];
double f(double x)
{
return sin(x)*exp(-x);
}
double integrate_exact(double a, double b)
{
return ((cos(a) + sin(a))*exp(-a) - (cos(b) + sin(b))*exp(-b))/2.0;
}
double generate_points(double x[], double y[], double a, double b, int n)
{
int i;
for(i = 1; i <= n; i++)
{
x[i] = a + (i-1) * (b - a)/(n - 1);
y[i] = f(x[i]);
}
}
double integrate_trapez(double x[], double y[], int n)
{
int i;
double s = 0.;
for(i = 2; i <= n; i++)
s += (y[i] + y[i-1]) * (x[i] - x[i-1]) * 0.5; // an naive approach
return s;
}
double integrate_simpson(double x[], double y[], int n)
{
int i;
double odd = 0., even = 0.;
for(i = 1 ; i < n ; i+=2) odd += y[i];
for(i = 2 ; i <= n ; i+=2) even += y[i];
return (2.0*odd-y[1]-y[n] + 4.0*even)*(x[3]-x[1])/6.0;
}
int main(int argc, char* argv[])
{
const double a = 0.;
const double b = 5.;
generate_points(x,y,a,b,NPTS);
printf("exact %.16lg\n", integrate_exact(a,b));
printf("tapez %.16lg\n", integrate_trapez(x,y,NPTS));
printf("simpson %.16lg\n", integrate_simpson(x,y,NPTS));
puts("press enter...");
getchar();
return 0;
}
Następne wpisy z tego wątku
- 11.11.12 14:33 slawek
- 11.11.12 16:42 e...@g...com
- 11.11.12 16:45 kenobi
- 11.11.12 17:29 s
- 11.11.12 17:58 Baranosiu
- 11.11.12 18:33 AK
- 11.11.12 18:44 AK
- 12.11.12 05:11 slawek
- 12.11.12 06:01 slawek
- 12.11.12 06:08 slawek
- 12.11.12 10:46 AK
- 12.11.12 10:48 AK
- 12.11.12 12:53 slawek
- 12.11.12 13:07 slawek
- 12.11.12 13:32 Michoo
Najnowsze wątki z tej grupy
- Alg. kompresji LZW
- Popr. 14. Nauka i Praca Programisty C++ w III Rzeczy (pospolitej)
- Arch. Prog. Nieuprzywilejowanych w pełnej wer. na nowej s. WWW energokod.pl
- 7. Raport Totaliztyczny: Sprawa Qt Group wer. 424
- TCL - problem z escape ostatniego \ w nawiasach {}
- Nauka i Praca Programisty C++ w III Rzeczy (pospolitej)
- testy-wyd-sort - Podsumowanie
- Tworzenie Programów Nieuprzywilejowanych Opartych Na Wtyczkach
- Do czego nadaje się QDockWidget z bibl. Qt?
- Bibl. Qt jest sztucznie ograniczona - jest nieprzydatna do celów komercyjnych
- Co sciaga kretynow
- AEiC 2024 - Ada-Europe conference - Deadlines Approaching
- Jakie są dobre zasady programowania programów opartych na wtyczkach?
- sprawdzanie słów kluczowych dot. zła
- Re: W czym sie teraz pisze programy??
Najnowsze wątki
- 2025-02-12 Dziwne zachowanie magistrali adresowej w 8085
- 2025-02-11 Mini pecet
- 2025-02-10 Spalił się spaliniak
- 2025-02-10 zarowka wifi - z sensowna apka lub lepiej albo lokalnie lub przez web. I zeby harmonogram miala
- 2025-02-10 Chrzanów => Programista NodeJS <=
- 2025-02-10 Kraków => DevOps Engineer (Junior or Regular level) <=
- 2025-02-10 Dlaczego takie preferencje banków?
- 2025-02-10 Białystok => iOS Developer (Swift) <=
- 2025-02-10 Mińsk Mazowiecki => Team Lead / Tribe Lead FrontEnd <=
- 2025-02-10 Białystok => System Architect (Java background) <=
- 2025-02-10 Współczesne mierniki zniekształceń nieliniowych THD audio, produkują jakieś?
- 2025-02-10 Szczecin => Senior Field Sales (system ERP) <=
- 2025-02-10 Gliwice => Business Development Manager - Dział Sieci i Bezpieczeńst
- 2025-02-10 Chrzanów => Specjalista ds. public relations <=
- 2025-02-10 Chrzanów => NodeJS Developer <=