12.1.1 まず $A'$ を作ってから $A$ を作る

定義に素朴に従って、$A'$ を作るところから始めてみよう

$A'$ を MATLAB で作るには (似た例を既にプリントに出してあるけれど)、 例えば

m=50;
Ap=4*eye(m,m)-diag(ones(m-1,1),1)-diag(ones(m-1,1),-1);

とすれば OK. それをもとに $A$ を作るにはどうするか？

B から $A$ を作る方法 (1)

c1=1; c2=10
A=[c1*Ap zeros(m,m); zeros(m,m) c2*Ap];

または

B から $A$ を作る方法 (2)

c1=1; c2=10;
n=2*m;
A=zeros(n,n);
A(1:m,1:m)=c1*Ap;
A(m+1:n,m+1:n)=c2*Ap;