function F=pFep(z,q)
h=.01;
F=(pFe(z-2*h,q)-8*pFe(z-h,q)+8*pFe(z+h,q)-pFe(z+2*h,q))/(12*h);