takes a list of operators and matrices and creates a block diagonal operator:</p><pre class="codeinput">A = opOnes(2,2);
B = 2*opOnes(3,2);
C = 3*opOnes(1,3);
D = blkdiag(A, B, C);
double(D)
</pre><pre class="codeoutput">
ans =

     1     1     0     0     0     0     0
     1     1     0     0     0     0     0
     0     0     2     2     0     0     0
     0     0     2     2     0     0     0
     0     0     2     2     0     0     0
     0     0     0     0     3     3     3

</pre><p>We can also hav