(*start*)
nn = 40;
T[n_, k_] :=
T[n, k] =
Which[n < 1 || k < 1, 0, n == 1 || k == 1, 1, k > n, T[k, n], n > k,
T[k, Mod[n, k, 1]],
True, -Sum[T[n, i], {i, n - 1}]];(*Michael Somos,Jul 18 2011*)
A = Chop[Table[
Table[N[T[n, k]/n^ZetaZero[1]], {k, 1, nn}], {n, 1, nn}]];
Total[MatrixExp[A t].(Range[nn]*0 + 1)];
Plot[Sign[Re[%]]*Log[Abs[Re[%]]], {t, -42, 42}]
(*end*)