Advertisement
Guest User

Untitled

a guest
Jul 25th, 2014
232
0
Never
Not a member of Pastebin yet? Sign Up, it unlocks many cool features!
text 0.83 KB | None | 0 0
  1. Clear[t, n, k, a, b, x, y, b1, at, T, AT, at];
  2. nn = 200;
  3. x = 2 + 3 I;
  4. y = 5 + 7*I;
  5. T[n_, 1] = 1;
  6. T[n_, k_] :=
  7. T[n, k] =
  8. If[n >= k,
  9. x*Sum[T[n - i, k - 1], {i, 1, n - 1}] +
  10. y*Sum[T[n - i, k], {i, 1, n - 1}], 0];
  11. a = Table[Table[Expand[T[n, k]], {k, 1, nn}], {n, 1, nn}];
  12. (*TableForm[a]*)
  13. AT = Total[Transpose[a]];
  14. AT[[1 ;; 20]];
  15. Table[Limit[((s + 1)^(n - 1) + s - 1)/s, s -> y + 1], {n, 1, 12}];
  16. c = Product[1/(1 - 1/AT[[i]]), {i, 2, nn}];
  17. N[c, 12]
  18.  
  19. Clear[t];
  20. t[n_, 1] = 1;
  21. t[n_, k_] :=
  22. t[n, k] =
  23. If[n >= k,
  24. x*Sum[t[n - i, k - 1], {i, 1, k - 1}] +
  25. y*Sum[t[n - i, k], {i, 1, k - 1}], 0];
  26. a = Table[Table[Expand[t[n, k]], {k, 1, nn}], {n, 1, nn}];
  27. aa = Table[
  28. Total[Table[t[n - k + 1, k], {k, 1, nn}]], {n, 1 + 1, nn + 1}];
  29. (*TableForm[a]*)
  30. at = Total[Transpose[a]];
  31. at[[1 ;; 20]];
  32. N[AT[[nn]]/at[[nn]], 12]
Advertisement
Add Comment
Please, Sign In to add comment
Advertisement