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); // Slow down the simulation
end