var theta = 90;

qc.reset(2);
qc.write(0);
qc.hadamard();

// Używamy dwóch CNOT i trzech PHASE…
qc.phase( theta / 2, 0x2);
qc.cnot(0x2, 0x1);
qc.phase(-theta / 2, 0x2);
qc.cnot(0x2, 0x1);
qc.phase( theta / 2, 0x1);

// Tworzymy tę samą operację jako dwukubitową CPHASE
qc.phase(theta, 0x1, 0x2);
