function g=eval_mpt(psi,x);
[r,c]=size(x);
g=sym(zeros(r,c));
for i=1:r
for j=1:c
g(i,j)=(eval_pt(psi,x(i,j)));
end
end