Question

我尝试阅读其他提到此问题的帖子，但我不明白如何将其应用于我的情况。这也是我第一次使用堆栈溢出，因此对于任何格式错误，我深表歉意。

from sympy import *
from math import sqrt
x = symbols('x')
a = symbols('a',positive=True)
f = 1 / ((x**4)*((x**2)-(a**2))**.5)
#Part A
g = f.subs(4*a,x)
h = f.subs(2*a,x)
area=g-h
#Part B
answ = solve(area-5,a)
print ('part a:' , str(area))
print ('part b:' , str(answ[0].evalf()))

我正在尝试获取要打印的值。

此代码块中也有相同的问题

from sympy import *
x = symbols('x')
# Part A
f = (x**3-4*x+3) / (((x-5)**2)*(x**2+3)*(x**2+5))
print("The partial fraction decomposition of f is :", apart(f), "\n")

# Part B
a = symbols('a', positive=True, constant=True)
g = (x**2+2*x+a) / (x-1)**3*(x**2+1)
print("The partial fraction decomposition of g is :", apart(g, x), "\n")

# Part C
g = apart(g, x)
G1 = integrate(x, (x, 2, 4))
G2 = integrate(2*(a + 3)/(x - 1)**3, (x, 2, 4))
G3 = integrate(2 * (a + 7) / (x - 1)**2, (x, 2, 4))
G4 = integrate((a + 13)/(x - 1), (x, 2, 4))
G5 = integrate(5, (x, 2, 4))
print(G1, G2, G3, G4, G5)
value = solve((G1 + G2 + G3 + G4 + G5 - 1), a)

print((value))

由于几个人似乎感到困惑，documentation是我要使用python解决的问题。

Answer 1

因此，在与我的助教交谈了一段时间后，我想出了答案，如果有人好奇，我将在此处发布解决方案和解释。

joinResolvers

问题的一部分是为sqrt（）函数导入数学运算，因为它与sympy的运算符冲突，从而使得无法解决任何问题。

from sympy import *
from sympy.plotting import (plot, plot_parametric,plot3d_parametric_surface, plot3d_parametric_line, plot3d)

下一组问题是我误用了subs（）函数，将要解决的变量放在了错误的一边。我还预先集成了功能，以便实际解决曲线下的面积。为了获得answ的值，我使用了与evalf（）一起创建的向量的第一个值来接近它。

# Insert code for problem #1 here.
x = symbols('x')
a = symbols('a', positive=True)
f = 1 / ((x**4)*(sqrt((x**2)-(a**2))))
#Part A
F = integrate(f, x)
g = F.subs(x, 4*a)
h = F.subs(x, 2*a)
area = g - h
#Part B
# print(area.evalf())
answ = solve(area-5, a)
print ('part a:', area)
print ('part b:', answ[0].evalf())

下一组代码的问题是g（x）函数上的一些简单语法错误。我还必须删除constant = True，这是因为它使解决不可能。

当我在sympy中使用solve（）时，变量的空列表？

1 个答案: