Python (sympy) TypeError: не может конкатенировать объекты 'str' и 'Add'

Question

Я пытаюсь использовать уравнение, которое я создал с помощью sympy.Eq(func, var) в лямбда. Кажется, он возвращает список объектов «Добавить», и я не уверен, как их использовать. Я попробовал typecasting до str перед использованием eval, и я попробовал from sympy import I, что, я думаю, имеет какое-то отношение к этой проблеме.

Вот код, и сообщение об ошибке:

py solid.py -v y = (x-1.0)(x-3.0)*2 1.0 3.0 -axis y --method disk

---

Вот бегущий пример:

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) 

Traceback (most recent call last):
File "solid.py", line 136, in print volume(args.var, args.equation, args.a, args.b, args.axis, args.m) File "solid.py", line 75, in volume func = eval('lambda y: ' + '(' + func[0] + ')**2') TypeError: cannot concatenate 'str' and 'Add' objects

и когда я пытаюсь типажей его на ул() сначала я получаю это:

File "solid.py", line 136, in print volume(args.var, args.equation, args.a, args.b, args.axis, args.m) File "solid.py", line 76, in volume integral = integrate.quad(func, a, b) File "/usr/lib/python2.6/dist-packages/scipy/integrate/quadpack.py", line 185, in quad retval = _quad(func,a,b,args,full_output,epsabs,epsrel,limit,points) File "/usr/lib/python2.6/dist-packages/scipy/integrate/quadpack.py", line 249, in _quad return _quadpack._qagse(func,a,b,args,full_output,epsabs,epsrel,limit) quadpack.error: Supplied function does not return a valid float.

также print func[0] дает:

2.33333333333333 + (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))

спасибо за помощь всем, извините за не предоставление полного сценария раньше.

Answer 1

«Я пытался типажей к ASTR перед использованием Eval, и я пытался от SymPy импорта I, который я думаю, что-то делать с этой проблемой»

Да, без сомнения, что именно эта проблема. Я уверен, что у Numpy есть документация о том, как использовать объекты Добавить. Преобразование их в строки, а затем их оценка, безусловно, не является правильным решением.

Кроме ознакомления с документацией, вы также можете осмотреть объекты с помощью dir (func), чтобы узнать, какие функции и атрибуты у него есть. В интерпретаторе есть также отличная функция справки, поэтому вы можете ввести help (func) и получить информацию.