NumPy: Fibonacci-like

02-numpy/fib-solve.py

+#!/usr/bin/python3
+
+import numpy as np
+
+n = 10
+
+a = -np.eye(n)/2 - np.eye(n, k=1)/2 + np.eye(n, k=2)
+a[-2:, -2:] = np.eye(2)
+
+b = np.zeros(n)
+b[-2] = 426
+b[-1] = 427
+
+x = np.linalg.solve(a, b)
+print(x[:2])

02-numpy/fib-test.py

+#!/usr/bin/python3
+
+a = [256, 512]
+for i in range(2, 10):
+    a.append((a[-1] + a[-2]) / 2)
+print(a)