NDSolve[{x'[t] == 0.1 x[t] - 0.1 x[t] x[t - 10.2914]/40 - 0.02 x[t] z[t],
y'[t] == 0.03 x[t - 10.2914] y[t] - 0.02 y[t] - 0.05 y[t] y[t],
z'[t] == 0.02 x[t - 10.2914] z[t] - 0.02 y[t] z[t] - 0.04 z[t] -
0.006 z[t] z[t], x[0] == 10, y[0] == 4 z[0] == 2}, {x, y, z}, {t,
0, 600}, MaxSteps -> Infinity]
怎么把这种图形输出出来

可以查看InterpolatingFunction的帮助，其实可以把式中的xyz当成函数名，后面使用方括号表示函数的自变量就可以把InterpolatingFunction当成函数使用了，如下
NDSolve[{x'[t] ==
0.1 x[t] - 0.1 x[t] x[t - 10.2914]/40 - 0.02 x[t] z[t],
y'[t] == 0.03 x[t - 10.2914] y[t] - 0.02 y[t] - 0.05 y[t] y[t],
z'[t] ==
0.02 x[t - 10.2914] z[t] - 0.02 y[t] z[t] - 0.04 z[t] -
0.006 z[t] z[t], x[0] == 10, y[0] == 4 z[0] == 2}, {x, y, z}, {t,
0, 600}, MaxSteps -> Infinity]
{x, y, z} /. % // Flatten
Plot[#[x] & /@ %, {x, 0, 600}]