我需要在一张图上画出8个贝壳。方程是x+y=0
。l**2 * y**2 = (x**2 + y**2)*(y+a)**2
等于l
。因此,我有一个代码:
x, y, a, l = symbols('x, y, a, l')
a = 1
equat = l**2 * y**2 - (x**2 + y**2)*(y+a)**2
equats = [None] * 8
for i in range(1, 9):
equats[i-1] = equat.subs({l : a*i*1/4})
print(equats)
plot_implicit(equats, (x, -5,5), (y, -5, 5))
这段代码引发了一个
{a/4, a/2, 3*a/4..., 2*a}
:---------------------------------------------------------------------------
SympifyError Traceback (most recent call last)
<ipython-input-31-26c783d64bc7> in <module>()
6 equats[i-1] = equat.subs({l : a*i*1/4})
7 print(equats)
----> 8 plot_implicit(equats, (x, -5,5), (y, -5, 5))
~\Anaconda3\lib\site-packages\sympy\plotting\plot_implicit.py in plot_implicit(expr, x_var, y_var, **kwargs)
316
317 elif not isinstance(expr, Relational):
--> 318 expr = Eq(expr, 0)
319 has_equality = True
320 elif isinstance(expr, (Equality, GreaterThan, LessThan)):
~\Anaconda3\lib\site-packages\sympy\core\relational.py in __new__(cls, lhs, rhs, **options)
290 from sympy.simplify.simplify import clear_coefficients
291
--> 292 lhs = _sympify(lhs)
293 rhs = _sympify(rhs)
294
~\Anaconda3\lib\site-packages\sympy\core\sympify.py in _sympify(a)
385
386 """
--> 387 return sympify(a, strict=True)
388
389
~\Anaconda3\lib\site-packages\sympy\core\sympify.py in sympify(a, locals, convert_xor, strict, rational, evaluate)
301
302 if strict:
--> 303 raise SympifyError(a)
304
305 try:
SympifyError: SympifyError: [0.0625*y**2 - (x**2 + y**2)*(y + 1)**2, 0.25*y**2 - (x**2 + y**2)*(y + 1)**2, 0.5625*y**2 - (x**2 + y**2)*(y + 1)**2, 1.0*y**2 - (x**2 + y**2)*(y + 1)**2, 1.5625*y**2 - (x**2 + y**2)*(y + 1)**2, 2.25*y**2 - (x**2 + y**2)*(y + 1)**2, 3.0625*y**2 - (x**2 + y**2)*(y + 1)**2, 4.0*y**2 - (x**2 + y**2)*(y + 1)**2]
我怎样才能解决这个问题?我在网上搜索,但没有找到任何解决办法。
最佳答案
如果一个一个地绘制方程而不是传递一个list
的方程,问题可能就解决了。
p = plot_implicit(equats[0], (x, -5,5), (y, -5,5), show=False)
for eq in equats[1:]:
p2 = plot_implicit(eq, (x, -5,5), (y, -5,5), show=False)
p.append(p2[0])
p.show()
关于python - 获取SimpifyError,我们在Stack Overflow上找到一个类似的问题:https://stackoverflow.com/questions/52974515/