function[xx,yy] = przyk1(x,y);
k=grand(1,1,'uin',1,6);
    if k==1
        xx=0.195*x-0.488*y+0.4431;
        yy=0.344*x+0.443*y+0.2452;
    end;
    if k==2
        xx=0.462*x+0.414*y+0.2511;
        yy=-0.252*x+0.361*y+0.5692;
    end;
    if k==3;
        xx=-0.058*x-0.070*y+0.5976;
        yy=0.453*x-0.111*y+0.0969;
    end;
    if k==4
        xx=-0.035*x+0.070*y+0.4884;
        yy=-0.469*x-0.022*y+0.5069;
    end;
    if k==5 then
        xx = -0.637*x+0.8532;
        yy = 0.501*y+0.2513;
    end
    if k==6 then
        xx = 0.56*x+0.56*y+10.0;
        yy = 0.44*x+0.44*y+10.0;
    end
endfunction;

x=0.6;
y=0.7;
for i=1:20000
[xx,yy] = przyk1(x,y);
x=xx;
y=yy;
xx0(i)=x;
yy0(i)=y;
end;


f=gdf();
f.background=12;
xset('background',12);
xset('color',7);
subplot(2,2,1)
//plot2d(xx1,yy1,style=0);
plot(xx0,yy0,"y.x");

//*********************************************

function[xx,yy] = przyk2(x,y);
k=grand(1,1,'uin',1,5);
    if k==1
        xx=0.195*x-0.488*y+0.4431;
        yy=0.344*x+0.443*y+0.2452;
    end;
    if k==2
        xx=0.462*x+0.414*y+0.2511;
        yy=-0.252*x+0.361*y+0.5692;
    end;
    if k==3;
        xx=-0.058*x-0.070*y+0.5976;
        yy=0.453*x-0.111*y+0.0969;
    end;
    if k==4
        xx=-0.035*x+0.070*y+0.4884;
        yy=-0.469*x-0.022*y+0.5069;
    end;
    if k==5 then
        xx = -0.5*x+0.5;
        yy = 0.5*y+0.5;
    end
endfunction;

x=0.6;
y=0.7;
for i=1:20000
[xx,yy] = przyk2(x,y);
x=xx;
y=yy;
xx1(i)=x;
yy1(i)=y;
end;


//g=gdf();
//g.background=2;
//xset('window',1)
xset('background',14);
xset('color',7);
subplot(2,2,2)
//plot2d(xx1,yy1,style=0);
plot(xx1,yy1,"b.x");

//*********************************************

function[xx,yy] = przyk3(x,y);
    k = grand(1, 1, 'uin', 1, 6);
    if k==1 then
        xx = -0.67*x+0.1*y-0.0;
        yy = -0.18*x+0.81*y-10.0;
    end;
    if k==2 then
        xx = 0.4*x-0.4*y+0.0;
        yy = -0.1*x-0.4*y+0.0;
    end;
    if k==3 then
        xx = -0.4*x-0.4*y-0.0;
        yy = -0.1*x+0.4*y;
    end;
    if k==4 then
        xx = -0.1*x+2.0*y;
        yy = -0.44*x + 0.44*y-2.0;
    end;
    if k==5 then
        xx = -0.5252*x + 0.5252*y ;
        yy = 0.55*x - 0.55*y ;
    end;
    if k==6 then
        xx= -0.601*x - 0.512*y+0.0;
        yy= -0.401*x + 0.712*y+0.0;
    end
endfunction

x=0.6;
y=0.7;

for i=1:20000
[xx,yy] = przyk3(x,y);
x=xx;
y=yy;
xx3(i)=x;
yy3(i)=y;
end;
//xset('window',2)
//clf();
xset('background',3);
subplot(2,2,3);
//plot2d(xx1,yy1,style=0);
plot(xx3,yy3,"r.*");

//**********************************************

function [xx,yy]=przyk4(x,y);
    k=grand(1,1,'uin',1,2);
    if k==1
      xx=0.824*x+0.281*y-1.88;
      yy=-0.212*x+0.86*y-0.11;
    end;
    if k==2
      xx=0.088*x+0.521*y+0.78;
      yy=-0.463*x-0.377*y+8.09;
    end;
endfunction;

x=0.6;
  y=0.7;

for i=1:20000
  [xx,yy] = przyk4(x,y);
  x=xx;
  y=yy;
  xx4(i)=x;
  yy4(i)=y;
end;
subplot(2,2,4);
plot(xx4,yy4,"o.x");
xset('background',5);