function z = stab_RK3(ar,ai)

for k=1:length(ai)
  a(k,:)=ar+i*ai(k);
end

z = abs(1+a+a.^2/2+a.^3/6+a.^4/24);