El problema está en que en la parte del cálculo en la que utilizas (1/2)
.
28: r=(1/2)*(ra+(num/ra));
Al realizar esta operación aritmética, el compilador interpreta la operación como una división entre enteros, dando como resultado 0.5. Por lo tanto, al arrojar 0
como resultado de esa operación y multiplicarla por el resto de la expresión, el valor de r
se convierte a 0
.
En esta pregunta de Stack Overflow en inglés podrás ver una mejor explicación del tema:
Why is (1/2)*x different from 0.5*x?
Para que el compilador lo tome como un número flotante, o con parte decimal, tienes que cambiar los valores por:
(1/2)
por (1.0/2.0)
ó
(1/2)
por 0.5
Este tema se explica en el siguiente capítulo del curso:
9. Conversión de tipos.
En estos casos, cuando los operandos de cada operación binaria
asociados a un operador son de distinto tipo, el compilador los
convierte a un tipo común. Existen reglas que rigen estas
conversiones, y aunque pueden cambiar ligeramente de un compilador a
otro, en general serán más o menos así:
- Cualquier tipo entero pequeño como char o short es convertido a int o unsigned int. En este punto cualquier pareja de operandos será int (con o sin signo), long, long long, double, float o long double.
- Si un operando es de tipo long double, el otro se convertirá a long double.
- Si un operando es de tipo double, el otro se convertirá a double.
- Si un operando es de tipo float, el otro se convertirá a float.
- Si un operando es de tipo unsigned long long, el otro se convertirá a unsigned long long.
- Si un operando es de tipo long long, el otro se convertirá a long long.
- Si un operando es de tipo unsigned long, el otro se convertirá a unsigned long.
- Si un operando es de tipo long, el otro se convertirá a long.
- Si un operando es de tipo unsigned int, el otro se convertirá a unsigned int.
- Llegados a este punto ambos operandos son int.
Este caso puede explicarse con un ejemplo de la cita del curso:
4. Operadores I.
Veamos, por ejemplo, la expresión 17 / 7, es decir 17 dividido entre
7. Cuando aprendimos a dividir, antes de que supiéramos sacar decimales, nos enseñaron que el resultado de esta operación era 2, y
el resto 3, es decir 2*7+3 = 17.
En C++, cuando las expresiones que intervienen en una de estas
operaciones sean enteras, el resultado también será entero, es decir,
si 17 y 7 se almacenan en variables enteras, el resultado será entero,
en este caso 2.
Por otro lado si las expresiones son en punto flotante, con decimales,
el resultado será en punto flotante, es decir, 2.428571. En este caso:
7*2.428571=16.999997, o sea, que no hay resto, o es muy pequeño.
Por eso mismo, calcular el resto, usando el operador %, sólo tiene
sentido si las expresiones implicadas son enteras, ya que en caso
contrario se calcularán tantos decimales como permita la precisión de
tipo utilizado.
Siguiendo con el ejemplo, la expresión 17 % 7 dará como resultado 3,
que es el resto de la división entera de 17 dividido entre 7.
Quedando así tu código corregido:
#include<stdio.h>
#include<math.h>
float calculo_raiz(float,float);
int main()
{
float rta;
float numero;
float tol;
printf("\n\t\t ¿¿¿ AVERIGUA LA RAIZ CUADRADA ???\n\n" );
printf("Usted desea averiguar la raiz del numero:");
scanf("%f",&numero);
printf("\nIngrese la tolerancia que admitira en el resultado:" );
scanf("%f",&tol);
rta=calculo_raiz(numero,tol);
printf("\nLa raiz de %.2f es %.2f\n",numero,rta);
return 0;
}
float calculo_raiz(float num,float tolerancia)
{
float r=0;
float ra=1;
int inc=1;
while(inc>=1)
{
r=(1.0/2.0)*(ra+(num/ra)); // ó r=0.5*(ra+(num/ra));
if((r-ra)<tolerancia)
{
return r;
}
else
{
inc++;
ra=r;
}
}
}