t = (0:0.01:31.415);
y = pi/2*cos(0*t);
axis ([0, 8, -0.5, 3.5]);  # x軸が0:32 , y軸が -0.5:3.5 の範囲を表示
for n = 1:10
	y = y - 4/pi/(2*n-1)/(2*n-1)*cos( (2*n-1)*t );
	plot (t, y );
	hold on;
endfor
print("hoge202.eps","-depsc2")