On Mar 26, 2009, at 6:58 PM, [email protected] wrote:

> This will bomb because numerical_integral doesn't allow vars for the
> limits of integration:
>                      Y(t)=2500+numerical_integral(S(u)-R(u),0,t)[0]
>
> I tried Y(t)=2500+numerical_integral(S(u)-R(u),0,x).substitute(x=t)[0]
>
> and
>
> Y(t) = ( 2500+numerical_integral(S(u)-R(u),0,x)[0]).substitute(x=t)
>
> but no luck

You can do

def Y(t):
     return 2500+numerical_integral(S(u)-R(u),0,t)[0]

but then it won't be a symbolic object. (It will be a Python function.)

- Robert


--~--~---------~--~----~------------~-------~--~----~
To post to this group, send email to [email protected]
To unsubscribe from this group, send email to 
[email protected]
For more options, visit this group at 
http://groups.google.com/group/sage-support
URLs: http://www.sagemath.org
-~----------~----~----~----~------~----~------~--~---

Reply via email to