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');
y=inv(eye(2)-a1)*b1
plot(y(1),y(2),'go');
y=inv(eye(2)-a2)*b2
plot(y(1),y(2),'ro');
y=inv(eye(2)-a3)*b3
plot(y(1),y(2),'bo');
y=inv(eye(2)-a4)*b4
plot(y(1),y(2),'mo');

