# -*- coding: utf-8 -*-
"""
Created on Mon Feb 08 09:56:21 2016

@author: hermanr
"""

from pylab import * 

N=10;
a=-2*pi;
b=2*pi;

t=linspace(a,b,100);
fs=pi/2;
for k in range(1,N): 
    fs=fs+(2/(pi*k**2))*((-1)**k-1)*cos(k*t);

plot(t,fs)
xlabel('x')
ylabel('f(x)')
title('Plot of Partial Sum' )
axis([a,b,0,pi])