16 votos

"El más antiguo" error en el sistema de álgebra computacional?

El objetivo de esta pregunta es para encontrar un error en un cálculo mediante un sistema de álgebra computacional, donde la 'correcta' de respuesta (completa con un correcto razonamiento para justificar la respuesta) puede encontrarse en la literatura. Tenga en cuenta que el sistema debe presentar la solicitud para ser capaz de realizar ese cálculo, no de la implementación de un pedazo de (muy viejo) la matemática es triste, pero es un tema diferente.

Por mi conocimiento en el campo, hay un montón de ejemplos del siglo 19 las matemáticas, donde hoy en día el sistema de álgebra computacional obtener la respuesta equivocada. Pero qué tan lejos podemos ir?


Permítanme ilustrar lo que quiero decir. James Bernoulli en las letras a Leibniz (circa 1697-1704) escribió que [en el día de hoy notación, donde voy a suponer que $y$ es una función de $x$ a lo largo] no podía encontrar una forma cerrada a $y' = y^2 + x^2$. En una carta de Noviembre. 15 de 1702, escribió a Leibniz, que fue sin embargo capaz de reducir esto a un 2º orden a la LODE, es decir,$y''/y = -x^2$. Arce puede encontrar (correcto) cerrado-formas para estas dos ecuaciones diferenciales, en términos de funciones de Bessel.

Un ejemplo de lo que es "triste", pero menos interesante es $$r^{n+1}\int_0^{\pi}\cos(r\rho \cos (\omega))\sin(\omega)^{2n+1}d\omega$$ con $n$ supone un entero positivo, $r>0$ $\rho$ real; esto puede ser evaluado como funciones de Bessel, pero, por ejemplo, Arce no. Poisson publicado este resultado en una larga memoria de 1823.

Uno puede quejarse de que (después de Schloemilch, 1857) que él sabía muy bien que $$J_n(z) = \sum_{0}^{\infty} \frac{(-1)^m(z/2)^{n+2*m}}{m!(n+m)!}$$ Arce parece pensar que esta suma está en lugar de $J_n(z)\frac{\Gamma(n+1)}{n!}$, que ningún matemático nunca iba a escribir de esta manera.

Otro ejemplo que pone más cerca de un bug real es que Lommel en 1871 mostró que el Wronskian de $J_{\nu}$$J_{-\nu}$$-2\frac{sin(\nu\pi)}{\nu z}$. Arce puede calcular el Wronskian, pero no se puede simplificar el resultado a $0$. Esto puede ser transformado en un error por el uso de la expresión resultante en un contexto donde nos de la fuerza de la CAS a dividir por él.

