修正'ImmutableDenseNDimArray‘对象没有属性'could_extract_minus_sign’

Question

from sympy import *
import numpy as np
import matplotlib.pyplot as plt

x = symbols('x')
y = symbols('y')
C = symbols('C')

# y' + 2/x (y) = sin(x) 

Intfact = exp((integrate(2/x,x)))
print(Intfact)
Zsoln = ((integrate((sin(x)/x)*Intfact,x)+C)/Intfact)
print('y(x) =', Zsoln.expand()) 

x = np.array(range(0,100))
y = (-2/x**2) - (cos(x)/x) + (sin(x)/x**2)
plt.plot(x,y)
plt.show()

错误出现在'y = (-2/x2) - (cos(x)/x) + (sin(x)/x2)‘，我需要在C等于-2，-1,0,1,2的情况下画出这个整型。

Answer 1

嗨，你可以这样写：

y = (np.divide(-1, np.power(x, 2))) - (np.divide(np.cos(x), x)) + (np.divide(np.sin(x), np.power(x, 2)))