Wolfram / 方形波の合成

s[x_, n_] = (4/Pi)*Sum[Sin[(2 k - 1) x]/(2 k - 1), {k, 1, n}];
Manipulate[Plot[Evaluate[s[x, n]], {x, -2 Pi, 2 Pi}, PlotRange -> {-3, 3}], {n, 1, 10, 1}]

参考: レクチャーズオンMathematica, pp.91-92