Discussion :

[goldensun]

J'essaie de calculer la racine carré en utilisant la méthode de héron, mais ça marche pas...la condition du if c'est que dans la suite des racines carrés xi+1-xi/xi doit être inférieur à une précision p que l'on aura fixé. Et la je comprends pas ce qui va pas. je veux une précision de 0.01 chiffre après la virgule je mets donc fabs((y-i)/i)<p mais ça donne rien, le résultat est absurde .

Code :

Sélectionner tout - Visualiser dans une fenêtre à part

1
2
3
4
5
6
7
8
9
10
11
12
13
14
15
 
#include <stdio.h>
#include <stdlib.h>
#include <math.h>
#define n 25.0
int main()
{float y,i,p;
 
p=0.001;
for(i=1;fabs((y-i)/i)<p;i=y)
y=(i+(n/i))/2;
printf("la racine carré de %.2f vaut %.2f",n,y);
 
return 0;
}

[Emmanuel Delahaye]

J'essaie de calculer la racine carré en utilisant la méthode de héron, mais ça marche pas...

Est-ce suffisamment clair ?

Code :

Sélectionner tout - Visualiser dans une fenêtre à part

1
2
3
4
5
6
7
8
9
10
11
12
 
 
-------------- Build: Debug in hello ---------------
 
