2.4.1 MATLAB の場合

MATLAB の場合は前回の資料⁵で 説明してある程度の機能で十分のはずである。 行列のブロックに小さな行列を埋め込むことだけ説明しておく。

>> n=5
>> n2=2*n
>> at=(1:n)'*(1:n)
>> a=zeros(n2,n2);
>> a(1:n,1:n)=at;
>> a(n+1:n2,n+1:n2)=10*at

とすると、次のようになる。

a =

     1     2     3     4     5     0     0     0     0     0
     2     4     6     8    10     0     0     0     0     0
     3     6     9    12    15     0     0     0     0     0
     4     8    12    16    20     0     0     0     0     0
     5    10    15    20    25     0     0     0     0     0
     0     0     0     0     0    10    20    30    40    50
     0     0     0     0     0    20    40    60    80   100
     0     0     0     0     0    30    60    90   120   150
     0     0     0     0     0    40    80   120   160   200
     0     0     0     0     0    50   100   150   200   250