sympy.Eq(func, var)ラムダで使用して生成した方程式を使用しようとしています。「追加」オブジェクトのリストを返しているようですが、使用方法がわかりません。strを使用する前に型キャストを試みましたevalがfrom sympy import I、問題と関係があると思います。
コードとエラーメッセージは次のとおりです。
py solid.py -vy = (x-1.0) (x-3.0) *2 1.0 3.0 -軸 y --メソッド ディスク
実行例を次に示します。
import sys
import argparse
import math
from scipy import integrate
from sympy import Symbol, Eq, solve, I
import numpy
import pylab
parser = argparse.ArgumentParser(description='Find the volume of the solid of rotation defined. Use * for multiplication and ** to raise a power. Bounds are set equal to the implicit variable. Trailing .0\'s must be used for constants, and parentheses must be escaped ( e.g. \(x-2\) ). Answers are accurate to 6 decimals')
parser.add_argument('var', choices='xy', help='x|y')
parser.add_argument('=', metavar='=', choices='=')
parser.add_argument('equation', help='continous function')
parser.add_argument('a', type=float, help='bound \'a\'')
parser.add_argument('b', type=float, help='bound \'b\'')
parser.add_argument('-axis', choices='xy', help='axis of revolution (default == x)')
parser.add_argument('-m', choices='ds', help='method (disk/shell)')
parser.add_argument('-g', action='store_true', help='show graph')
parser.add_argument('-v', action='store_true', help='verbose (for debugging)')
args = parser.parse_args()
y = Symbol('y')
x = Symbol('x')
def solve_for_implicit(func, var):
if var == 'x':
equation = Eq(eval(func), x)
func = solve(equation, y)
else:
equation = Eq(eval(func), y)
func = solve(equation, x)
return func
def volume(var, func, a, b, axis=None, method=None):
if axis == None: axis = 'x'
if method == 's': method = 'shell'
if method == 'd': method = 'disk'
if var == axis and axis == 'x':
if args.v: print 'x = y about x'
if not method == 'disk':
pass
else:
pass
elif var == axis and axis == 'y':
# SHELL METHOD
if args.v: print 'y = x about y'
if not method == 'disk':
# this is what should be used for y = x about y
if args.v: print 'using SHELL method'
func = eval('lambda x: ' + 'x*(' + func + ')')
integral = integrate.quad(func, a, b)
if args.v: print 'integral from', a, 'to', b, '=', integral[0]
answer = 2 * math.pi * integral[0]
if args.v: print '2*pi*', integral[0], '=', answer
else:
func = solve_for_implicit(func, var)
print func[0]
func = eval('lambda y: ' + '(' + str(func[0]) + ')**2')
integral = integrate.quad(func, a, b)
if args.v: print 'integral=', integral[0]
answer = math.pi * integral[0]
elif not var == axis and axis == 'y':
# DISK METHOD
if args.v: print 'x = y about y -- '
if not method == 'shell':
pass
else:
pass
elif not var == axis and axis == 'x':
# DISK
if args.v: print 'y = x about x --',
if not method == 'shell':
pass
else:
pass
return answer
print volume(args.var, args.equation, args.a, args.b, args.axis, args.m)
if args.g: graph(args.equation, args.a, args.b)
トレースバック (最新の呼び出しが最後):
ファイル "solid.py"、136 行目、印刷ボリューム (args.var、args.equation、args.a、args.b、args.axis、args.m) 内 ファイル "solid. py"、75 行目、ボリューム func = eval('lambda y: ' + '(' + func[0] + ')**2') TypeError: 'str' オブジェクトと 'Add' オブジェクトを連結できません
そして、最初にそれを str() に型キャストしようとすると、次のようになります:
ファイル「solid.py」、136 行目、印刷ボリューム (args.var、args.equation、args.a、args.b、args.axis、args.m) ファイル「solid.py」、76 行目、ボリュームintegral = integrate.quad(func, a, b) ファイル "/usr/lib/python2.6/dist-packages/scipy/integrate/quadpack.py"、185 行目、quad retval = _quad(func,a,b) ,args,full_output,epsabs,epsrel,limit,points) ファイル "/usr/lib/python2.6/dist-packages/scipy/integrate/quadpack.py", 行 249, in _quad return _quadpack._qagse(func,a ,b,args,full_output,epsabs,epsrel,limit) quadpack.error: 指定された関数は有効な浮動小数点数を返しません。
また、print func[0]与えます：
2.3333333333333 + (0.296296296296298 - y/2 + (-0.0877914951989024 + (0.592592592592595 - y) 2/4) (1/2)) (1/3)*(1/2 - I*3 (1/2)/2) + 0.444444444444444/((1/2 - I*3* (1/2)/2) (0.296296296296298 - y/2 + (-0.0877914951989024 + (0.592592592592595 - y) 2/4) (1/2))**( 1/3))
以前に完全なスクリプトを提供できなかったことをお詫び申し上げます。