| New Reply |
Numeric integration in python |
Share Thread | Thread Tools |
| Nov14-12, 02:24 PM | #1 |
|
|
Numeric integration in python
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? |
| Nov15-12, 06:38 AM | #2 |
|
Recognitions:
|
[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] |
| Nov15-12, 06:55 AM | #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. |
| Nov15-12, 07:33 AM | #4 |
|
Recognitions:
|
Numeric integration in pythonQuad 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
|
| Nov15-12, 09:14 AM | #5 |
|
|
Your explanation makes lots of sense...
I'd defined a function but had difficulty making it a function of t.... Thanks. |
| New Reply |
| Tags |
| numeric integration, python, scipy |
| Thread Tools | |
Similar Threads for: Numeric integration in python
|
||||
| Thread | Forum | Replies | ||
| Compile Python, Matlab and Python | Programming & Comp Sci | 3 | ||
| Numeric Integration | Engineering, Comp Sci, & Technology Homework | 0 | ||
| Many-Body Numeric Integration Algorithm | Differential Equations | 1 | ||
| Numeric keypad | Computers | 0 | ||
| Numeric Nonsense | Brain Teasers | 2 | ||