Numeric integration in python

  • Python
  • Thread starter sola maths
  • Start date
  • #1

Main Question or Discussion Point

I'm trying to write a python program that is able to numerically execute functions of the form:

y(t) = exp(Integrate[A(x),x]) within the bounds of 0 and t

I tried using quad from scipy.integrate but it seems not to be able to evaluate expressions of this form.

Any other suggestions on appropriate packages or commands?
 

Answers and Replies

  • #2
uart
Science Advisor
2,776
9
I'm trying to write a python program that is able to numerically execute functions of the form:

y(t) = exp(Integrate[A(x),x]) within the bounds of 0 and t

I tried using quad from scipy.integrate but it seems not to be able to evaluate expressions of this form.

Any other suggestions on appropriate packages or commands?
Can you be a bit more specific sola maths? Do you mean:

[tex]y(t) = \exp(\int_0^t A(x) \, dx)[/tex]

You have some specific function A(x) that you can evaluate numerically?

Also, is that just you integrand? Do you need to evaluate,

[tex]\int_a^b \, e^{\int_0^t A(x) \, dx} \, dt[/tex]
 
  • #3
Hi uart,

Yes, I have some function A(x) that I need to evaluate numerically as x changes. The first expression you wrote is what I meant.
 
  • #4
uart
Science Advisor
2,776
9
Hi uart,

Yes, I have some function A(x) that I need to evaluate numerically as x changes. The first expression you wrote is what I meant.
Ok so it's just [itex]y(t) = \exp(\int_0^t A(x) \, dx)[/itex] that you need to evaluate.

Quad can do this easily, but only for one particular value of "t" at a time. However you could call it (quad) from within a function if you wished to properly make a function of "t". For example,

Code:
def functA(x):
     return x*x/2.0

def functY(t):
     return exp(integrate.quad(functA,0,t)[0])

functY(3)
90.017131300521896
 
  • #5
Your explanation makes lots of sense...

I'd defined a function but had difficulty making it a function of t.... Thanks.
 

Related Threads for: Numeric integration in python

Replies
50
Views
2K
  • Last Post
Replies
3
Views
490
  • Last Post
Replies
8
Views
3K
Replies
21
Views
964
Replies
6
Views
6K
Replies
6
Views
2K
  • Last Post
Replies
6
Views
887
Replies
5
Views
658
Top