Estoy tratando de evaluar al jacobiano en (x, y) = (0,0) pero no puedo hacerlo. import sympy as sp from sympy import * import numpy as np x,y=sp.symbols(‘x,y’, real=True) J = Function(‘J’)(x,y) f1=-y f2=x – 3*y*(1-x**2) f1x=diff(f1,x) f1y=diff(f1,y) f2x=diff(f2,x) f2y=diff(f2,y) J=np.array([[f1x,f1y],[f2x,f2y]]) J1=J(0,0) print J1 El error correspondiente a —> 16 J1=J(0,0) es TypeError: ‘numpy.ndarray’ object […]
Cuando se utiliza solve para calcular las raíces de una ecuación cuadrática, SymPy devuelve expresiones que podrían simplificarse pero no puedo lograr que las simplifique. Un ejemplo mínimo se ve así: from sympy import * sqrt(-24-70*I) Aquí, SymPy simplemente devuelve sqrt(-24-70*I) mientras que Mathematica o Maple responderán con el equivalente de 5-7*I Soy consciente de […]
Soy bastante nuevo en SymPy y tengo lo que podría ser una pregunta básica. O simplemente podría estar malinterpretando cómo se debe usar SymPy. ¿Hay una manera de crear una expresión que no esté representada por átomos, sino por una combinación de otras expresiones? Ejemplo: >>> from sympy.physics.units import * >>> expr1 = m/s >>> […]
Me gustaría lambdificar una expresión simbólica que contiene la función erf con SymPy. Esto se puede hacer para los argumentos escalares de la siguiente manera: log_normal = 0.5 + 0.5 * sym.erf((sym.log(x) – mu) / sym.sqrt(2 * sigma**2)) F = sym.lambdify([x, mu, sigma], log_normal) F(1.0, 0.0, 1.0) Me gustaría vectorizar lo anterior. Normalmente haría lo […]
Para mi proyecto, necesito resolver para una matriz X dadas matrices Y y K. (XY = K) Los elementos de cada matriz deben ser números enteros módulo una prima aleatoria de 256 bits. Mi primer bash de resolver este problema usó la función mod_inv(n) . El problema con esto es que me estoy quedando sin […]
Anteriormente he usado la impresión bonita de las matemáticas en el cuaderno de ipython. Después de actualizar a jupyter (también actualiza muchos otros paquetes relacionados con ipython), la impresión bastante ya no funciona como antes. Utilizo este código en la parte superior de mis cuadernos para configurarlo: import sympy as sp sp.init_printing() También he intentado […]
Tengo un progtwig en el que un usuario ingresa una función, como sin(x)+1 . Estoy usando ast para tratar de determinar si la cadena es ‘segura’ al incluir en la lista blanca los componentes como se muestra en esta respuesta . Ahora me gustaría analizar la cadena para agregar signos de multiplicación ( * ) […]
R (teta, phi) = cos (phi ^ 2), teta [0, 2 * pi], phi [0, pi] ¿Cómo dibujar una gráfica de esta función (R (teta, phi)) en coordenadas esféricas con la ayuda de matplotlib? La documentación no he encontrado coordenadas esféricas.
Quiero analizar las fórmulas de LaTeX y usarlas directamente como expresiones de SymPy. En otras palabras, lo que necesito es algo similar a sympify: from sympy import sympify f = sympify(‘x^2 + sin(y) + 1/2’) print f pero eso puede tomar expresiones LaTeX (cadenas) como entrada, por ejemplo: f = latex_sympify(‘\frac{x}{1+x}’) Dado que sympify es […]
Estoy tratando de encontrar la solución óptima para el siguiente sistema de ecuaciones en Python: (x-x1)^2 + (y-y1)^2 – r1^2 = 0 (x-x2)^2 + (y-y2)^2 – r2^2 = 0 (x-x3)^2 + (y-y3)^2 – r3^2 = 0 Dados los valores, un punto (x, y) y un radio (r): x1, y1, r1 = (0, 0, 0.88) x2, […]