Pi Berechnen

  • Pi Berechnen

    Hallo Liebes Forum!

    Ich brauche dringend hilfe^^ xD mal wieder :<

    also..
    Ich hab einen Algorithmus gefunden um Pi auf die Nte Stelle auszurechnen wie schon in diesem Thread gesagt..
    Das Ergebnis jeder Nten Stelle wird im Hexadezimal angegeben, also Basis 16.

    Der Algorithmus ist die BBP Reihe; Umgestellt.
    Hab dies hier her.

    Ich habe mich mal dran Probiert dies zu Programmieren allerdings bekomm ich nicht das richtige Ergebnis(kann man sich denken xD).

    Ich kommentier leider nicht also bitte nicht meckern deswegen^^
    btw bin C++ boon also anfänger :D

    Mein Code bis Jetzt:

    Quellcode

    1. #include <iostream>
    2. #include <string>
    3. #include <conio.h>
    4. #include <cmath>
    6. using namespace std;
    8. #define Unendlich 999999
    10. double D_Mod(double a, double b)
    11. {
    12. int result = static_cast<int>( a / b);
    13. return a - static_cast<double>(result) * b;
    14. }
    17. int main()
    18. {
    19. //Variablen deklaration
    20. int Anzahl = 0;
    21. string Pi = "3.";
    22. char Addieren = 'z';
    25. double Summe1 = 0, Summe2 = 0, Summe3 = 0, Summe4 = 0;
    27. double dZahl1, dZahl2, k1, n1, Zwischensumme;
    28. int temp, iErgebnis;
    30. cout << "Wie viele Stellen sollen Berechnet werden?" << endl;
    31. cin >> Anzahl;
    33. for(int n=0;n<Anzahl;n++)
    34. {
    35. Summe1 = 0;
    36. Summe2 = 0;
    37. Summe3 = 0;
    38. Summe4 = 0;
    39. dZahl1 = 0;
    40. dZahl2 = 0;
    41. Zwischensumme = 0;
    42. temp = 0;
    43. iErgebnis = 0;
    45. //Summe1
    46. for(int k=0;k<n;k++)
    47. {
    48. k1 = k;
    49. n1 = n;
    50. temp = pow(16,n1-k1);
    51. dZahl1 = temp % (8*k+1);
    52. dZahl2 = 8*k+1;
    53. Summe1 += dZahl1 / dZahl2;
    54. }
    56. for(int k=n+1;k<Unendlich;k++)
    57. {
    58. k1 = k;
    59. n1 = n;
    60. dZahl1 = pow(16,n1-k1);
    61. dZahl2 = 8*k+1;
    63. Summe1 += dZahl1 / dZahl2;
    64. }
    66. Summe1 *= 4;
    69. //Summe2
    70. for(int k=0;k<n;k++)
    71. {
    72. k1 = k;
    73. n1 = n;
    74. temp = pow(16,n1-k1);
    75. dZahl1 = temp % (8*k+4);
    76. dZahl2 = 8*k+4;
    77. Summe2 += dZahl1 / dZahl2;
    78. }
    80. for(int k=n+1;k<Unendlich;k++)
    81. {
    82. k1 = k;
    83. n1 = n;
    84. dZahl1 = pow(16,n1-k1);
    85. dZahl2 = 8*k+4;
    87. Summe2 += dZahl1 / dZahl2;
    88. }
    90. Summe2 *= 2;
    92. //Summe3
    93. for(int k=0;k<n;k++)
    94. {
    95. k1 = k;
    96. n1 = n;
    97. temp = pow(16,n1-k1);
    98. dZahl1 = temp % (8*k+5);
    99. dZahl2 = 8*k+5;
    100. Summe3 += dZahl1 / dZahl2;
    101. }
    103. for(int k=n+1;k<Unendlich;k++)
    104. {
    105. k1 = k;
    106. n1 = n;
    107. dZahl1 = pow(16,n1-k1);
    108. dZahl2 = 8*k+5;
    110. Summe3 += dZahl1 / dZahl2;
    111. }
    113. //Summe4
    114. for(int k=0;k<n;k++)
    115. {
    116. k1 = k;
    117. n1 = n;
    118. temp = pow(16,n1-k1);
    119. dZahl1 = temp % (8*k+6);
    120. dZahl2 = 8*k+6;
    121. Summe4 += dZahl1 / dZahl2;
    122. }
    124. for(int k=n+1;k<Unendlich;k++)
    125. {
    126. k1 = k;
    127. n1 = n;
    128. dZahl1 = pow(16,n1-k1);
    129. dZahl2 = 8*k+6;
    131. Summe4 += dZahl1 / dZahl2;
    132. }
    134. //Berechnung
    135. Zwischensumme = Summe1 - Summe2 - Summe3 - Summe4;
    136. Zwischensumme = D_Mod(Zwischensumme,1);
    137. Zwischensumme *= 16;
    138. iErgebnis = Zwischensumme;
    140. switch(iErgebnis)
    141. {
    142. case 0:
    143. Addieren = '0';
    144. break;
    145. case 1:
    146. Addieren = '1';
    147. break;
    148. case 2:
    149. Addieren = '2';
    150. break;
    151. case 3:
    152. Addieren = '3';
    153. break;
    154. case 4:
    155. Addieren = '4';
    156. break;
    157. case 5:
    158. Addieren = '5';
    159. break;
    160. case 6:
    161. Addieren = '6';
    162. break;
    163. case 7:
    164. Addieren = '7';
    165. break;
    166. case 8:
    167. Addieren = '8';
    168. break;
    169. case 9:
    170. Addieren = '9';
    171. break;
    172. case 10:
    173. Addieren = 'A';
    174. break;
    175. case 11:
    176. Addieren = 'B';
    177. break;
    178. case 12:
    179. Addieren = 'C';
    180. break;
    181. case 13:
    182. Addieren = 'D';
    183. break;
    184. case 14:
    185. Addieren = 'E';
    186. break;
    187. case 15:
    188. Addieren = 'F';
    189. break;
    190. }
    192. Pi += Addieren;
    193. }
    195. cout << endl;
    196. cout << Pi;
    201. _getch();
    202. return 0;
    203. }
    Alles anzeigen


    Findet ihr meinen Fehler / meine Fehler?

    DANKE FÜR EURE HILFE! :D

  • Ich such jetzt zwar nicht deinen Fehler, kann dir aber sagen, wie du ihn wahrscheinlich eher findest ....
    1. Benutz Methoden, schreib nicht alles in die Main-Methode rein
    2. Versuch dich bei den Variablen Namen an ein Schema zu halten, nicht teils ungarische Notation, teils gar nichts
    3. Summe1 - Summe4, warum kein Feld ?
    3. Ein ganz klein wenig mehr Kommentare ...

    So wird es sich imho nämlich ziemlich sicher fast keiner Antun deinen Code zu entschlüsseln ...


    MrSpider

  • Normalerweise lager ich auch in funktionen aus aber hier das ist eher noch pseudo code und mit den variablen ist mir klar^^ mach ich ehrlich gesagt normalerweise auch anders..

    Was ich ja So komisch finde ist, wenn ich mir den Innentel mit den Schleifen rauskopiere und mache das so das ich für n 5 sage berechnet er mir die Stelle auch richtig, erhöhe ich aber n um 1 ists falsch :O .

  • Ist gar nicht so unwahrscheinlich - du solltest das Rechnen mit double forcieren indem du konsequent bspsw. 8.0 anstatt 8 schreibst. Also immer double-Werte verwenden, auch bei Zahlen die im Code stehen. Ansonsten solltest du dir die Operatorenprioritätentabelle und die implizite Konvertierungstabelle neben den PC legen oder auswendig lernen. Ich glaube das möchtest du nicht.
    ~ mfg SeBa

    Ich beantworte keine PMs zu Computer-/Programmierproblemen. Bitte wendet euch an das entsprechende Forum.

    [Blockierte Grafik: http://i.creativecommons.org/l/by-sa/3.0/80x15.png]