[TPW] Méthode de Wallis

Version imprimable

Salut les gars ;
Bon voila j'ai un petit probléme dans mon code source :aie: !!de l'aide sera la bienvenue surtout que ma tete me fait trés mal :lol: et je vous remercie d'avance pour vos réponses !!;)
voila l'enoncé

Citation:

On se propose de determiner une valeur approchée de pi par la methode de Wallis definie par la formule suivante pi/2=(2/1)*(2/3)*(4/3)*(4/5)*(6/5)*(6/7)*(8/7)*(8/9)etc Ecrire un programme pascal qui permet de determiner et d'afficher une valeur approchée de pi .Le calcul s'arrete lorsque la difference entre 2 valeurs consecutifs de cette formule devient strictement < 10 puissance -6

Voila ce que j'ai pu faire en utilisant l'algorithme de récurrence
Code:

1 2 3 4 5 6 7 8 9 10 11 12 13 14 15 16 17 18 19 20 21 22 23 24 25 26 27 28 29 30 31 32 33 34 35 36 37 38 39 program exercice ; uses wincrt ; var x,s :real ; function calcul : real ; var a,b,c : integer ; s1 : real ; begin a:=2 ; b:=1 ; s:=0 ; c:=0 ; repeat c:=c+1 ; s1:=s ; s:=s+(a/b) ; if c mod 2 <> 0 then begin b:=b+2 ; end else if c mod 2 = 0 then begin a:=a+2 ; end; until abs(s-s1) < 0.000001 ; calcul:=s*2 ; end; begin x:=calcul ; writeln(x); end.

Bonjour,

Ton code ne suit pas l'algorithme:
Code:

1 2 3 ... s:=s+(a/b) ; ...
devrait être
Code:

1 2 3 ... s:=s*(a/b) ; ...
Ce qui interdit d'initialiser s = 0.

Pour échapper à ça, tu peux initialiser s avec la valeur du 1er terme (=2/1 = 2),
et donc aussi changer la valeur d'origine de ton c, qui contrôle la manière d'incrémenter les variables a et b:
Code:

1 2 3 4 5 6 7 8 9 10 11 12 13 14 15 16 17 18 19 20 21 22 23 24 25 26 27 28 29 30 31 32 33 34 35 36 program exercice ; uses wincrt ; var x :real ; function calcul : real ; var a,b,c : integer; s,s1 : real ; begin a:=2 ; b:=3 ; s:= 2; c:=1 ; repeat c:=c+1 ; s1:=s ; s:=s*(a/b) ; if c mod 2 <> 0 then begin b:=b+2 ; end else if c mod 2 = 0 then begin a:=a+2 ; end; until (abs(s-s1) < 0.000001); calcul:=s*2 ; end; begin x:=calcul ; writeln(x); readln; end.
Tu auras noté au passage que j'ai supprimé la variable globale s, pas du tout justifiée, et que j'ai ajouté "readln;" en dernière ligne, pour éviter que la console se ferme sans que tu aies le temps de lire l'affichage.

22/01/2007, 12h39
WhiteTigerZ

salut , ;)
ça marche pas !! il n'affiche pas le Pi lol
Merci quand même
22/01/2007, 12h43
WhiteTigerZ

ça marche

Bon ça marche il faut modifier la condition

Code:

