[C] Errore calcolo della radice quadrata con sqrt()

[as91]

ciao a tutti, sto cercando di imparare a programmare in c, ho scritto un programma che dati i valori di a,b,c stampa il tipo di equazione e i risultati.

#include <stdio.h>
#include <stdlib.h>

int main()
{
int a;
int b;
int c;
int d2;
double d3;
double x;
double y;
double y2;
//risoluzioni espressioni di primo e secondo grado
printf("inserisci il valore di a \n");
scanf("%d", &a);
printf("inserisci il valore di b \n");
scanf("%d", &b);
printf("inserisci il valore di c \n");
scanf("%d", &c);
			printf("l'espressione inserita e': %d", a);
			printf("x^2 + %d", b);
			printf("x + %d", c);
			printf("=0 \n");

switch (a) {
case 0:
printf("l'equazione e' di 1^ grado \n");
x=(c*-1);
y=x/b;
printf("x=%f\n", y);
break;
default:
printf("l'equazione e' di secondo grado \n");
   switch (b) {
   case 0:
   printf("si tratta di un equazione pura \n");
   x=(c*-1);
39---   y=sqrt(x/a);
   printf("i risultati sono: x1=%f", -y);
   printf("; x2=%f\n", y);
   break;
   default:
   	  	switch (c) {
         case 0:
         printf("si tratta di un' equazione spuria \n");
         x=(b*-1);
         y=(x/a);
         printf("i risultati sono: x1=0; x2=%f \n", y);
         break;
         default:
         d2=(b*b-4*a*c);
53---         d3=sqrt(d2);
         printf("d: %f", d3);
         y=(-b-d3)/2*a;
         printf("%f -b", -b);
         printf("%f y:", y);
         y2=(-b+d3)/2*a;
         printf("i risultati dell'equazione sono: x1=%f", y);
         printf("; x2=%f\n", y2);
         break;
         }
   break;
   }
break;
}

fflush(stdin);
system("pause");
return 0;
}

il programma funzione, ma nelle righe 39 e 53

switch (b) {
      case 0:
      printf("si tratta di un equazione pura \n");
      x=(c*-1);
      y=sqrt(x/a);
      printf("i risultati sono: x1=%f", -y);
      printf("; x2=%f\n", y);
      break;

default:
         d2=(b*b-4*a*c);
         d3=sqrt(d2);
         printf("d: %f", d3);
         y=(-b-d3)/2*a;
         printf("%f -b", -b);
         printf("%f y:", y);
         y2=(-b+d3)/2*a;
         printf("i risultati dell'equazione sono: x1=%f", y);
         printf("; x2=%f\n", y2);
         break;

il programma sbaglia a calcolare la radice. inoltre ho dovuto dividere tutti i calcoli altrimenti i risultati erano completamente sbagliati.

x=(c*-1);
y=x/b;
printf("x=%f\n", y);

!=

x=(c*-1)/b;
printf("x=%f\n", x);

suggerimenti?

[google01]

dal momento che gli operandi (b e c) sono di tipo int, nella versione "corta" dell'operazione il risultato sarà uguale a quello degli operandi, perdendo così di precisione.
il problema non si presenta nella versione estesa perchè, prima della divisione, il risultato di (c*-1) viene convertito automaticamente (nell'assegnazione ad x) in un double e reso quindi pienamente compatibile per la successiva divisione.
le due possibili soluzioni sono:
- cast esplicito

	x=(double) (c*-1)/b;

- cast implicito

	x=(c*-1)/(b*1.0);

[as91]

il primo metodo l'avevo già provato ma il risultato era sbagliato, ma non stò parlando di errori di approsimazioni: sqrt(400)=36500; sqrt(4)=4 || = 0;

sono risultati impossibili...

[linoma]

Se hai questi errori, la FPU avra' generato un errore che non gestisci (vedo solo divisioni), hai la gestione abilitata?
Cmq la vedo strano perche le librerie standard dovrebbero resettare la FPU, se non ricordo male.

[oregon]

Originariamente inviato da as91

x=(c*-1)/b;
printf("x=%f\n", x);

suggerimenti?

Questo suggerimento

x=(c*-1)/(double)b;
printf("x=%f\n", x);

[YuYevon]

Anche includere math.h non sarebbe male comunque.

[as91]

Questo suggerimento

x=(c*-1)/(double)b;
printf("x=%f\n", x);

potrebbe risolvere anche il problema della radice? perchè ho provato a fare un programma semplice

double a=4;
double b=sqrt(a);
printf("%f", b);

Originariamente inviato da YuYevon
Anche includere math.h non sarebbe male comunque.

io uso borland e a volte non lo capisco, pensa che in questo programma faceva funzionere sytem("pause"); senza <stdlib.h>

[oregon]

Cosa ti succede con questo codice?

double a=4;
double b=sqrt(a);
printf("%f", b);

[Zorgatone]

Anchio sto imparando, ho sentito che c'è un header chiamato <math.h>, potrebbe aiuterti?

[oregon]

Come ha già scritto YuYevon, si presuppone che l'include math.h sia usato ...