Guest User

Untitled

a guest
Apr 21st, 2018
80
0
Never
Not a member of Pastebin yet? Sign Up, it unlocks many cool features!
text 0.72 KB | None | 0 0
  1. import numpy as np
  2. import fractions as fr
  3.  
  4. N = 4
  5. arr = np.empty((N * 2 + 1, N * 2 + 1, 2), dtype='object')
  6.  
  7. def plus(i, j):
  8. arr[i + 1, j] = 4 * arr[i, j] - arr[i - 1, j] - arr[i, j + 1] - arr[i, abs(j - 1)]
  9.  
  10. def even(i):
  11. arr[i, i] = arr[i - 1, i - 1] + [0, fr.Fraction(2, 2 * i - 1)]
  12. for k in range(1, i + 1): plus(i + k - 1, i - k)
  13.  
  14. def odd(i):
  15. arr[i + 1, i] = 2 * arr[i, i] - arr[i, i - 1]
  16. for k in range(1, i + 1): plus(i + k, i - k)
  17.  
  18. arr[0, 0] = 0
  19. arr[1, 0] = [fr.Fraction(1, 2), 0]
  20.  
  21. for i in range(1, N):
  22. even(i)
  23. odd(i)
  24.  
  25. even(N)
  26.  
  27. for i in range(0, N + 1):
  28. for j in range(0, N + 1):
  29. a, b = arr[max(i, j), min(i, j)]
  30. print('(', a, ')+(', b, ')/π', sep='', end='t')
  31. print()
Add Comment
Please, Sign In to add comment