function z=eval_pt(psi,f)
     f=(f+1)/2;
     v1=subs(psi,sym('p'),sym(f));
     %z=v1;
     dd=digits;
     z=vpa(v1,dd);