function y=f3(x,m)
global r
if m==0
  p=1;
else
  p=ones(size(x,1),1);
  for k=1:m/2
    p=(x.^2+(2*r*k*pi)^2).*p;
  end
end
y=p.*(einstein(x/r)).^r;
