1function summary = mtrace_summary(T,C)
4summary.M = [mean(T), mean(T.^2), mean(T.^3), mean(T.^4), mean(T.^5)];
5summary.ACF = trace_acf(T, 1:100);
6summary.F1 = mtrace_forward_moment(T, C, 1);
7summary.F2 = mtrace_forward_moment(T, C, 2);
8summary.B1 = mtrace_backward_moment(T, C, 1);
9summary.B2 = mtrace_backward_moment(T, C, 2);
10summary.C1 = mtrace_cross_moment(T, C, 1);
11summary.C2 = mtrace_cross_moment(T, C, 2);
12summary.Pc = mtrace_pc(T, C);
13summary.Pab = mtrace_sigma(T, C);