Page 1 of 1

[WetProg] Snelheid Oefening 6 2e Zittijd (trapeziumregel)

Posted: Wed Aug 30, 2006 12:12 pm
by Valthonis
Heeft iemand oefening 6 al gemaakt, waarbij Pi moet worden berekend via de integraal op interval [-1,1] van de functie

Je deelt het interval op in k deelintervallen en berekent zo elke keer de T(k) zijnde de berekende benadering, je voert dan k op tot dus de relatieve fout tussen 2 opeenvolgende stappen vrij tot zeer klein is.

Nu heb ik hem een hele tijd laten lopen tot k = +- 70000 en dan heb ik voor pi =3.14159247644715 wat dus correct is tot op 6 cijfers na de komma en mijn relatieve fout naar de volgende stap is -1,20... e-12 en zou kleiner moeten worden dan 9,09... -13 dus we zijn er wel bijna maar is dit wel correct?

Posted: Thu Aug 31, 2006 9:18 pm
by Sond
Bij mij stopt hij ergens in de 30 of de 40 (weet nie meer correct en geen zin om te rebooten om te checken) dus ik denk dat ge iets verkeerd hebt gedaan...

Ik heb da gedaan in c++ met long doubles
Schijnt da hoge precisies (net zoals in oefening 2) hier wel handig zijn....

Posted: Sat Sep 02, 2006 8:15 am
by Valthonis
Ik denk dak al weet wa dak fout heb gedaan. Het zijn 2^k deelintervallen en niet in k intervallen. Dat kan natuurlijk wel de fout zijn want dan gaat dat wel iets sneller gaan.

Woeps :lol:

Posted: Sat Sep 02, 2006 10:11 am
by Valthonis
ok 2 vectoren van elks enkele miljoenen doubles is dus duidelijk teveel voor matlab.
Dan maar even zelf implementeren en dan haal ik het gevraagde na 29 iteraties. Eindelijk af.

Nu het verslag nog en al dik geen goesting meer :lol: