That works! And it's faster than our implementation, do you know why?
I think its the datatype of the number h cause type(1/10) is <type
'sage.rings.rational.Rational'>
Anyway thanks again.
On 2 oct, 18:58, "David Joyner" <[EMAIL PROTECTED]> wrote:
> This works:
>
> var('t')
> h = 1/10
> f = exp(t)*cos(t) - exp(t)*sin(t) #la ecuacion diferencial
>
> def x(n):
> if n == 0:
> return 1 #punto conocido
> else:
> return (x(n-h) + h*f(n-h)) #funcion de Euler
>
> for i in range(10):
> print RR(x(i*h))
>
> On Thu, Oct 2, 2008 at 12:35 PM, Janzo <[EMAIL PROTECTED]> wrote:
>
> > Hi guys I have a problem with this recursive function, I am getting a
> > infinite loop.
>
> > The thing is that it is working fine for x(h), x(2*h) but not for any
> > upper value.
>
> > Any suggestions? Thanks a lot.
>
> > var('h,t')
> > h = 0.1
> > f = exp(t)*cos(t) - exp(t)*sin(t) #la ecuacion diferencial
>
> > def x(n):
> > if n == 0:
> > return 1 #punto conocido
> > else:
> > return (x(n-h) + h*f(n-h)) #funcion de Euler
>
> > for i in range(10):
> > print x(i*h)
--~--~---------~--~----~------------~-------~--~----~
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
-~----------~----~----~----~------~----~------~--~---