Compiling: main.c
Linking console executable: bin\Debug\hello.exe
C:\dev\hello\main.c:7: warning: function declaration isn't a prototype
C:\dev\hello\main.c: In function `main':
C:\dev\hello\main.c:7: warning: 'y' might be used uninitialized in this function
Output size is 18.55 KB
Process terminated with status 0 (0 minutes, 1 seconds)
0 errors, 2 warnings

[goldensun]

Nan pas trop, le soucis c'est que j'ai pas initialisé y?
Pour le warning sur le prototype j'ai pas comprit, j'utilise qu'une fonction c'est fabs, et elle est dans la bibliothèque math.

[Emmanuel Delahaye]

Nan pas trop, le soucis c'est que j'ai pas initialisé y?

Oui.

Pour le warning sur le prototype j'ai pas comprit,

Ici, main() n'a pas de paramètres, il faut le préciser avec void

Code :

Sélectionner tout - Visualiser dans une fenêtre à part

int main (void)

[goldensun]

Je peux initialiser y à n'importe quelle valeur, puisque de toute façon au premier passage de la boucle for, la valeur de y va être réaffecter?
Mais ça marche toujours pas.

[diogene]

Code :

Sélectionner tout - Visualiser dans une fenêtre à part

1
2
for(i=1;fabs((y-i)/i)<p;i=y)
...

A l'entrée dans la boucle, l'expression fabs((y-i)/i)<p est évaluée, et à ce moment, y n'est pas initialisé et sa valeur est quelconque (et peut être que de ce fait, tu n'entres jamais dans la boucle et y n'est alors jamais initialisé)

[goldensun]

Si j'initialise y avec une valeur quelquonque ça doit donc marcher, par exemple si je prends un y qui vaut 4 on aura alors bien fabs((y-i)/i))<p peut importe la valeur avec laquelle j'initialise y du moment que l'on vérifie ça, puisque de toute manière y va être réaffecter dans la boucle.

for(i=1;fabs((y-i)/i)<p;i=y)
...

Le soucis c'est que même si je fais ça, ça me renvoie une valeur absurde.

for(i=1,y=4;fabs((y-i)/i)<p;i=y)
...

[3DArchi]

Code :

Sélectionner tout - Visualiser dans une fenêtre à part

1
2
3
 
p=0.001;
for(i=1,y=4;fabs((y-i)/i)<p;i=y)

Quand je calcule à la main fabs(4-1)/1=3>0.001
Donc ça ne rentre pas dans la boucle.
On n'initialise pas y (ou quoique ce soit d'autre) avec "n'importe quelle valeur". Il faut que ton initialisation, ton test et ta boucle est un sens les uns avec les autres.
Peut être doit-tu utiliser un do/while si ton test de sortie n'est pas strictement lié à ton itérateur de boucle.

[goldensun]

Merci en effet comme ça ça marche (je m'étais déjà trompé dans la condition de la boucle for)

Code :

Sélectionner tout - Visualiser dans une fenêtre à part

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
 
#include <stdio.h>
#include <stdlib.h>
#include <math.h>
#define n 25.00
 
int main(void)
{float  y,i,p;
p=0.01;
i=1.00;
y=4.00;
 
 
do
{
    i=y;
    y=(i+(n/i))/2;
}
while(fabs((y-i)/i)>p);
 
 
printf("la racine carré de %.2f vaut %.2f",n,y);
 
    return 0;
}

Par contre la j'initialise y avec une valeur quelconque. Et je suis obligé car y et i désignent les éléments d'une même suite, il y a juste un décalage d'un élément entre eux deux.

[Emmanuel Delahaye]

Merci en effet comme ça ça marche (je m'étais déjà trompé dans la condition de la boucle for)

Ton algo est bizarre :

Code :

Sélectionner tout - Visualiser dans une fenêtre à part

1
2
3
4
5
6
7
 
  i = 1.00;
   y = 4.00;
 
   do
   {
      i = y;

pourquoi affecter 2 fois de suite (de plus, avec des valeurs différentes) ?

[goldensun]

Parce qu'il faut faire l'affectation en dehors de la boucle nan?

[goldensun]

Donc je devrai faire ça

Code :

Sélectionner tout - Visualiser dans une fenêtre à part

1
2
3
4
5
6
7
8
9
10
11
12
13
14
15
16
17
18
19
20
21
22
23
24
 
#include <stdio.h>
#include <stdlib.h>
#include <math.h>
#define n 25.00
 
int main(void)
{float  y,i,p;
p=0.01;
 
y=1.00;
 
do
{
    i=y;
    y=(i+(n/i))/2;
}
while(fabs((y-i)/i)>p);
 
 
printf("la racine carré de %.2f vaut %.2f",n,i);
 
    return 0;
}

nan j'ai rien dire ça marche impec

[3DArchi]

OK,
Repartons de la méthode de héron:
Pour le calcul de la racine carré de A, on utilise la suite:
x(n+1) = (x(n) + A/x(n))/2.
On demande que x(0) soit assez proche de racine(A).
Évidement, tu veux stopper ton calcul quand x(n+1)-x(n) devient négligeable.
Donc ton algo serait quelque chose comme ça:

Initialiser A (c'est, si je comprend ton code ton #define n 25.00)
Initialiser y avec une valeur proche de racine(A) (et pas au hasard!). y sera notre x(n).
Initialiser p avec ton seuil de précision
Initialiser x = 0 (car c'est une bonne pratique d'initialiser toute variable déclarée). x sera notre x(n+1)
Initialiser delta=0 (idem précédent)
Faire
x = (y + A/y)/2
delta = abs(y-x)
y = x
Tant que (delta>p)

La vrai question est comment initialiser y avec une valeur proche de racine(A) puisque c'est ce qu'on cherche! Probablement en cherchant un algo plus rapide mais plus imprécis, ou en utilisant une seule récursion d'un autre algo. Bref, quelques recherches s'imposent.

[goldensun]

Nan x(o) peut être n'importe quelle valeur du moment qu'elle est différent de 0. Il y a pas besoin qu'elle soit proche de racine de A, c'est l'algo qui nous rapproche de racine de A à partir d'une valeur quelconque de x(0).

[3DArchi]

Nan x(o) peut être n'importe quelle valeur du moment qu'elle est différent de 0.

Effectivement. Ne connaissant pas la méthode, j'ai été sur wiki, mais il y a la boulette. Enfin, disons plutôt que plus x(0) sera proche de la solution, plus la convergence sera rapide.