41Z Routine of the Week: Complex Fibonacci

[Ángel Martin]

Very nice JM, thanks for your contribution.

here's the 41Z version to use for complex orders - which requires a modified Binet formula using an interesting correction term:

input is the complex order z in {Y, X}
41 bytes / SIZE 000

Code:

01 LBL "ZFIB"
02  ZRPL^
03  PHI
04  X^Z   
05  ZENTER^
06  ZINV
07  ZRUP
08  PI
09  ST*  Z
10  *
11  ZCOS
12  Z*
13  Z-
14  5
15  SQRT
16  ST/  Z
17  /
18  ZAVIEW
19  END

Example: Fi(3+4i) ~ -5,248.511-J14,195.962

Cheers,
ÁM