function fern3(num)
a1=[.85 .04;-0.04 .85];
b1=[0;1.6];
a2=[0.2 -0.26;0.23 0.22];
b2=[0;1.6];
a3=[-0.15 0.28; 0.26 0.24];
b3=[0;0.44];
a4=[0 0; 0 0.16];
b4=[0 ;0];
p=[0.85; 0.92; 0.99; 1];
x=[1;1];
s(1,1)=x(1);
s(1,2)=x(2);
figure;
axis([-6,6,-1,11]);
hold on
plot(s(1,1),s(1,2),'ko');
for i=2:num
x=a1*x+b1;
s(i,1)=x(1);
s(i,2)=x(2);
plot(s(i,1),s(i,2),'go');
pause(0.5)
end
x=[1;1];
for i=2:num
x=a2*x+b2;
s(i,1)=x(1);
s(i,2)=x(2);
plot(s(i,1),s(i,2),'ro');
pause(0.5)
end
x=[1;1];
for i=2:num
x=a3*x+b3;
s(i,1)=x(1);
s(i,2)=x(2);
plot(s(i,1),s(i,2),'bo');
pause(0.5)
end
x=[1;1];
for i=2:num
x=a4*x+b4;
s(i,1)=x(1);
s(i,2)=x(2);
plot(s(i,1),s(i,2),'mo');
pause(0.5)
end

end

