errors in the integrate function

[AmplitudeGravity]

[mzaffalon]

It looks like you are missing a * before (1+x).

In the future, could you please post code, not screenshots, to make it easier for others to reproduce your report?

[AmplitudeGravity]

Thanks. That seems more errors. integrate(x^(1/2)*(1+x)^(-23/10),(x,0,1)) PyError ($(Expr(:escape, :(ccall(#= /Users/gangchen/.julia/packages/PyCall/ttONZ/src/pyfncall.jl:44 =# @pysym(:PyObject_Call), PyPtr, (PyPtr, PyPtr, PyPtr), o, pyargsptr, kw))))) <class 'sympy.integrals.meijerint._CoeffExpValueError'> _CoeffExpValueError('expr not of form ax**b: 0.857109621959463') File "/Users/gangchen/.julia/conda/3/lib/python3.7/site-packages/sympy/integrals/integrals.py", line 1477, in integrate return integral.doit(**doit_flags) File "/Users/gangchen/.julia/conda/3/lib/python3.7/site-packages/sympy/integrals/integrals.py", line 541, in doit function, xab[0], **eval_kwargs) File "/Users/gangchen/.julia/conda/3/lib/python3.7/site-packages/sympy/integrals/integrals.py", line 1012, in _eval_integral h = meijerint_indefinite(g, x) File "/Users/gangchen/.julia/conda/3/lib/python3.7/site-packages/sympy/integrals/meijerint.py", line 1621, in meijerint_indefinite res = _meijerint_indefinite_1(f.subs(x, x + a), x) File "/Users/gangchen/.julia/conda/3/lib/python3.7/site-packages/sympy/integrals/meijerint.py", line 1646, in _meijerint_indefinite_1 gs = _rewrite1(f, x) File "/Users/gangchen/.julia/conda/3/lib/python3.7/site-packages/sympy/integrals/meijerint.py", line 1569, in _rewrite1 g = _rewrite_single(g, x, recursive) File "/Users/gangchen/.julia/conda/3/lib/python3.7/site-packages/sympy/integrals/meijerint.py", line 1472, in _rewrite_single exponents_only=True), x) File "/Users/gangchen/.julia/conda/3/lib/python3.7/site-packages/sympy/integrals/meijerint.py", line 337, in _get_coeff_exp raise _CoeffExpValueError('expr not of form ax**b: %s' % expr)

Stacktrace: [1] pyerr_check at /Users/gangchen/.julia/packages/PyCall/ttONZ/src/exception.jl:60 [inlined] [2] pyerr_check at /Users/gangchen/.julia/packages/PyCall/ttONZ/src/exception.jl:64 [inlined] [3] macro expansion at /Users/gangchen/.julia/packages/PyCall/ttONZ/src/exception.jl:84 [inlined] [4] __pycall!(::PyCall.PyObject, ::Ptr{PyCall.PyObject_struct}, ::PyCall.PyObject, ::Ptr{Nothing}) at /Users/gangchen/.julia/packages/PyCall/ttONZ/src/pyfncall.jl:44 [5] _pycall!(::PyCall.PyObject, ::PyCall.PyObject, ::Tuple{Sym,Tuple{Sym,Int64,Int64}}, ::Int64, ::Ptr{Nothing}) at /Users/gangchen/.julia/packages/PyCall/ttONZ/src/pyfncall.jl:29 [6] _pycall!(::PyCall.PyObject, ::PyCall.PyObject, ::Tuple{Sym,Tuple{Sym,Int64,Int64}}, ::Base.Iterators.Pairs{Union{},Union{},Tuple{},NamedTuple{(),Tuple{}}}) at /Users/gangchen/.julia/packages/PyCall/ttONZ/src/pyfncall.jl:11 [7] #call#111(::Base.Iterators.Pairs{Union{},Union{},Tuple{},NamedTuple{(),Tuple{}}}, ::PyCall.PyObject, ::Sym, ::Vararg{Any,N} where N) at /Users/gangchen/.julia/packages/PyCall/ttONZ/src/pyfncall.jl:89 [8] (::PyCall.PyObject)(::Sym, ::Vararg{Any,N} where N) at /Users/gangchen/.julia/packages/PyCall/ttONZ/src/pyfncall.jl:89 [9] #integrate#212(::Base.Iterators.Pairs{Union{},Union{},Tuple{},NamedTuple{(),Tuple{}}}, ::Function, ::Sym, ::Tuple{Sym,Int64,Int64}) at /Users/gangchen/.julia/packages/SymPy/tgSHx/src/importexport.jl:102 [10] integrate(::Sym, ::Tuple{Sym,Int64,Int64}) at /Users/gangchen/.julia/packages/SymPy/tgSHx/src/importexport.jl:102 [11] top-level scope at In[4]:1

[AmplitudeGravity]

Although the int power function works.

integrate(x^(2)*(1+x)^(-3),(x,0,1)) −5/8+log(2)

[jverzani]

In general, you should use rational exponents to avoid floating point round off. However, in this case once doing so, the integral takes forever to compute. Even on sympygamma it fails. The integration with Wolfram Alpha is just better.