Advertisement
MatsGranvik

A155585 Bernoulli numbers matrix inverse sum

Dec 23rd, 2017
173
0
Never
Not a member of Pastebin yet? Sign Up, it unlocks many cool features!
text 0.42 KB | None | 0 0
  1. nn = 20;
  2. A = Table[
  3. Table[If[n >= k, Binomial[n - 1, k - 1], 0]*
  4. If[And[n > k, Mod[k, 2] == 0], 0, 1], {k, 1, nn}], {n, 1,
  5. nn}]; B =
  6. Table[Table[
  7. If[Mod[n - k, 2] == 1, 0, 1]*If[n >= k, Binomial[n - 1, k - 1], 0]*
  8. If[And[n > k, Mod[k, 2] == 0], 0, 1], {k, 1, nn}], {n, 1, nn}];
  9. W = Inverse[B] - Inverse[A] + IdentityMatrix[nn];
  10. Inverse[B][[All, 1]]
  11. Inverse[A][[All, 1]]
  12. %% - %
  13. https : // oeis.org/A155585
Advertisement
Add Comment
Please, Sign In to add comment
Advertisement