function y=f5(x,m)
if m==0
  p=1;
else
%  p=ones(size(x,1),1);
  p=(x.^2+(2*m*pi)^2).^m;
end
y=p.*(einstein(x/m)).^m;
