Advertisement
Not a member of Pastebin yet?
Sign Up,
it unlocks many cool features!
- ilist[l_] := Range[Length[l]];
- applyN[f_, a_] := f /@ Zip[ilist[a], a];
- uArrTest[n_, m_] := Module[{lAm, arr1, arr2, arr3},
- lAm = ConstantArray[1, m];
- arr1 = NestList[ListConvolve[lAm, #, {1, -1}, 0] &, lAm, n - 1];
- arr2 = applyN[PadRight[ ArrayPad[#[[2]], {#[[1]] - 1, 0}], n m] &,
- arr1];
- arr3 = Accumulate /@ arr2;
- N[# / #[[-1]]] & /@ arr3
- ];
- For[m = 56, m <= 100, m++, u100tmp = uArrTest[100, m];
- Export["~/un100m" <> ToString[m] <>
- ".dat", u100tmp]; Print[m]]
Advertisement
Add Comment
Please, Sign In to add comment
Advertisement