Para un verdadero error, considerar $$\int_{0}^{\infty} t^{-\lambda} J_{\mu}(at) J_{\nu}(bt)$$ como investigado por Weber en 1873. Arce devuelve una respuesta incondicional, que a priori parece bien. Si, sin embargo, se formula la misma pregunta pero con $a=b$, no hay respuesta se devuelve! ¿Qué está pasando? Bueno, en realidad esa respuesta sólo es válida para uno de $0\lt a\lt b$ o $0\lt b \lt a$. Pero resulta que (como Watson explica lúcidamente en las páginas 398-404 de su maestro, tratado sobre las funciones de Bessel, esta integral es discontinua para $a=b$. En realidad, la respuesta también es problemático para $\lambda=\mu=0, \nu=1$. Y para los curiosos, la respuesta dada es $$\frac{2^{-\lambda}{a}^{\lambda-1-\nu}{b}^{\nu} \Gamma \left( 1/2\nu+1/2\mu-1/2\lambda+1/2 \right)} { \Gamma\left( 1/2\mu+1/2\lambda+1/2-1/2\nu\right) \Gamma \left( \nu+1 \right)} {F(1/2-1/2\mu-1/2\lambda+1/2\nu,1/2\nu+1/2\mu-1/2\lambda+1/2;\nu+1;{\frac {{b}^{2}}{{a}^{2}}})} $$


EDIT: me preguntó por primera vez a esta pregunta cuando el MO de la comunidad era mucho más pequeño. Ahora que ha crecido mucho, creo que se necesita una segunda vuelta. Una gran cantidad de matemáticos de Casos de uso de manera rutinaria en su trabajo, por lo que no iban a estar interesados en saber la 'edad' de la brecha entre las matemáticas y (confiable) CAS matemáticas?

10voto

Paul Puntos 555

Si recuerdo correctamente de ~30 años, en el Apple ][ el cálculo 7^2 regresaría 49.0001. Hace más de 100 años (o incluso 100 años antes), los matemáticos ya sabía que el cuadrado de un número entero es un número entero.

9voto

Matthew Murdoch Puntos 11530

No estoy seguro de que estoy totalmente de entender la pregunta, pero muchas versiones de Arce dar la respuesta equivocada al contar el número de particiones de $n$, para algunas de las $n$. Obviamente, los matemáticos han sabido cómo hacerlo, ya que en menos de Euler. (Uno podría argumentar que los matemáticos han sabido contar por un tiempo muy largo, de hecho.)

7voto

user1593 Puntos 24

Según Wolfram Alpha y las tablas en [2], $\pi(10^{10}) = 455, 052, 511$. Sin embargo, en Don Zagier de papel enumerados a continuación nos encontramos con que $\pi(10^{10}) = 455, 052, 512$.

Pregunto si alguien ya lo ha señalado esta discrepancia entre las fuentes en otros lugares. Naturalmente, la discrepancia implica la existencia de un error en cualquiera de las rutinas de Zagier o en WA la implementación de la primer función de conteo. No creo que se trata de un error tipográfico en Zagier' nota porque, si mi memoria no me falla, a la derecha, hay algunos otros textos en la literatura que avalan los cálculos de Zagier (por ejemplo, ver [1, página 7].).

Referencias

[1] A. E. Ingham. La distribución de los números primos. Cambridge Biblioteca Matemática, 1934 (Reeditado en 1990).

[2] H. Riesel. Números primos y los Métodos Informáticos para la Factorización. Birkhäuser, Segunda Edición, 1994.

[3] D. Zagier. Los primeros 50 millones de primers. De matemáticas. Intelligencer, 0 (1977).

1voto

Joe Freeman Puntos 133

En mathematica, si usted mira los caracteres de dirichlet módulo 4, que en realidad no llega a los personajes.

0voto

Zurab Silagadze Puntos 5695

En el espacio de Minkowski-tiempo de espera $$\epsilon_{ijkl}\epsilon_{i^\prime j^\prime k^\prime l^\prime}g^{ii^\prime}g^{jj^\prime}g^{kk^\prime}g^{ll^\prime}=-24, \tag{1}$$ donde $\epsilon_{ijkl}$ es la de Levi-Civita tensor y $g^{ij}$ representa la métrica. Sin embargo, si usted (acríticamente) calcular el l.h.s de (1) en la manipulación simbólica del sistema utilizando el FORMULARIO de FixIndex declaración para asignar valores específicos a determinados elementos de la diagonal de la delta de Kronecker, que, por defecto, representa la métrica, se obtiene +24, no -24.

De hecho, este no es un bug, sino una sutileza de funcionamiento interno de la FORMA y de los usuarios fueron advertidos de que uno puede tratar de cambiar el comportamiento de la delta de Kronecker un poco, pero "esto es peligroso y necesidades, además de un buen la comprensión de lo que está sucediendo, buena prueba para asegurarse de que lo que el el usuario quiere es en realidad lo que sucede": http://www.nikhef.nl/~formulario/maindir/documentación/referencia/online/

Curiosamente, esta sutileza de manipulaciones con la de Levi-Civita tensor en FORMA condujo a una larga duró señal de error en los cálculos de la pion-polo dominante plazo en el hadrónica luz-por-dispersión de la luz contribución para el muón momento magnético anómalo, y el origen de este error fue descubierto sólo cuando después de algunos cálculos de tiempo de la misma cantidad, basado en la reducción de Sistema de Álgebra computacional, dio un opuesto signo: http://arxiv.org/abs/hep-ph/0112102 (Comentario sobre el signo de la pseudoscalar polo contribución a la muon g-2, por Masashi Hayakawa y Toichiro Kinoshita).

Creo que el siguiente párrafo en la citada FORMA manual da un muy buen consejo de cómo moderna álgebra computacional los sistemas deben ser utilizados:

"Como en el dicho Zen:

Para el estudiante que comienza las montañas son montañas y el agua es agua. Para el estudiante avanzado de las montañas dejan de ser montañas y el agua deja de ser agua. Para el maestro de las montañas son montañas de nuevo y el agua es el agua de nuevo.

Por supuesto, el maestro moderno también se comprueba que lo que él espera que el sistema, de hecho es lo que el sistema hace."

i-Ciencias.com

I-Ciencias es una comunidad de estudiantes y amantes de la ciencia en la que puedes resolver tus problemas y dudas.
Puedes consultar las preguntas de otros usuarios, hacer tus propias preguntas o resolver las de los demás.

Powered by:

X