function e=exyz(p)
% function e=exyz(p)  calls function q=xyz(p), then sets e=q'*q
q=xyz(p);
e=q'*q;
