Advertisement
Guest User

Untitled

a guest
Jan 30th, 2015
172
0
Never
Not a member of Pastebin yet? Sign Up, it unlocks many cool features!
text 0.49 KB | None | 0 0
  1. ilist[l_] := Range[Length[l]];
  2. applyN[f_, a_] := f /@ Zip[ilist[a], a];
  3. uArrTest[n_, m_] := Module[{lAm, arr1, arr2, arr3},
  4. lAm = ConstantArray[1, m];
  5. arr1 = NestList[ListConvolve[lAm, #, {1, -1}, 0] &, lAm, n - 1];
  6. arr2 = applyN[PadRight[ ArrayPad[#[[2]], {#[[1]] - 1, 0}], n m] &,
  7. arr1];
  8. arr3 = Accumulate /@ arr2;
  9. N[# / #[[-1]]] & /@ arr3
  10. ];
  11. For[m = 56, m <= 100, m++, u100tmp = uArrTest[100, m];
  12. Export["~/un100m" <> ToString[m] <>
  13. ".dat", u100tmp]; Print[m]]
Advertisement
Add Comment
Please, Sign In to add comment
Advertisement