Matplot1

clf();

ax=gca();//カレントの軸のハンドルを取得

ax.data_bounds=[0,0;10,10];//data_boundsを設定

ax.box='on'; //ボックスを描画 

a=5*ones(11,11); a(2:10,2:10)=4; a(5:7,5:7)=2;

// 矩形 [1,1,3,3]に最初の行列をプロット

Matplot1(a,[1,1,3,3])

a=ones(10,10); a= 3*tril(a)+ 2*a; 

// 矩形 [5,6,7,8]に二番目の行列をプロット

Matplot1(a,[5,6,7,8])

//--- アニメーション)

n=100;

clf();

f = gcf();//カレントの図のハンドルを取得

ax = gca();//カレントの軸のハンドルを取得

ax.data_bounds = [0, 0 ; 10, 10];//data_boundsを設定

ax.box = 'on'; //ボックスを描画 

o = ones(n, n);

a = 3 * tril(o, -n) + 2*o;

a = a + a';

Matplot1(a, [0, 2, 7, 9])

mp = gce();

for k = (-n + 1):n,

  a = 3 * tril(o, k) + 2 * o;

  a = a + a';

  k1 = 3 * (k + 100) / 200;

  mp.data = a;

  mp.rect = [k1, 2, k1 + 7, 9];

  sleep(8); // シミュレーションをスローダウン

end