(abs(s-s1) < 0.000001 en (abs(s-s1) < 0.0001 {10 a la puissance -4 et non -6}

peut etre il s'agit d'un probléme dans l'enoncé de l'exercice
:)
22/01/2007, 12h52
thewho

Bonjour,

Chez moi, ça marche avec (abs(s-s1) < 0.000001)

Mais ça peut venir du compilateur, en particulier de la précision qu'il utilise pour le type real.
22/01/2007, 13h08
WhiteTigerZ

Citation:

Envoyé par thewho

Bonjour,

Chez moi, ça marche avec (abs(s-s1) < 0.000001)

Mais ça peut venir du compilateur, en particulier de la précision qu'il utilise pour le type real.

j'utilise TP for Windows mais il m'affiche 0.00000000 quand je modifie la difference a 10 puissance -6 !!
qu'est ce que je dois faire pour modifier la précision du type réel ??
22/01/2007, 14h31
thewho

Bonjour,

Citation:

Envoyé par WhiteTigerZ

j'utilise TP for Windows mais il m'affiche 0.00000000 quand je modifie la difference a 10 puissance -6 !!
qu'est ce que je dois faire pour modifier la précision du type réel ??

Il y a si longtemps que je n'utilise plus ce truc !

Vas voir dans l'aide les types numériques supportés, et prend celui qui correspond à tes besoins.

J'ai vaguement le souvenir que le type "double" existe.
Si c'est bien le cas, choisis le, car il présente tous les avantages par rapport à "real" :
- meilleure précision (8 octets au lieu de 6)
- beaucoup plus rapide pour les calculs, car il utilise le coprocesseur, ce qui n'est pas le cas avec le type "real"
23/01/2007, 16h25
WhiteTigerZ

j'ai essayé le type double mais il parait qu'il n'a rien fait
toujours le même probléme !! a propos j'ai pas trouvé la page des types numériques lol ! sinon qu'est ce que tu utilises comme compilateur :)
23/01/2007, 17h29
thewho

Bonjour,

Je fais tous ces petits trucs dans une console créée par Delphi.

Pour le double, as-tu tout simplement essayé de remplacer tous les "Real" par des "double" ?
En pricipe, ça devrait aller.

Côté aide, il doit bien y avoir un endroit où sont donnés tous les types basiques du compilateur, normalement avec leur nom, taille, limites, en particulier, pour les types numériques.

Pour savoir la taille d'un type, "real", par exemple, tu peux tout simplement ajouter à ton programme

Code:

writeln('sizeof(real) = ',sizeof(real));
23/01/2007, 18h41
WhiteTigerZ

j'ai essayé meme le type "extended" mais ça marche pas (laisse tomber) lol
Merci quand même ;)
23/01/2007, 23h02
thewho

Bonjour,

Avec la précision des extended, ça devrait marcher, puique celle des double suffit déjà.

As-tu bien corrigé le programme comme je te l'ai donné ?

Il ne faut pas baisser les bras, et essayer de comprendre pourquoi ça ne veut pas marcher chez toi, alors que c'est ok pour moi.
C'est, entre autres, en résolvant de petites choses comme celle-ci qu'on apprend à trouver d'où peuvent venir les problèmes (pas tous, hélas).

salut :) ,
voila le programme à propos peut être il suffit de modifier les options du Turbo pascal 1.5 for Windows (ce que j'utilise) !

Code:

1
2
3
4
5
6
7
8
9
10
11
12
13
14
15
16
17
18
19
20
21
22
23
24
25
26
27
28
29
30
31
32
33
34
35
36
37
38
39
40
 
 
 
program exercice ;
uses wincrt ;
var
   x : extended ;
 
function calcul : extended ;
var
  a,b,c : integer;
  s,s1 : extended ;
begin
  a:=2 ;
  b:=3 ;
  s:= 2;
  c:=1 ;
  repeat
    c:=c+1 ;
    s1:=s ;
    s:=s*(a/b) ;
    if c mod 2 <> 0 then
    begin
      b:=b+2 ;
    end
    else if c mod 2 = 0 then
    begin
      a:=a+2 ;
    end;
  until (abs(s-s1) < 0.000001);
 
  calcul:=s*2 ;
end;
 
begin
  x:=calcul ;
  writeln(x);
 
 
end.

24/01/2007, 10h45
thewho
Bonjour,

Ça a l'air bon, et ça devrait marcher.

Côté options de TPW, je n'ai plus du tout en mémoire ce qui est disponible, et quels sont les réglages possibles, mais je ne pense pas que cela puisse changer la précision des calculs pour un type donné.

Au passage, je n'avais pas fait attention, mais la partie de code mise en rouge
Code:

1 2 3 4 5 6 7 8 if c mod 2 <> 0 then begin b:=b+2 ; end else if c mod 2 = 0 then begin a:=a+2 ; end;
est absolument inutile, étant redondante avec le simple "else"