x=-pi/2:.01:pi;
y=a1f(x);
plot(x,y)
xlabel('x')
ylabel('f(x)')