Estoy intentando hacer una función de seno, mediante las series de Taylor. Este algoritmo funciona en Java y python sin embargo en c++ me produce errores ya que al calcular con n términos grandes el numero se convierte en nan o inf.
Este es el algoritmo en C++
#include<iostream>
#include<math.h>
#include<iomanip>
#include<string>
#include<cmath>
#include<limits.h>
using namespace std;
#define Math_PI 3.14159265358979323846
class Math{
public:
int factorial(int n){
int factor = 1;
for(int i=n;i>0;i--){
factor*=i;
}
return factor;
}
float radianes(float d){
return (d*Math_PI)/180;
}
double sen(double x){
double seno = 0;
int sentido = 1;
for(int i=1;i<100;i+=2){
seno+=sentido*(pow(x,i)/factorial(i));
sentido*=-1;
}
return round(seno);
}
};
int main(){
Math calcu;
float degrees;
cout<<"Digita el angulo: ";cin>>degrees;
degrees = calcu.radianes(degrees);
cout<<"El angulo en radianes es: "<<degrees<<endl;
cout<<"El seno del angulo es: "<<roundf(sin(degrees))<<endl;
cout<<"El seno del angulo es: ";
cout<<calcu.sen(degrees);
return 0;
}
Este es el mismo algoritmo y procedimiento en python, sin embargo en este no genera error alguno en el calculo
import math
def factorial(n):
factor = 1
for i in range(1,n+1):
factor*=i
return factor
def sen(degree):
seno = 0
sentido = 1
i = 1
while(i<100):
seno+=sentido*(pow(degree ,i)/factorial(i))
i+=2
sentido*=-1
seno = round(seno)
return seno
def cos(degree):
coseno = 0
sentido = 1
i = 0
while(i<100):
coseno+=sentido*(pow(degree, i)/factorial(i))
i+=2
sentido*=-1
coseno = round(coseno)
return coseno
def euler(exp):
e = 1
i = 1
for i in range(101):
e+=pow(5, i)/factorial(i)
e = round(e)
return e
def main():
numero = float(input("Digita un angulo: "))
degree = (numero*math.pi)/180
print(f'El seno de {degree} es : {sen(degree)}')
print(f'El coseno de {degree} es: {cos(degree)}')
exponente_euler = int(input("Digita el exponente de e: "))
print(f'El resultado de elevar e^{exponente_euler} es: {euler(exponente_euler)}')
if __name__ == "__main__":
main()
Al ejecutar el codigo en C++ obtengo esta salida:
Digita el angulo: 90
El angulo en radianes es: 1.5708
El seno del angulo es: 1
El seno del angulo es: nan
Y quisiera saber como puedo reducir la longitud del double seno, ya que en cierto punto del buble for se elevan demasiado los decimales y se transforma en inf, y posteriormente en un nan.
round
redondea al número entero más cercano y la claseMath
es innecesaria ya que no guarda ningún estado, solo la usas para contener las funciones. Si lo que querías es "agrupar" de alguna manera las funciones, puedes ponerlas en un namespace.