005
16.09.2004, 20:05 Uhr
Bruder Leif
dances with systems
(Operator)
|
Moin!
Nur zur Ergänzung noch ein bißchen Mathematik...
| C++: |
#include <stdio.h>
main()
{
int iCounter, iVorzeichen;
double fPi, fOldValue, fFaktor;
double fInterpol, fOldInterpol;
double fdoubleInterpol, fOlddoubleInterpol;
double fTripleInterpol, fOldTripleInterpol;
double fQuadInterpol, fOldQuadInterpol;
double fQuintInterpol, fOldQuintInterpol;
double fHexInterpol, fOldHexInterpol;
double fHeptInterpol, fOldHeptInterpol;
double fOktInterpol, fOldOktInterpol;
double fNonInterpol, fOldNonInterpol;
double fDekInterpol, fOldDekInterpol;
double fErgebnis;
printf( "Berechne PI, bitte warten... " );
fPi = fInterpol = fdoubleInterpol = 0;
fTripleInterpol = fQuadInterpol = 0;
fQuintInterpol = fHexInterpol = 0;
fHeptInterpol = fOktInterpol = 0;
fNonInterpol = fDekInterpol = 0;
fFaktor = iVorzeichen = 1;
for( iCounter=0; iCounter<=40; iCounter++ )
{
fOldDekInterpol = fDekInterpol;
fOldNonInterpol = fNonInterpol;
fOldOktInterpol = fOktInterpol;
fOldHeptInterpol = fHeptInterpol;
fOldHexInterpol = fHexInterpol;
fOldQuintInterpol = fQuintInterpol;
fOldQuadInterpol = fQuadInterpol;
fOldTripleInterpol = fTripleInterpol;
fOlddoubleInterpol = fdoubleInterpol;
fOldInterpol = fInterpol;
fOldValue = fPi;
fPi = fPi + iVorzeichen * ( 1 / fFaktor );
fFaktor += 2;
if( iVorzeichen == 1 ) iVorzeichen = -1;
else iVorzeichen = 1;
fInterpol = (fPi+fOldValue)*2;
fdoubleInterpol = (fInterpol+fOldInterpol)/2;
fTripleInterpol = (fdoubleInterpol+fOlddoubleInterpol)/2;
fQuadInterpol = (fTripleInterpol+fOldTripleInterpol)/2;
fQuintInterpol = (fQuadInterpol+fOldQuadInterpol)/2;
fHexInterpol = (fQuintInterpol+fOldQuintInterpol)/2;
fHeptInterpol = (fHexInterpol+fOldHexInterpol)/2;
fOktInterpol = (fHeptInterpol+fOldHeptInterpol)/2;
fNonInterpol = (fOktInterpol+fOldOktInterpol)/2;
fDekInterpol = (fNonInterpol+fOldNonInterpol)/2;
}
fErgebnis = (fDekInterpol+fOldDekInterpol)/2;
printf( "\rPi ist ungef„hr %.14lf\n", fErgebnis );
}
|
--
Mit 40 Fieber sitzt man nicht mehr vor dem PC.
Man liegt im Bett.
Mit dem Notebook. |
Anders ist es natürlich, wenn man eine "astronomische Genauigkeit" braucht, also wenigstens 30 Stellen nach dem Komma. Aber dann arbeitet man eh mit speziellen Bibliotheken, und nicht mit der Arithmetik des Compilers...
