%Initialize: c=0.38+0.34*i; %c=-0.5-0.56*i; z=0; h=[]; for j=1:10 z=z^2+c; end for j=1:15 z=z^2+c; h(j)=z; end plot(h)