Me acordé que python tiene en sympy una libaría para generar series de Taylor . Esto sirve para cualquier función que se escriba en la variable $$f$$ en este caso el arcoseno hiperbólico
import sympy as sp
a,b = map(float, input("Ingresa el punto donde expandir la serie y el punto donde evaluar separado por comas para P").split(","))
x = sp.symbols("x")
f = sp.asinh(x)
# Taylor de asinh(x) alrededor de x=a, hasta grado 5 cambiarlo para mas grados .. remove elimina el error
polinomio = f.series(x, a, 6).removeO()
print()
print(polinomio)
print()
# Evaluarlo en x=b
valor = polinomio.subs(x, b)
print(float(valor))