Recursive sequences in CAS on prime.

Recursive sequences in CAS on prime. - Printable Version

+- HP Forums (https://www.hpmuseum.org/forum)
+-- Forum: HP Calculators (and very old HP Computers) (/forum-3.html)
+--- Forum: HP Prime (/forum-5.html)
+--- Thread: Recursive sequences in CAS on prime. (/thread-4185.html)

Recursive sequences in CAS on prime. - wawa - 06-19-2015 07:22 AM

Hello

I can't figure out how to define a recursive sequence in CAS on the HP prime.
I want to define:
u(1)=3/2
u(2)=5/3
u(n)=2003-6002/u(n-1)+4000/(u(n-1)u(n-2))

Of course, I can define this sequence in the sequence app in HOME mode but it gives wrong answers (the sequence is made for that purpose Smile

) so I want to have exact values in CAS.

I made a CAS program with a FOR loop and it works fine but I'm looking for a more straightforward way to do this.
In XCAS, I can type in

Code:

u(n):={si n==1 alors 3/2; sinon si n==2 alors 5/3; sinon 2003-6002/u(n-1)+4000/(u(n-1)*u(n-2)); fsi; fsi;}

and it works.
I am looking for a similar syntax working on the prime too.

thank you for your help.

RE: Recursive sequences in CAS on prime. - Arno K - 06-19-2015 09:01 AM

Here is a solution:

Code:

#cas

mf(n):=

BEGIN

IF n==1 THEN

RETURN(3/2);

ELSE

IF n==2 THEN

RETURN(5/3);

END;

RETURN(2003-6002/mf(n-1)+4000/(mf(n-1)*mf(n-2)));

END;

END;

#end

which seems to work in your desired way. It provides exact values at least until n=10, now I am trying 25 and this runs now the fifth minute on the emulator which is due to the 3 recursive calls in the code, so I cancelled it.
Arno

RE: Recursive sequences in CAS on prime. - wawa - 06-19-2015 09:20 AM

Thank you. This is a way to do this. But is there another way, without making a CAS program ? with a declaration like u(n):= ... ?

RE: Recursive sequences in CAS on prime. - ww63 - 06-19-2015 10:33 AM

Hello Arne,

maybe your program runs faster, if you create two local variables, wich get the values of mf(n-1) and mf(n-2), so you have to do the recursive call only twice ind not three times.

regards

Wolfgang

RE: Recursive sequences in CAS on prime. - Arno K - 06-19-2015 02:12 PM

(06-19-2015 10:33 AM)ww63 Wrote: Hello Arne,

maybe your program runs faster, if you create two local variables, wich get the values of mf(n-1) and mf(n-2), so you have to do the recursive call only twice ind not three times.

regards

Wolfgang

Yes, you are right, a reduction of recursive calls is necessary. But your solution, I think, will not work with local variables as they are created each time when the function is called. Here additional Parameters are necessary to pass the last two values and then the solution with a for-loop seems better.
Arno

RE: Recursive sequences in CAS on prime. - wawa - 06-19-2015 04:37 PM

I managed to declare the sequence directly in CAS:

Code:

u:=(n)->CASE IF (n==1) THEN 3/2 END IF (n==2) THEN 5/3 END IF n>2 THEN 2003-6002/u(n-1)+4000/(u(n-1)*u(n-2)) END END

It works fine until n>10. After this, the recursion is too deep.

I made a CAS program too with a FOR loop :

Code:

#cas

u(n):=

BEGIN

LOCAL r,r1,r2

CASE

IF n==1 THEN RETURN 3/2; END;

IF n==2 THEN RETURN 5/3; END;

DEFAULT

r1=3/2

r2=5/3

FOR I FROM 3 TO n DO

r := 2003-6002/r2+4000/(r1*r2);

r1:=r2;

r2:=r;

END;

RETURN r;

END;

#end

But after all this, I think my favorite method is to use the spreadsheet with cells in CAS mode. It is fast and quite easy.