c33ae98Companion to 'Coding closed and open quantum systems in MATLAB' (arXiv:1911.04906)Jeremy Magland 1function Sci = getSci(sc,i,Nspins)
2Is = eye(2);
3Op_total = cell(1,Nspins);
4for site = 1:Nspins
5 Op_total{site} = Is+double(eq(i,site))*(sc-Is);
6end
7Sci = Op_total{1};
8for site = 2:Nspins
9 Sci = kron(Sci,Op_total{site});
10end
11end