function z= eval_a_term(fracnum,fracden,pi_pow,deriv)
syms p
psi=cos(2*pi*(p*p-p-1/16))/cos(2*pi*p);
z=diff(psi,p,deriv)*fracnum/(fracden*sym(pi)^pi_pow);