syms p
psi=cos(2*pi*(p*p-p-1/16))/cos(2*pi*p);
c0=psi;

m1=[ 1 -96 2 3];
m2=[ 1 18432 4 6; 1 64 2 2];
m3=[1 -5308416 6 9; 1 -3840 4 5; 1 -64 2 1];
m4=[1 2038431744 8 12; 11 5898240 6 8; 19 24576 4 4; 1 128 2 0];
m5=[-5 3072 4 3; -901 82575360 6 7; -7 849346560 8 11; -1 978447237120 10 15];
c1= sum_many(m1);
c2= sum_many(m2);
c3= sum_many(m3);
c4= sum_many(m4);
c5= sum_many(m5);