a=linspace(-pi,pi,100);
plot(a,f(a))
