[X Y] = meshgrid(linspace(-5,5),linspace(-5,5));
Z=X+i*Y;

figure(1);

for s=1:5
    R = 1+Z.^s/factorial(s);
    contour(X,Y,abs(R)); %subplot(5,1,s);
    pause
end