Programme racine carrée

Version imprimable

Bonsoir,

Sauriez-vous me dire pourquoi ce programme ne fonctionne pas?

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
41
42
43
44
#include<stdio.h>
#include<stdlib.h>
 
long double racine(long double x);
 
main()
{
   long double y,m;
   printf("x=");
   scanf("%lf",&y);
   printf("\aracine carree = %lf\n",racine(y));
   system("PAUSE");
   return 0;
}
 
long double racine(long double x)
{
   long double n,a,b,eps;
   int Trouve=0;
   printf("eps=");
   scanf("%lf",&eps);
   if (x>=1)
   {
      a=1;
      b=x;
   }
   else
   {
      a=x;
      b=1;
   }
   do
   {
      n=(a+b)/2;
      if (n*n==x)
         Trouve=1;
      else if(n*n>x)
         b=n;
      else a=n;
   }
   while(Trouve==0 && (b-a)>eps );
 
   return n;
}

15/01/2009, 20h56
tonton fred

Bonjour,

A vue de nez ca va bloquer la

Code:

if (n*n==x)

n et x sont des double. Il ne faut jamais tester l'egalite avec des nombres flottants a cause de l'imprecision de la machine. Ceci devrait marcher (incertitude a definir selon ce que tu veux):

Code:

if (n*n - x =< incertitude)

Salut,

Citation:

Envoyé par yo_haha

Bonsoir,

Sauriez-vous me dire pourquoi ce programme ne fonctionne pas?

Il faudrait peut-être résoudre aussi ces warnings:

Code:

1
2
3
4
5
6
7
8
9
  === dummy, Debug === 
C:\data\sources\C\dummy\main.c 7 warning: return type defaults to `int' 
C:\data\sources\C\dummy\main.c  In function `main': 
C:\data\sources\C\dummy\main.c 10 warning: double format, long double arg (arg 2) 
C:\data\sources\C\dummy\main.c 11 warning: double format, long double arg (arg 2) 
C:\data\sources\C\dummy\main.c 8 warning: unused variable `m' 
C:\data\sources\C\dummy\main.c  In function `racine': 
C:\data\sources\C\dummy\main.c 20 warning: double format, long double arg (arg 2) 
  === Build finished: 0 errors, 5 warnings ===

A+

Pfeuh

15/01/2009, 21h54
Invité(e)
Bonsoir,

Pour entrer un long double, il faut utliser "%Lf" ou "%llf" (voir pages man).

Il serait aussi intéressant de vérifier que l'utilisateur n'a pas rentré de valeurs aberrantes (0 ou valeurs négatives).

Sinon, le programme a l'air de fonctionner (une fois les warning corrigés) :
Code:

1 2 3 4 5 6 7 >./a.out x=15 (x)=15.000000 eps=.1 (eps)=0.100000 racine carree = 3.898438 sh: line 1: PAUSE: command not found
15/01/2009, 22h56
yo_haha

Merci pour vos réponses.
16/01/2009, 07h18
Invité(e)

Citation:

Envoyé par yo_haha

Merci pour vos réponses.

Peux tu nous poster le code final ?

Merci.

Le voilà!!!

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
#include<stdio.h>
#include<stdlib.h>
 
double racine(float x);
 
main()
{
float y;
printf("x=");
scanf("%f",&y);
printf("\aracine carree = %f\n",racine(y));
system("PAUSE");
return 0;
}
double racine(float x)
{
float a, b, n, eps;
printf("eps=");
scanf("%f",&eps);
if (x>1)
{
a=1;
b=x;
}
else
{
a=x;
b=1;
}
n=(float)(a+b)/2;
while (n*n!=x && (b-a)>eps )
{
if (n*n<x)
a=n;
else
b=n;
n=(float)(a+b)/2;
}
return n;
}

19/01/2009, 17h35
Porteuris

Pardon j'ai testé ce programme et je ne vois pas à quoi correspond eps
merci.
19/01/2009, 18h16
ram-0000

Citation:

Envoyé par Porteuris

Pardon j'ai testé ce programme et je ne vois pas à quoi correspond eps
merci.

epsilon (une toute petite valeur)
19/01/2009, 20h00
yo_haha

eps est la précision de la solution (la racine du nombre entré). Si la racine est réelle mais n'est pas entière, je dois préciser à la machine à quel point elle doit arrêter la recherche de la racine sinon elle recherchera infiniment la racine. Ainsi, elle me renvoie la racine du nombre que j'ai entré une fois que la longeur b-a de l'intervalle [a, b] qui contient cette racine est inférieure à eps. La racine renvoyé sera donc, le milieu de cet intervalle (ce n'est pas la racine exacte mais c'est une valeur approchée à eps près).
La racine est plus exacte si eps est très petite.
Pour plus d'information, se renseigner sur la recherche dichotomique d'une racine carrée.
J'espère avoir été clair et compréhensible.