function morse=morse_erk(e0,r0,k,ref=0)
  a=sqrt(k/2/e0);
  morse=@(x) morse_aux(x,e0,r0,a)+ref;
end

function res=morse_aux(x,e0,r0,a)
  expp=exp(-a*(x-r0));
  res=e0*expp.*(expp-2);
end
